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Abstract: This thesis is divided in two parts, each one addressing problems that can be 
relevant in the study of compact objects. The first part deals with the study of a magnetized 
and self-gravitating gas of degenerated fermions (electrons and neutrons) as sources of a 
Bianchi-I space-time. We solve numerically the Einstein-Maxwell field equations for a large 
set of initial conditions of the dynamical variables. The collapsing singularity is isotropic 
for the neutron gas and can be anisotropic for the electron gas. This result is consistent 
with the fact thatelectrons exhibit a stronger coupling with the magnetic field, which is 
the source of anisotropy in the dynamical variables. In the second part we calculate the 
entropy of extremal black holes in 4 and 5 dimensions, using the entropy function formalism 
of Sen and taking into account higher order derivative terms that come from the complete 
set of Riemann invariants. The resulting entropies show the deviations from the well know 
Bekenstein-Hawking area law. 
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Intro duccion 



Los objetos compactos: Enanas Blancas (EBs), Estrellas Neutronicas (ENs) y Agujeros 
Negros (ANs) son objetos muy densos que invaden el Universo y son el ultimo estadio por 
el que transita una estrella antes de apagarse. Que una cstrella se convierta en uno u otro 
tipo de objeto compacto dcpcndc decisivamentc de su masa. Por ejemplo para estrellas 
con una masa entre 9 ~ IOMq, que lian agotado su combustible nuclear evolucionan hacia 
una EB. Sin embargo, estrellas con masa mayor que IOMq forman ANs. 

En particular las EBs y ENs son laboratorios astrofisicos naturales en los cuales se 
puede probar la veracidad de las teorias fisicas mas actuales ya que en ellas viven fenomenos 
en los que estan presentes parti'culas fuertemente interactuando, en presencia de grandes 
campos gravitacionalcs y clectromagneticos. 

El modelo mas aceptado para la descripcion de parti'culas en el micro mundo es el Mod- 
elo Estandar (ME). Por otra parte, la Teorfa General de la Relatividad (TGR) describe 
fielmente los fenomenos que ocurren a una mayor escala. Por tanto, un modelo fi'sico que 
pretenda describir estrellas de este tipo, debe incluir la gravedad y las interacciones entre 
particulas para que sea valido estrictamente hablando. Sin embargo en nuestro cuadro fi'sico 
del Universo, todavfa no tenemos una teoria unificada y experimentalmente demostrada 
que pueda juntar las cuatro fuerzas fundamentales de la naturaleza (fuerte, debil, electro- 
magnetica y gravitatoria) de modo que se puedan realizar predicciones teoricas medibles 
experimentalmente. De las fuerzas mencionadas, la linica que se niega a ser unificada 
es la gravitatoria, aunque en cstos liltimos tiempos se viene trabajando fuertemente en 
teori'as prometedoras (como las "teori'as de cuerdas") [1]. Las Teori'as de Cuerdas no tienen 
por aliora muchas esperanzas de ser probadas experimentalmente debido a que se requiere 
para verificarlas, de experimentos que produzcan energi'as superiores a las que se espera, 
se generen en el LHC^. El LHC fue disenado y puesto en marcha para encontrar el boson 
de Higgs que es un pilar fundamental del ME [2] . Sin embargo, existen casos dentro de la 
misma teoria de cuerdas en que el espacio de configuracion es tan grande que las cantidades 
de energfa para su deteccion no es tan alta, si este fuera el caso real, entonces con el LHC 
se podn'an ver indicios de estos objetos. 

Por otro lado, desde hace varias decadas se le da gran importancia al estudio de 
campos magneticos intensos en el universo y en particular en objetos compactos. En los 
liltimos tiempos, con la aparicion de nuevas y mas precisas observaciones, se avanza en las 
mediciones astroffsicas y se introducen nuevas cotas para estos campos magneticos. Asi' en 
los alios 80 campos magneticos dc 10^ Tesla (T) ^ eran considerados maximos. Hoy en dia, 
segiin observaciones mas refinadas se estiman que existen campos 1000 veces mayores que 
el conocido campo limite de Schwinger 4,4 x 10^ T [3]. 

Las EBs son estrellas con una masa de 1,2Mq, (o sea < 1,44M0), densidades entre 
10^ — radios de 10^ m y campos magneticos del orden de 10'* T. 

Sin embargo, una Estrella Neutronica tiene una masa M entre 1.5 y 2 masas solares 
Mq, radios entre 10 — 12 km y una densidad central de 10^^ kg /m^. Es por tanto una 

^LHC significa "Large Hadron Collider" o sea el Gran Colisionador de Hadrones. 
^ITesla = W^Gauss conversion para la induccion del campo magnetico. 
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de las formas mas densas dc materia encontradas en el Universe observable. Aunque los 
neutroncs dominan la componente barionica de la Estrella de Neutrones, existen tambien 
algunos protones (y suficientes electrones y muones) que garantizan la neutralidad de carga 
en la estrella. Ademas, engendran un super campo magnetico de unos 10^° T [4]. 

La primera EB encontrada, fue la companera de la brillante estrella Sirio. Sirio y su 
companera estan en una orbita mutua, una alrcdcdor dc la otra, y esto pcrmitio que se 
determinaran las masas dc cada una. A partir del brillo y temperatura de la compafiera, 
podemos determinar su tamano, que es de unos 10^ m de diametro, menos que el de la 
Tierra [5]. 

Las EBs son intrmsecamente muy poco brillantes, y son por tanto difi'ciles de detectar. 
Aun asi, ellas son el estado final de todas las estrellas de masa mediana, y por ello podrfamos 
esperar hallar muchas EBs en el Universo. 

Los astronomos han logrado encontrar muchas, usando tecnicas que dependen ya sea 
de que son companeras de otras estrellas, o dc que son estrellas calientes con grandes 
movimientos relatives respecto a otras estrellas (indicando que ellas estan mucho mas 
cercanas que otras estrellas de la Secuencias Principal con la misma temperatura), y a 
partir de su emision de radiacion de alta energfa, tal como luz ultravioleta. La presencia 
de EBs en sistemas binaries ha sido muy importante para entender violentos estallidos en 
sistemas estelares. Supernovas del Tipo I, Novas, y estrellas variables catacKsmicas, son 
todos casos en los que la compafiera dc una EB ha alcanzado un punto en su evolucion en 
el que esta aumentando cn diametro y pcrdicndo masa hacia la EB. 

La deposicion de material en un disco de acrecion alrededor de la EB, o en la superficie 
de la EB, determinara la naturaleza de cualquier estallido. La estrella aparentemente mas 
brillante en el cielo, Sirio, fue observada por Bessel en 1844 mostrando un bamboleo en su 
movimiento a traves del cielo. Bessel atribuyo esto a la presencia de una companera, pero 
no se observe alguna hasta que Alvan Clark, mientras prebaba un nuevo telescopio, vie la 
tenue estrella compafiera. En 1925, el espectre dc la estrella compafiera confirmo que era 
una estrella con aproximadamente la misma temperatura que Sirio A. 

Ella tiene un perfodo de 50 afios, con una maxima separacion en el cielo de 7,6 arcose- 
gundos. La diferencia de luminosidad entre Sirio A y B llega a un factor de mas de 8 000. 
La solucion de su movimiento orbital arrojo para A y B, masas de 2,3 y 1 veces la masa 
del Sol. Sirio A tiene un radio de cerca de 10^ m, mientras que Sirio B tiene un radio de 
solo 10^ m. 

En 1926 R. H. Fowler junto con Paul Dirac trabajaron en aplicar la mccanica cstadfstica 
para explicar la estabilidad de una estrella EB y en 1930 Chandrasekhar descubre la masa 
maxima de 1.44 para estas estrellas [6, 7]. 

Para las estrellas neutronicas los resultados teoricos antecedieron a los observacionales. 
En 1932 Landau hablo por primera vez de la EN. En el afio 1934 se predijo por primera 
vez, basandose en calculos teoricos, la existencia de una ENs y a partir de esa fecha se 
estudiaron muchos modelos para dcscribirlas, pero no es hasta el aho 1967 en que se verifica 
observacionalmente por azar, la existencia de un objeto con las caracterfsticas de una EN 
[8]. 
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Las ENs no se habi'an observado antes debido a que no emiten en el visible, por lo cual 
fue necesario que pasaran mas de 30 aiios y que la tecnologi'a permitiera la construccion 
de radiotelescopios capaces de detectar las primeras ENs. El descubrimiento de la primera 
EN provoco un aumento en el niimero de trabajos teoricos encaminados a describir los 
procesos fi'sicos que ocurren en el interior de estos objetos [9, 10, 11, 12, 13]. 

Desde hace veinte anos los estudios de las Estrellas de Neutroncs se nutren de cuantiosas 
observaciones que sc obtienen gracias a observatories que detectan emisiones de rayos X 
y radiacion gamma colocados en satclites. Entre los observatorios mas conocidos estan 
el Telescopio Espacial Hubble, el Observatorio de Rayos X Chandra y el Observatorio de 
Rayos Gama Compton. 

Actualmente tambien se habla de la posible existencia de Estrellas de Quark, o Es- 
trellas Extranas. Estas deben aparecer para super altas densidades, donde la materia 
presenta caracten'sticas todavia mas exoticas, tales como bariones con extraiieza, mesones 
condensados (piones, kaones), o quarks no confinados [14, 15, 16]. 

Estos fcrmiones en forma de hadrones o quark no confinados se espera que tambien 
exhiban propiedades de superfluidez o superconductividad. 

Una EN normal, tiene materia hadronica en su exterior y tanto la presion como la 
densidad de materia barionica se anula en la frontera exterior (en el interior puede con- 
tener cualquier combinacion de particulas exoticas permitidas por la fisica de interacciones 
fuertes). Sin embargo, una Estrella de quarks, puede tener tambien una superficie de quark 
desnuda cuya presion sc desvanece cn la frontera y una siiper alta densidad dc materia, o 
una fina capa de materia normal soportada por las fuerzas Coulombianas sobre la superficie 
de materia de quark. Tales objetos estan constituidos por la materia en el ultimo estado 
fundamental cuya energi'a es menor que la del micleo de Hierro. Esta materia comprimida 
a una densidad suficientemente alta prodria espontaneamente convertirse en materia de 
quark desconfinada. 

Hasta el presente ningiin experimento ha encontrado quark no confinados, aunque se 
supone que en las Estrellas Extranas los quarks estan libres. A diferencia de una estrella 
normal, la estrella de materia extraha, si tiene este estado auto-ligado de energfa, no 
requeriria de gravedad para unirse. 

Los ANs son los remanentes estelares mas controversiales que habitan el Universe, 
ellos aparecen tras la explosion de una Supernova cuya estrella progenitora es muy masiva 
(M > 10 Mq) o como resultado de la coalescencia de un sistema binario dc dos ENs. El 
colapso gravitacional no puede ser cvitado por la presion del gas dcgcncrado dc particulas 
que constituye la estrella y por ello colapsa apareciendo una singularidad en el espacio- 
tiempo. 

Se presupone que existan en el Universe infinidad de ellos pero observacionalmente no 
han sido detectados hasta hoy. Sin embargo la existencia de ellos ayudaria a comprender 
observaciones astronomicas, emision de rayos X por estrellas binaria y galaxias activas. 
Actualmente, muy pocos astronomos ponen en duda la existencia de los ANs, aunque 
la evidencia de su existencia es indirecta, a partir del comportamiento de otros objetos 
cercanos a los ANs, tales como estrellas brillantes. 
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En 1963, Roy Kerr [17] demostro que en un espacio-ticmpo dc cuatro dimensiones todos 
los ANs, debian tener una geometria cuasi-esfericas determinada por tres parametros: su 
masa M, su carga electrica Q y su momento angular J. 

Con posterioridad, la teon'a de los ANs se nutrio de los trabajos de Hawking, Ellis y 
Penrose de los anos 70 que demostraron varios teoremas importantes sobre la ocurrencia y 
geometn'a de los ANs [18]. 

Como la presencia de campos magneticos intensos cs relevante en los objetos com- 
pactos, su inclusion en las Ecuaciones de Estado (EE) se convierte en un importante tema 
a tener en cuenta para explicar las observaciones, ademas de considerar que tenemos un 
sistema en presencia de gravedad. 

El gas magnetizado de electrones, neutrones y quarks ha sido estudiado previamente 
en las Rcf. [19, 20, 21]. 

Entrc muchos resultados, el mas importante es que se descubre una anisotropi'a entre las 
presiones paralela y perpendicular al campo magnetico, siendo la presion paralela (py > p±) 
mayor. En el caso extreme de campos magneticos muy intensos {B = lO^^T) la presion 
P|| puede hacerse cero y entonces el sistema tiende a colapsar anisotropicamente. Esto 
trac como consecuencia una logica conclusion cualitativa sobre la forma que va a tener 
dicho sistema: el debera tomar la forma de un objcto alargado cn la direccion del campo 
magnetico y a mcdida que cstc campo sea mas intcnso pucs mas alargado se hara cl objcto. 

Es importante seiialar que los estudios citados antcriormcnte se realizaron sin incluir 
el efecto de la gravedad. Ahora, si queremos describir sistemas con alta densidad de 
materia, los cuales puedan producir su propia gravedad (sistemas auto-gravitantes) en una 
proporcion tan elevada, que no podamos ignorar su efecto, entonces podrfamos referirnos 
a la posibilidad de que este sistema pueda describir un objeto estelar. Debido a que los 
sistemas con este tipo de anisotropfa, cn general no son estables en el universo, entonces 
un objeto de este tipo debe colapsar [22] o realizar una transicion a un estado estable [23]. 

Incluir la presencia de gravedad, ahora nos permite obtener informacion sobre el papel 
que ella juega cuando aparecen estos tipos de anisotropias en las presiones, esta es unas 
de las principales motivaciones de este trabajo y nos orienta a responder las siguientes 
preguntas: 

• ^Aparecera o no el colapso cuando la gravedad esta presente en un sistema magneti- 
zado de Fermi? 

• En caso que colapsaran, ^Como sen'a el tipo de singularidad que aparecen'a en estos 
sistemas? 

• ^Como es en general la dinamica de un gas magnetizado de Fermi autogravitante? 

Con el proposito de responder estas interrogantes, el primer objetivo de esta tesis esta 
dirigido a estudiar una fuente de materia densa y magnetizada de fermiones (electrones y 
neutrones), y con el uso de una metrica no estacionaria obtener la evolucion en el tiempo 
del sistema. 

Escogeremos un espacio-tiempo Bianchi I que ha sido utilizado con cxito para la de- 
scripcion de modelos cosmologicos [24]. Esta vez sin embargo estudiaremos un escenario 
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astrofisico, consideraremos un gas degenerado dc fcrmiones, es decir, la temperatura se 
toma igual a cero. La utilidad de estos espacio-tiempos se debe a que nos permiten tomar 
en cuenta las anisotropi'as debidas al campo magnetico. Una extension de este formalismo 
que tome en cuenta la temperatura podria describir una situacion cosmologica [25] . 

Por supuesto la formulacion de un modelo para una estrella EB o una ENs magnetizada 
requiere del use de una metrica mas complicada, que se traduce en un problema numerico 
mas engorroso. Nuestro modelo aunque no muy realista puede dar informacion cualitativa 
muy interesante sobre los procesos dinamicos que ocurren en el nucleo de tales objetos 
astrofi'sicos, en particular para una EB o una ENs. Tal estudio cualitativo podn'a dar luz 
para entender estos sistemas tan complejos. 

Por otro lado a pcsar dc la dcscripcion clasica dc los ANs y sus "pelos" o caracteristicas 
principalcs dcrivadas de la dcscripcion de su M, Q, J (Masa, Carga y Momento Angular 
de Rotacion) [26, 27], descripcioncs de estos entcs a partir dc modcrnas tcorfas de cuerdas, 
le atribuyen situaciones intermedias, donde otras magnitudes pueden ser consideradas, 
como por ejemplo campos escalares, o cargas magneticas. Por otra parte, el estudio de 
la entropia y la perdida de informacion de los ANs es un tema abierto. Ha sido probado 
teoricamente que la entropia de un AN es proporcional al area dc la garganta del AN, 
lo cual ha mantcnido a los Cosmologos y Ffsicos Teoricos a la expectativa y en busca de 
nuevas intcrprctacioncs. 

A pesar de que los objetos astrofisicos deben de ser neutros electricamente existen 
soluciones matematicas que caracterizan a los ANs como objetos cargados electricamente. 
Son soluciones conocidas por el nombre de Reissner-Nordstrom (RN) y tiene el Ifmite 
extremal cuando la masa se iguala a la carga eletrica. La conjetura de Penrose objeta a 
todo AN con carga electrica mayor que la masa, ya que significaria una densidad de energfa 
negativa y por lo que conocemos, la densidad de energfa de todo sistema ffsico conocido 
siempre es positiva. 

La teon'a de cuerdas surgida con la pretension de unificar las interacciones fundamen- 
tales de la naturaleza, ha contribuido a que se desarrollen numerosas tecnicas que con exito 
logran ser aplicadas a otras areas de la Ffsica. Uno de los exitos mayores lo ha logrado con 
la gravitacion. En particular el metodo de Sen o formalismo entropico, fue disciiado origi- 
nalmente como mecanismo entropico en teorfas de Supcrgravedad [28, 29] y extendido con 
exito al estudio de la entropia de ANs con diferentes topologfas. Este formalismo simplifica 
los procedimientos de calculo de entropia conocidos. 

Por tanto el segundo objetivo de esta tesis sera entonces utilizar este mecanismo de 
Sen para estudiar ANs con constante cosmologica y campo electromagnetico en cuatro y 
cinco dimensiones teniendo en cuenta los invariantes de Riemann como conjunto completo 
introducido por Carminati-McLenagha [30, 31]. 

Esta tesis consta de dos partes que se asocian con los dos objetivos que nos hemos 
propuesto con este trabajo. La primera parte aborda el estudio de las fuentes magnetizadas 
autogravitantes dc clcctrones y neutrones. En la segunda parte nos centramos en el calculo 
de la entropia de ANs estudiada a partir del metodo dc Sen. 

Aqui sc estudian ANs en cuatro y cinco dimensiones que ticnen campo electromagnetico 
y constante cosmologica. Este tipo de objeto tiene un interes fundamental en la Fisica de 
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Particulas dondc algunas teon'as de campo predicen la aparicion de soluciones de este tipo 
debido a que la curvatura se modifica por la energi'a y no por la masa del objeto. De igual 
manera conocer estas soluciones, podn'a ser interesante para la Cosmologfa y la Astrofisica. 

La primer a parte de las tesis consta de 4 capitulos. El primero es introductorio y los 
capitulos 2,3,4, contienen la contribucion original del autor a esta tematica. 

En el capftulo introductorio describimos las ecuaciones de la TGR, dos de sus soluciones 
exactas, se explica que son las singularidades y la dinamica de un fluido covariante. 

El segundo capftulo se dedica a estudiar el gas de electrones magnetizado y degenerado 
autogravitante, se escriben las ecuaciones de estado, se estudia la dinamica y se discuten 
los resultados numericos obtenidos. 

El tercer capftulo se dedica a estudiar el gas de neutrones magnetizado y degenerado 
autogravitante, se escriben las ecuaciones de estado, se estudia la dinamica y se discuten 
los resultados numericos obtenidos. 

El cuarto capftulo se dedica a estudiar cualitativamente las condiciones que frenarfan 
el colapso. Esto se hace perturbando el espacio-tiempo Bianchi I y analizando la ecuacion 
de Raychaudhuri [27]. 

La segunda parte de la tesis sc divide en dos capitulos. El primero es introductorio 
y en el se exponen los resultados dc la entropfa de ANs a partir dc la TGR y se explica 
el mecanismo dc Sen. El segundo capftulo conticnc cl aportc original del autor a esta 
tematica. O sea, la aplicacion del mecanismo de Sen al estudio de la entropfa de ANs, 
obteniendo soluciones en 4 y 5 dimensiones teniendo en cuenta el conjunto de invariantes 
de Riemann. Se discuten ademas los resultados obtenidos. 

Finalmente nos dedicamos a enumerar las conclusiones de la tesis y las perspectivas 
del trabajo futuro. 
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Part I 

GASES AUTO-GRAVITANTES 
MAGNETIZADOS 
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1. RELATIVIDAD GENERAL 



1.1 Fuentes de campo en Relatividad Especial 

En esta seccion obtenemos en forma heun'stica las ecuaciones de Einstein, las cuales son 
las ecuaciones de campo de la TGR. Como punto de partida consideramos la Tcoria de 
Relatividad Especial (TRE), la cual sc basa en suponer un espacio-tiempo seudo-cuclideano 
de 4-dimensiones de Minkowski, cuya metrica en coordenadas cartesianas = (x, y, z, t) 
es, 

ds2 ^ r}abdx"-dx^ = -dt^ + dx^ + dy'^ + dz^ , (1.1) 

donde r/ab = diag[— 1, 1, 1, 1]. Las transformaciones de Lorentz, representadas por las ma- 
trices de Lorentz A"j,, describen las transformaciones de coordenadas entre sistemas de 
referenda inerciales y Sfi arbitrarios en el espacio de Minkowski. La forma explicita de 
dichas transformaciones es, 

x"=A"bX^ x"=A%x^ (1.2) 

donde las matrices satisfacen I^\l^^c — '^"c- ^as transformaciones de Lorentz describen 
rotaciones arbitrarias (ver forma explicita en pag 67 de [26]). Toda teoria fisica debe 
ser invariante bajo transformaciones de Lorentz al menos en una region local del espacio- 
tiempo. Esta invariancia asegura el cumplimiento de las leyes de conservacion locales 
(de la energia, momento, etc) que debe cumplir toda teoria fisica en sistemas de referenda 
inerciales arbitrarios. Las transformaciones de Lorentz entre dos sistemas de referenda para 
la 4-velocidad -u", el 4-momento y el tensor de Maxwell F°'^ se expresan en terminos de 
estas matrices de la forma, 

^x" = AV^ p'' = A\p\ F"^ = K^^K^F"^. (1.3) 

El tensor de energia-momento T"'', es un mapeo bilineal simetrico que se transforma 
como un tensor de segundo orden, 

T"^ = M'^K^^f'"^. (1.4) 

Este tensor cuantifica el flujo local de 4-momento p'^ de una fuente dada sobre su- 
perficies = constante, y por lo tanto proporciona toda la informacion dinamica sobre 
las fuentes, las cuales pueden ser particulas masivas, neutrinos, fotones, de ondas electro- 
magneticas, etc. Dada una clase arbitraria de observadores caracterizados por un campo 
vectorial de 4-velocidad u", la forma matematica mas general de T"^ es (ver pag 7, [32]), 

Tab = PUaUb+phab + 2q(^aUb) + T^ab, (1-5) 

donde hab = gab + UaUb es el tensor de proyeccion ortogonal a la hipersuperficie ortogonal 
a u"-. 

Como toda fuente de campo puede ser matematicamente reducida a la forma (1.5), se 
utiliza el concepto formal de "fluido" para describir genericamente a las fuentes. Las compo- 
nentes de T"^ en el marco dc referenda dado por u"' proporcionan las cantidadcs dinamicas 
del fluido (para cualquier fuente) que los observadores asociados a esta 4-velocidad de- 
tectan. 
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• p = TabU^-u^, es la densidad de energia del sistema. 

• p = Tabh"-^ /'^ es, la presion isotropica del sistema. 

• qa = —h^TbcvF es el flujo de energia en la direccion a;". Este tensor es ortogonal a 
la 4-velocidad: QaU^ = 0. 

• T^ab = ^(a^(,)^^c<i ^ describe la presion anisotropica y es tambien ortogonal a la 4- 
velocidad: TTabu"' = 0. 

Un caso caso particular importante es el llamado "fluido perfecto" , el cual se obtiene 
de (1.5) cuando Qa y iTab se anulan. Dada una fuente cualquiera caracterizada por un tensor 
de energia-momento T"^, las ecuaciones de conservacion, 

r"5, = o, (1.6) 

contienen informacion dinamica importante. Por ejemplo, proyectando (1.6) en direccion 
paralela y ortogonal a obtenemos, 

Ecuacion de Euler: UaT"'' = 0, (1.7) 
Ecuaciones de Navier-Stokes: habT^^ ,c = ^- (1-8) 

Para el caso del campo electromagnetico el tensor energia-momento tiene la forma (ver 
pag 89, [26]), 

T"^ = 1-(^f'"'fJ' - jTj'^FikF^''), (1.9) 

y es tambien reducible a la forma (1.5). Las ecuaciones dc Maxwell se obtienen aplicando 
la ley de conservacion (1.6) a (1.9), 

J'' = F''\b, (1-10) 

Fab,c + Fbc,a + Fca,b = 0, (l-H) 

donde el si'mbolo indica derivacion parcial con respecto a .x". La conservacion dc la 
carga viene dada por la proyeccion (1.7) que conduce a (1.9), a la forma dc la ecuacion 
de continuidad J" ^ = 0, donde J" = (/?, J) es el 4-vector de la corriente electromagnetica 
(con p = UaJ"' la densidad de carga y J la corriente de carga). 

1.2 Sistemas de referencia no-inerciales 

Los sistemas dc referencia no-incrcialcs correspondcn a clascs dc obscrvadorcs acelerados 
(por ejemplo, sujctos a rotacion). Las Imeas dc universe dc dichos obscrvadorcs (desde 
un marco incrcial) son curvilmeas, por lo que pucden ser utilizadas para dcfinir un sis- 
tema de coordenadas curvilmeas relacionado con las coordenadas cartesianas mediante la 
transformacion general de coordenadas, 

/=/(x^), (1.12) 

^Los simbolos de corchetes angulares ( ) indican simetrizacion, traza nula y proyeccion ortogonal mediante 
hab para el par de indices correspondiente. Ademas los parentesis en (1.5) indican sirnetn'a. 
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de modo que un campo vectorial se transfer ma como, 

V'' = J''yV^\ (1.13) 

donde, 



J^y^i^^ (1-14) 



es la matriz jacobiana dc la transformacion dc coordenadas. Es importante notar que 
J" h' generaliza las matrices de Lorentz A" ^ a coordenadas generales, de modo que la 
transformacion (1.13) se reduce a las formas (1.3) cuando involucra a dos sistemas de 
referenda inerciales. La metrica del espacio-tiempo de Minkowski en las coordenadas y"' 
se obtiene, a partir de la forma estandar (1.1), mediante la transformacion de un tensor de 
segundo orden, 

ga'iy = Ja'" Jb''^Vcd, (1-15) 

donde Ja' es la matriz inversa de Jc Notese que las componentes de la metrica ga' b' son, 
en general, funciones de las coordenadas curvilmeas y*^ . Es importante que las leyes de la 
fisica pucdan ser expresadas tambien en forma covariante entre observadores no-inerciales, 
de modo que un 4-vector en un sistema de referenda sea tambien 4-vector en cualquier 
otro (y por ende, para tensores arbitrarios y escalares que surgen de sus contracciones) . 

Consideremos una partfcula libre en el espacio de Minkowski con 4-velocidad = 
dx^/dr, donde donde r es el tiempo propio. El 4-momento asociado a dicha partfcula es, 

p'' = mju'' = mJ{l,v^'), (1.16) 



donde /U = 1,2, 3, 7 = 1/-^1 — {v/c)"^, m es la masa en reposo y v"^ = S^i.v'^v'^ con 6^1, = 
diag[l, 1, 1]. Es facil mostrar que la 4-velocidad y el 4-momento de la partfcula libre, desde 
sistemas de referenda no-inerciales, se transferman de acuerdo a la ley (1.13), 

dx°' dy^ y 

u = -5— = J b'-j— = b'U ■ (1-1'j 
dr dr 

Al tratarse de una partfcula libre, su 4-aceleraci6n debe ser nula, lo cual es evidente si 
consideramos el 4-momento (1.16) desde un sistema de referenda inercial en los que es 
un 3-vector constante, 

^j" = ^=m7«" = 0. (1.18) 
dr 

Es evidente que ii"' en (1.18) se transforma como 4-vector bajo transformaciones de 
Lorentz entre sistemas inerciales arbitrarios (al igual que u°' y p"'), de modo que si se 
anula en un sistema de referenda inercial, se debera anular en todos. Sin embargo, la 
4-aceleraci6n no se transforma como un 4-vector entre sistemas de referenda no-inerciales, 
lo cual es evidente al derivar (1.17) con respecto a r y comparar con (1.13), 

du" _ u_ dl/>\hf ' du''' _ „&'„c' , ja -b' n 

u = ——=u hU =J b',c'—^ ; \- J b'—, — — J b',c'U u +J b'U . (,1-iyj 

dr dr d.T dr 

Este hecho es sumamcnte importante, ya que so cspcra que una partfcula libre carezca 
de aceleracion (ya que no esta sujeta a interaccion alguna). Por lo tanto, para poder 



-14- 



generalizar la invariancia de Lorentz en forma adecuada a sistemas de referencia arbitrarios 
(por ejemplo, no-inerciales) , debemos introducir un operador diferencial que generalice a la 
derivada parcial ordinaria, de modo que obtengamos una definicion alterna de aceleracion 
que se transforme como (1.13). Dicho operador es la llamada "derivada covariante", la 
cual denotamos con el si'mbolo ; (punto y coma) . La derivada covariante de un 4- vector F" 
toma la forma, 

V'';b = V\b + T\eV'', (1.20) 
donde F^^ son los sfmbolos de Christoffel, 

r" be = ^g"-'^ {gbd,c + gcd,b - gbc,d) , (i-2i) 

y gab es la metrica (1.15) asociada a las coordenadas curvilmeas y"" . Los sfmbolos de 
Christoffel no se transforman como tensores. Notese que la derivada covariante se reduce 
a la derivada ordinaria cuando consider amos sistemas inerciales, ya que que estos sistemas 
estan asociados a la metrica estandar de Minkowski (1.1) en la que las derivadas de los 
componentes son cero y por tanto 5^ = 0. 

Tomando en cuenta (1.20), es natural generalizar la 4-aceleraci6n dada en (1.19) por, 

u« = u" -bu^ = (u« ,6 + r" bcu'')u\ (1.22) 

la cual podemos demostrar facilmente (aunque es algebraicamente engorroso) que se trans- 
forma como (1.13). 

Por lo tanto, con la mieva definicion (1.22) podemos categorizar a una particula libre 
por {1"" = en sistemas de referencia arbitrarios. Como consecuencia de este hecho, la 
derivada covariante -b se transforma como un tensor de segundo orden, 

V'^.b = J\'J'''bV'',d,, (1.23) 

lo cual es sumamente relevante, ya que las leyes de la fisica se expresan mediante tensores 

y derivadas de estos. 

Esto generaliza la invariancia de Lorentz a sistemas de referencia no-inerciales, sim- 
plemente sustituyendo en las ecuaciones tensoriales la derivada parcial ordinaria por la 
derivada covariante: coma pasa a ser punto y coma. Es importante enfatizar que hemos 

introducido el uso de coordenadas curvilmeas (1.12) para dcscribir observadores acclcrados 
(no-inerciales) en el espacio-tiempo seudo-euclideano de Minkowski. Es necesario recalcar 
que las coordenadas curvilfneas no implican que el espacio-tiempo es curvo, ya que siempre 
es posible encontrar una transformacion de coordenadas que reduzca la metrica (1.15) a la 
forma estandar de Minkowski (1.1). 

Una vez asumida la validez de la invariancia de Lorentz, contenida en la TRE, se 
pensaba hacia principios del siglo XX que el espacio de Minkowski era cl csccnario natural 
para formular las ecuaciones de campo de toda teorfa fi'sica (ya sea en sistemas de referencia 
inerciales o no-incrcialcs) . Como ejemplo importante, sc considcro dcntro dc cstc csqucma 
a la Teorfa Electromagnetica, cuyas ecuaciones de campo (Maxwell) toman la siguiente 
forma en un sistema de referencia arbitrario, 

r = F'^l , Fab;c + Fbc;a + F^a;b = 0. (1.24) 
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El siguiente paso seria elaborar una teoria de gravitacion en el espacio-tiempo de 
Minkowski, que generalizara la teon'a newtoniana y que cumpliera con la invar iancia de 
Lorentz. Sin embargo, los esfuerzos de Einstein por elaborar dicha teon'a fracasaron (ver 
capitulo 7, [26]). 

1.3 Teon'a de la Relatividad General 

Es evidente que un campo gravitacional local puede ser emulado o "cancelado" mediantc 
observadores no-inerciales (como, por ejemplo, un ascensor en caida libre). Sin embargo, 
existen serias dificultades para la descripcion no-local de campos gravitacionales mediante 
los sistemas de coordenadas curvilmeas asociados a observadores no-inerciales. A modo 
de ejemplo, podemos citar el caso de movimientos acelerados en traycctorias hiperbolicas 
(ver pag 166, [26]), dondc se puede dcmostrar que en general cxistc una ruptura de la 
comunicacion entre los observadores fundamentales para distancias mayores que el inverso 
de la aceleracion (ver pag 170, [26]). 

Dadas las inconsistencias en el intento de incorporar la gravedad en el marco teorico 
de la TRE (espacio-tiempo seudo-euclideano de Minkowski), Einstein planteo lo que sera 
considcrado como una de las propuestas teoricas mas innovadoras y rcvolucionarias en la 
historia dc la ciencia: La gravitacion es la manifestacion dindmica de la curvatura del 
espacio-tiempo. Esta sugcrcncia no invalida a la TRE, ya que todo cspacio curvo cs lo- 
calmente piano, por lo que el espacio-tiempo curvo en el cual se manifiesta la gravedad se 
reduce localmente al espacio seudo-euclideano (sin gravedad) de la TRE. Esta propiedad 
es consistente con el poder "anular" localmente a todo campo gravitacional mediante ob- 
servadores no inerciales. De hecho, Einstein llamo a esta propiedad el "Principio de Equiv- 
alencia", segun el cual, la dinamica dc todo sistema, a escala local, puede ser descrita en 
ausencia de gravedad por interacciones que suceden en el espacio-tiempo de Minkowski y 
entre observadores inerciales. 

Bajo este esquema, la gravedad a nivel local se manifiesta (y es equivalente) a acel- 
eraciones detectadas por observadores no-inerciales. Sin embargo, en escalas mas extensas 
es imposible ignorar la gravedad (curvatura del espacio-tiempo), la cual necesariamente 
influye en las ecuaciones de campo de toda teon'a al sustituir derivadas ordinarias por 
covariantes. Este hecho es muy relevante, ya que diversos campos que aparentemente 
no "producen" gravedad (como el campo electromagnetico) son efectivamente fuentes de 
curvatura del espacio-tiempo, y por lo tanto, fuentes de gravitacion. 

La nueva teon'a de la gravitacion basada en asumirla como una manifestacion de la cur- 
vatura del espacio-tiempo es conocida como la "Teon'a General de la Relatividad" (TGR). 
Existen muchas evidencias experiment ales que han validado a la TGR: las traycctorias 
curvas de los rayos de luz, el corrimiento hacia el rojo de un foton bajo la influencia del 
campo gravitatorio, la precesion en la orbita de Mercurio, etc [26] . Todos estos efectos son 
predicciones de la TGR, y son imposibles de predecir o justificar a partir de la TRE. 

En base a los razonamientos anteriores, y tomando en cuenta que todo fiujo de materia- 
energia es una fucntc de curvatura del espacio-tiempo, las ecuaciones de campo dc la TGR 
deben tener como fuente al tensor de energi'a-momento T"^, el cual es un tensor dc scgundo 
orden. Por lo tanto, la curvatura del espacio-tiempo que sera asociada a estos flujos de 
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materia-energia debe ser expresada por un tensor de curvatura que sea del mismo orden. 
Si llamamos G"^ a este tensor, podemos proponer las siguientes ecuaciones de campo, 

Qab ^ i^rj^ah^ ^^ 25) 

donde k es una constante de acoplamiento que sera derivada mas adelante. El tensor G"^ 
debe cumplir con las siguientes propiedades: 

• Debe estar relacionado mcdiante contracciones con el tensor de Riemanian, i,cd, el 
cual es el tensor fundamental de curvatura. 

• Debe ser una combinacion lineal del tensor de Ricci (que resulta de contraer el tensor 
de Riemann) y de la metrica gab- 

• Debe satisfacer la condicion de conservacion G"** -b = 0, ya que esta condicion es 
satisfecha por T"^. 

Es facil demostrar que el tensor Gab que cumple con estas propiedades es llamado 
tensor de Einstein, el cual tiene la forma, 

Gab = Rab - ^gabR + ^Qab, (1-26) 

donde A es la constante cosmologica y la constante de acoplamiento toma la forma k, = 
8irG/c^, la cual sc obticnc rccurricndo al Kmite clasico newtoniano al asumir que las ve- 
locidadcs asociadas a los flujos dc materia-energia son mucho menores que la de la luz. Las 
ecuaciones dc campo (1.25) con G"'^ dado por cl tensor de Einstein (1.26) son conocidas 
como "Ecuaciones de Einstein" , y son las ecuaciones de campo de la TGR. 

Los tensores que aparecen en las ecuaciones de Einstein son simetricos, de modo que 
en 4 dimensiones tienen 10 componentes independientes. Dada la libertad de eleccion de 
las cuatro coordenadas del espacio-tiempo, las ecuaciones independientes se reducen a 6. 
Las ecuaciones de Einstein son un sistema de ecuaciones difercnciales parciales no-lineales 
de alta complejidad, por lo que es dificil encontrar soluciones exactas de ellas. 

SOLUCIONES EXACTAS 

En la siguientes secciones presentamos dos clases de soluciones exactas de las ecuaciones de 
Einstein. Asimismo, se discute el concepto de singularidad de curvatura, y se derivan los 
parametros cinematicos que describen localmente la evolucion de elementos de un fluido 
arbitrario. 

1.4 Solucion de Schwarzschild 

La sohicion de Schwarzschild fue la primcra solucion exacta de las ecuaciones de Einstein 
(1.25). Esta es una solucion de vacio (T"^ = 0) que asumc simetria esferica, por lo que es 
adecuada para describir el exterior dc un objeto csfcrico. Es posiblc demostrar que dicha 
solucion es ademas estatica (existen coordenadas en que la metrica no depende del tiempo) 
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y linica (Teorema dc Birkhoff). La solucion de Schwarzschild en su maxima extension 
analitica permite describir un AN esferico (por lo tanto, sin momento angular) y sin carga 
electrica. 

Suponiendo coordenadas esfericas = {t, r,0,4>), la metrica de Schwarzschild tiene la 
forma, 



= ^ {2H- S3){Tm - 2Tu,p)P - m{Tp + Tu) , (1.27) 

^3,^ = 2PTm - 353^, (1.28) 

n,T = l5TM + \{Tu-Vp)-m^, (1.29) 

= 2^(53 - 2H), (1.30) 

ds^ = -(l- de + , t"L ^ + {d^^ + sill' ^ (1-31) 



r J (l-2M/r) 

donde M es la masa del objeto esferico en unidades de longitud (ya que suponemos Gjc? = 
1). Notese que a distancias lejanas de la fuente (r S> M) la metrica (1.31) tiende al espacio 
piano de Minkowski (metrica (1.1)). 

Es cvidcnte que la metrica de Schwarzschild (1.31) se torna problematica cuando r — >■ 
2M, ya que y ^^r ^ oo. Como M es una longitud, se le llama el "radio de 

Schwarszschild" . ^ Sin embargo, es facil mostrar que el comportamiento problematico 
de (1.31) no corresponde a condiciones fi'sicas singulares, ya que todos los escalares de 
curvatura son finitos en r = 2M, como por ejemplo. 



,ahcd 48 M2 



RabcdR'''''' = —^. (1.32) 



^6 

Ademas, toda clase de observadores en caida libre alcanzan el valor r = 2M, y Uegan a 
r = en un tiempo propio finito. El problema con r = 2M se debe a que las coordenadas 
(t, r) no son adecuadas para la dcscripcion dc rcgioncs del cspacio-tiempo cercanas a este 
radio. Debido a que la metrica (1.31) no es valida para r < 2M, esta no da una descripcion 
completa del espacio-tiempo. Podriamos omitir por el momento el valor problematico 
r = 2M y extender estas coordenadas al rango < r < 2M, de modo que la metrica 
completa es (1.31) para r > 2M y para < r < 2M tomar la forma, 

(]r^ ( 2M \ 

ds^ = -TTTTTT -T + 1 di^ + {dO"^ + sin^ 9 d<t)'^), < r < 2M, (1.33) 



(2M/r - 1) V ^ 

con lo que el caracter causal de (t,r) se invierte: para r > 2M, la coordenada t es de 
tipo tiempo y r de tipo espacio, con lo cual todo observador puede seguir una Imea de 
universo dada por r = constante, t = arbitrario (observador estatico). Sin embargo, para 
< r < 2M, r pasa a ser una coordenada tipo tiempo, por lo que ningiin observador pucdc 
seguir la h'nea de universo r = constante (la cual queda fuera del cono de luz). En otras 



*EI radio de Schwarzschild corresponde al radio que deberia tener un objeto esferico para que la velocidad 
de escape (newtoniana) fuera igual a la de la luz. 
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palabras, para un observador en esta region r dcbc disminuir hasta r = 0. Evidentemente, 
el radio de Schwarzschild marca este cambio en la causalidad de r y i. 

Debido a que la extension de las coordenadas discutida en el parrafo anterior es singu- 
lar (matematicamente hablando) en r = 2M, se utilizan otras coordenadas para entender el 
comportamiento del espacio-tiempo alrededor de este valor. Por ejemplo, Kruskal y Szck- 
eres propusieron el sistema de coordenadas dado por las relaciones {t, r, 6, (f>) (T, R, 6, (f>) 
dadas por, 

las cuales transforman la metrica (1.31) a la conocida forma de Kruskal-Szekeres, 

ds" = exp (^) {-dT^ + dR^) + {d9^ + sin^ 9 dcj>^), (1.36) 

donde ahora r = r{T, R) se debe obtener de la solucion del sistema de ecuaciones (1.34) y 
(1.35). 

En csta forma de la metrica, las coordenadas (T, R) cstan bicn definidas para todo 
valor real. Las gcodcsicas milas son lincas rcctas a ib 45 grades dc inclinacion, por lo que 
podemos examinar la estructura causal de este espacio-ticmpo. Inmediatamente podemos 
identificar cuatro regiones de interes, 

• Region I o region exterior de la solucion de Schwarzschild, caractcrizada por R > \T\; 
que corresponde a la metrica original (1.31), y que describe al campo gravitatorio 
externo (vacio) de un objeto esferico. Vease en el diagrama del Apendice B, Figura 
6. 

• Region II o region de AN de la solucion de Schwarzschild, caracterizada por Vl + -R^ > 
T > \R\. La direccion de las curvas r = constante son ahora tipo tiempo, por lo que 
todo observador dentro del cono de luz debe evolucionar hacia r = 0. A esta region se 
le denomina "Agujero Negro" (AN), ya que los observadores dentro de ella no pueden 
transmitir causalmente informacion alguna a observadores en la region I (sus senales 
luminosas no Uegan a r > 2M). 

• Region III o region de Agujero Blanco de la solucion de Schwarzschild, caracterizada 
por -VTTW < T < -\R\. Es analo ga a la region II, excepto que la direccion 
causal es hacia r creciente (desde r = hasta r = 2M). Observadores en esta region 
no pueden ser infiuenciados por los de las demas regiones, pero ellos pueden afectar 
causalmente a estas. 

• Region IV o region exterior paralela, caracterizada por R < —\T\. Esta region es 
identica a la region I, pero esta desconectada causalmente. 

El valor r = 2AI corresponde a una 2-esfera gencrada por segmentos de gcodcsicas 
nulas (vease diagramas del Apendice B), y se le denomina "horizonte de Schwarzschild", 
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ya que encierra a la region I (el AN). Como ya se comento, el espacio-tiempo es regular 
en r = 2M, sin embargo, los escalares de curvatura, como por ejemplo (1.32), divergen 
en el Iimite r — > 0, lo cual se puede asociar a fuerzas de marea colosales que destrozarian 
a todo objeto fisico. ^ Como los escalares de curvatura divergen, no es posible extender 
el espacio-tiempo a r = 0, por lo que toda curva causal (geodesica temporal o nula) debe 
terminar en r = 0. 

El valor r = corresponde entonces a una singularidad de curvatura, la cual puede 
ser definida como un conjunto de valores de las coordenadas = C^, tales que todos los 
escalares del tipo (1.32) divergen conforme — >■ a lo largo de las geodesicas causales 
en valores finitos del parametro afm de dichas curvas. Bajo este criterio, la solucion de 
Schwarzschild en su maxima extension analitica tiene una singularidad de curvatura en 
r = 0. 

1.5 Solucion de FLRW 

Al asumir la validez del Principio Cosmoldgico ^, la geometi'a del espacio-tiempo viene a 
ser descrita por la metrica de Robertson- Walker (RW), 

dondc k = 0, ±1 rcprcscnta la 3-curvatura (cero, positiva o negativa) de las secciones 
espaciales t = constante. 

Los espacio-tiempos asociados a esta metrica, para diferentes fuentes del campo, son 
los llamados modelos Priedman-Lemaitre-Robertson- Walker (FLRW), los cuales son am- 
pliamente utilizados en la Cosmologia contemporanea ^. 

El tensor de energia-momento compatibc con la metrica (1-37) en una representacion 
comovil ® es necesariamente de la forma del fluido perfecto, 

T" I, = p u" vi,+ph\ = diag[-p(i), Pit), p{t), p{t)], (1.38) 



^Las fuerzas dc marca se denominan en ingles "tidal forces". 

®EI Principio Cosmologico es la hipotesis principal de la Cosmologia moderna, la cual Ueva a una inter- 
pretacion teoricamente consistente a un niimero creciente de evidencias observacionales. El principio implica 
que el universe a escalas lo suficiciitcrncnte grandes (del orden de cientos de megaparsecs) es isotropico y 
homogeneo. La isotropi'a exige que las cantidades observables sean independientes de la direccion de ob- 
servacion. El Principio Cosmologico sugiere que esta propiedad sea valida (a grandes escalas) para todo 
observador cosmico, lo cual implica a su vez la "homogeneidad" , o sea, que las variables fi'sicas y geometricas 
sean independientes de la posicion de los observadores, dependiendo solo del tiempo cosmico (el cual es el 
tiempo propio de estos). Estas propiedades solo son satisfechas por la metrica Robertson- Walker. 

^Las iniciales FLRW corresponden a Friedman, Lemaitre, Robertson y Walker, que son los flsicos y 
matematicos involucrados en el estudio de cstos modelos. El ffsico ruso Alexander Friedman obtuvo entre 
1922-1924 la primera solucion de las ecuaciones de campo (para polvo) con la metrica RW. Esta solu- 
cion fue generalizada extensamente en 1927, para diversas fuentes de campo y valores de la constante 
cosmologica, por el jesuita belga Georges Lemaitre. Entre 1935 y 1936, los matematicos ingles y norteamer- 
icano, Robertson y Walker, estudiaron las propiedades formales de la metrica RW, probando rigurosamente 
que esta representa a la linica geometn'a del espacio-tiempo seudo-riemanniano que satisface al Principio 
Cosmologico. 

®En una representacion comovil, la 4-velocidad asociada a la metrica RW es = (5" o- 
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donde p = UaUhT"'^ Y P = (l/3)/ia6^"'' pueden ser identificados con la dcnsidad de masa- 
energia total y la presion medidas por un observador comovil (la metrica RW en esta 
representacion no admite presiones anisotropicas ni flujos de energi'a). Las ecuaciones de 
Einstein son, ^ 

87rp = ^^^^ = -GS, (1.39) 



i e + k _ 2R 



^^P = ^ - = G\ = G% = G'^^. (1.40) 



Las ecuaciones (1.39)-(1.40) pueden ser complementadas con la ley de conservacion UaVbT"'^ 
0, 

p + 3(p + p)| = 0, (1.41) 

la cual provee una condicion de integrabilidad. Por tanto, tenemos dos ecuaciones efectivas 
(cualesquiera dos de las tres ecuaciones (1.39)-(1.41)) para tres incognitas {p,p,R), por 
lo que hace falta especificar un vinculo extra, o "ecuacion de estado", para obtener un 
sistema de ecuaciones determinado. La ecuacion de estado no puede ser proporcionada por 
la TGR, sino que debe ser sugerida por consideraciones fisicas concernientes a la fuente de 
materia. Es litil reescribir las ecuaciones (1.41) y (1.39) como, 

^^{pR')+p^^{R')=0, (1.42) 
r2 = ^p{R) i?2 - k. (1.43) 

Dada una relacion entre py p, podemos resolver (1.42) y obtener una relacion funcional 
del tipo p = p{R), con la cual la ecuacion (1.43) se convierte en una ecuacion de evolucion 
efectiva, conocida como "ecuacion de Friedmann" cuya integracion sera de la forma R = 
R{t). 



1.5.1 Fuentes de materia— energia 

Una ecuacion de estado de uso sumamente comiin es la ecuacion barotropica p = p{p), y 
en particular la conocida "ley gamma" dada por p = (7 — l)p, donde 7 es una constante. 
Ciertos valores especificos de 7 se pueden asociar a situaciones fi'sicas, por ejemplo (7 = 4/3) 
modela radiacion incoherente, (7 = 1) modela polvo (gas de presion nula). El caso 7 = 
modela la energi'a del vacio (constante cosmologica) y lleva a los modelos de de Sitter. Otros 
tipos de fuentes de materia, tales como campos escalares o cuerdas cosmicas, pueden ser 
acomodados en el tensor de energia-momento (1.38), aunque en estos casos la 4-vclocidad 
solo tiene sentido formal. Diversas fuentes de materia son tratadas en el Apendice C. 

^Consideramos A, la constante cosmologica, como una forma adicional en la densidad de energi'a 

cumpliendo con el vinculo pa + pA = 0, lo que se conoce como "energi'a del vacfo" o "vacfo cuantico" 
y puede ser incorporado en las cantidades totales p y p reemplazando p + A y p — A. 
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1.5.2 Singular idades de curvatura 

De todas las rclaciones tipo p = p{R) obtenidas cn los ejemplos del Apendice C, es evidente 
que p diverge para cuando i? — t- . Esto indica que los modelos FLRW, con las ecuaciones 
de estado consideradas, inician su evolucion en una singularidad de curvatura escalar, 
marcada per i? = 0, y caracterizada per la divergencia de todos los escalares construidos a 
partir del tensor de Riemann y de Ricci. Los modelos FLRW comienzan por este estado de 
singularidad y se caracterizan por estados de evolucion primordialcs a muy altas dcnsidades 
y tcmpcraturas (en los casos en que la temperatura pueda ser definida). Esta situacion 
es conocida con el termino ingles dc "Hot Big Bang scenario" y la singularidad inicial cs 
llamada singularidad de "Big Bang". El factor de escala es una medida de las distancias 
locales relativas entre observadores fundamentales, de modo que si i? ^ estas distancias 
empiezan a disminuir y terminan por anularse en el limite. Las Imeas de universo de 
observadores fundamentales convergen y asi se caracteriza el Big Bang como un punto 
caustico. 

1.6 Forma covariante de la dinamica de fluidos y colapso gravitacional 

Podemos definir, para toda congruencia de observadores cosmicos asociados a un campo de 
4-velocidades it"(x^), parametros covariantes que determinan el comportamiento cinematico 
de elementos de volumen de un fluido dados por las mismas Imeas de universo de los 
observadores en las hipersuperficies 3-dimensionales espaciales y ortogonales a ■u". Dichos 
parametros surgen de la descomposicion irreducible y covariante de la derivada covariante, 

Q 

Ua;b = -^hab + Cfab + ^ab " UaUi, , (1.44) 

donde, 

• ^ = ;a es el escalar de expansion, el cual describe el cambio isotropico del volumen 
de los elementos de fluido. 

• Ua = Ua-cU"^ Gs la 4-aceleraci6n, la cual determina que tan difercntc cs la evolucion de 
las Imeas de universo de los observadores en caida libre, ya que en esta ultima dichas 
Imeas de universo serian geodesicas {iia = 0). 

• Cafe = h'^ ah'^ bU(^c;d) ~ (^/3)/ia6 + UaUb CS el tensor dc dcformaciones o cizalladura, el 
cual determina la dcformacion en las direccioncs dadas por sus autovalores (sin cambio 
de volumen) de los elementos del fluido. La notacion indica simetrizacion sobre 
los indices a, b. 

• ^ab = ah'^ b U[c;d] es el tensor de vorticidad, cl cual determina la rotacion dc ele- 
mentos del fluido (a volumen constante) en un eje definido por sus autovalores. La 
notacion [^j.^j indica anti-simetrizacion sobre los indices a, b. 

^"E1 unico caso de fuente que logra escapar a este analisis, es el caso de energia del vacio, en donde la 
densidad de energi'a es igual a la constante pt'''"^^") = A/Stt. 
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Supondremos que la dinamica dc los obscrvadores cosmicos viene dada por las ecua- 
ciones de Einstein G"'' + Kg""^ = SirT''^^, donde A es la constantc cosmologica y el tensor 
de energia-momento toma su forma mas general posible dada por (1.5). Es posible de- 
mostrar que el tensor de Ricci Rat asociado al espacio-tiempo de los observadores satisface 
la siguiente ecuacion, 

02 

e + — = -Rab + au" ;6 + 2(a;2 - a^) + u\ (1.45) 

donde ijp- = (l/2)a;a5 a;"^ y cr^ = (l/2)(7ttb cr"^. Si consideramos el tensor de energia- 
momento (1.5), y las ecuaciones de Einstein que implican la relacion, 

Rab ll" + ^ i?" a = StT Tab = ^TT p - A. (1.46) 

Dado que —R"' a = —T"' a = P + ^P, podemos transformar (1.45) en su forma conocida 
como "ecuacion de Raychaudhuri" , 

y = -47r(p + 3p) + A + 2(a;2 _ a^) + /i& ^u'^ + u"", (1.47) 
donde hemos definido al factor de escala L como el escalar que satisface, 



L 
L 



(1.48) 



La ecuacion de Raychaudhuri (1.47) ofrece una simple validacion general de nuestras 
expectativas intuitivas de como el efecto dinamico de la gravitacion como curvatura del 
espacio-tiempo (dado por Rabu"'u^) se traduce tambien en el efecto de una "fuerza atractiva" 
que se manifiesta por la aceleracion L/L de las distancias caracteristicas (L) entre dos 
elementos de fluido definidos por u". 

• Los terminos en (1.47) que favorecen el colapso para r > tq son p + 3p > y la 

deformacion o"^ . Vale la pena remarcar que a la condicion p + 3p > sole dcnomina 
"condicion de energi'a fuerte" , y es satisfecha por la mayon'a de las fuentes de masa- 
energia conocidas. 

• Los terminos que en (1.47) que no favorecen (se oponcn) al colapso para r > tq son 
p + 3p < (violacion de la condicion de energi'a fuerte), vorticidad no nula w^, 4- 
aceleracion (Imeas de universo que no son geodesicas) y una constante cosmologica 
A no nula y positiva. 

El simple analisis cualitativo de la ecuacion de Raychaudhuri que hemos hecho es la 
base de los llamados "teoremas de singularidad" (ver capitulo 9, [27]). En la practica, 
los elementos de fluido, cuyos terminos en (1.47) favorecen el colapso, necesariamente 
evolucionan en un tiempo propio, r, finito hacia un estado singular caracterizado por 



^^Aquf hemos tornado k = SirG/c*, con G/c^ = 1. 

^^Dicha ecuacion fue descubierta independientemente por Amal Kumar Raychaudhuri y por Lev D Lan- 
dau. Para ver demostracion a partir de las geodesicas cercana revisar Apendice D. 
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9 —7- — oo y volumen espacial nulo (colapso gravitacional) . Debido a (1.45), los escalares 
de curvatura (asi como la densidad y presion) tenderan a infinite para ese tiempo propio 
finite, lo cual define a una singularidad de curvatura. 

En la TGR, existen muchas versiones de los teoremas de singularidad de Penrose- 
Hawking. Muchas versiones establecen que si existen superficies nulas atrapadas y la 
densidad de energia es no negativa, cntonces alli existen geodesicas de longitud finita las 
cuales no pueden ser extendidas. Estos teoremas, estrictamente hablando, prueban que al 
menos una geodesica temporaloidea o nula puede ser finitamente extensible, solamente 
hacia el pasado, pero existen casos en los cuales las condiciones de estos teoremas se 
obtienen de tal forma que todos los caminos en el espacio-tiempos hacia el pasado terminan 
en una singularidad. Por otro lado la conjctTira dc Penrose asegura como linica singularidad 
desnuda, la singularidad inicial dc Big Bang al principio del origcn de nuestro Univcrso, 
en otro caso se violarian las condiciones de encrgi'a para cicrtos sistemas fi'sicos. Si nos 
guiamos por esta conjetura, que no ha sido demostrada, entonces solo nos queda pensar 
en singularidades de colapso donde la singularidad esta escondida dentro del horizonte de 
eventos del sistema. 
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2. ESTUDIO DE UN GAS MAGNETIZADO Y AUTO-GRAVITANTE 
DE ELECTRONES 



2.1 Introduccion 

Indiscutiblemente, los objetos compactos como lo son las EBs pueden ser interpretados 
como laboratories astroffsicos donde la existencia de altas densidades y campos magneticos 
fuertes plantean un reto para su estudio. En principio suponer que una seccion de volumen 
de un objeto compacto se encuentre inmerso en un campo magnetico fuerte no implica 
estar muy lejos de la realidad. Ademas, cuando decimos campos magneticos fuertes nos 
referimos a campos muy superiores a los obtenidos en laboratories terrestres. Es muy 
interesante estudiar dicho sistema, las Imeas del campo magnetico y su relacion con la 
presencia de gravedad. El posible vinculo con las condiciones de colapso del sistema en 
general aumentan'a la importancia de dicho estudio. 

Por otra parte, gran cantidad de modelos ya sean teoricos, fenomenologicos y niimericos 
han sido propuestos para la descripcion de estrellas EBs magnetizadas. Entre los modelos 
mas simples (no magnetizados), esta el modelo de Landau (ver [4, 33, 34]) que fija la 
masa en un Hmite finito para cstos objetos. Gcneralmente un modelo consistente para 
aplicarlo sobrc una cstrclla EE magnctizada debe requerir considerar codigos numcricos 
hidrodinamicos modclando las ecuaciones de Einstein-Maxwell con, al menos, condiciones 
de simetricas y de frontera. 

Es conocido que la teoria de Einstein no produce ecuaciones de estado, por tanto dichas 
ecuaciones se deben importar desde otra teoria, por ejemplo, la estadistica, la teoria de 
campo, o tambien pueden tener origen fenomenologico. En el caso especial de las EBs el 
caso mas simple de ecuacion de estado scn'a suponer un gas ideal de fermiones magnetizado 
y a temperatura cero, pues estas estrellas en su etapa final se enfn'an. Por tanto, considerar 
un gas degenerado no resulta ser una mala aproximacion. 

Sin embargo, se puede obtener informacion importante sobre la dinamica local de un 
gas de Fermi magnetizado y auto-gravitante a partir de la evolucion de esta fuente en un 
contexto mucho mas simple de la geometrfa espacio-tiempo. Lidiar con espacios-tiempo mas 
realistas nos harian dcscmbocar en sistemas de ecuaciones de compleja solucion numerica. 
Por tanto, en este capi'tulo nos proponemos estudiar la dinamica de un gas magnetizado y 
auto-gravitante de electrones a temperatura cero. El presente capi'tulo generaliza el trabajo 
previo hecho en [22]. 

La ecuacion de estado que consideramos para describir un volumen local de una estrella 
EB magnetizada es la obtenida en [35, 36, 37] y posteriormente es desarrollada para fuentes 
mas diversas en [19, 20, 21, 38, 39]. En estos trabajos se discuten ecuaciones de estado para 
gases fuertemente magnetizados que contienen presiones anisotropas. En estos casos las 
presiones, p\\ paralela y p± perpendicular al campo son diferentes. La presencia del campo 
magnetico provoca estas anisotropfas y la posibilidad de "colapso" se sugiere en la direccion 
perpendicular al campo. Estos resultados nos Uevan a pensar en extender tales estudios 
al marco de la TGR para configuraciones magnetizadas, prcguntandonos si se mantendra 
dicho "colapso", y si este podn'a dar paso a una singularidad tipo anisotropa a lo largo 
de la direccion del campo magnetico. A este tipo de singularidades se les define como 
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singularidades tipo "Imea" o "cigarro", en contraste con las singularidades tipo "punto" 
de sistemas isotropos y altamente simetricos^^. 

En aras de investigar la dinamica local y el colapso de un gas magnetizado en el 
contexto de la TGR, nosotros consideramos una de las configuraciones compatibles con 
presiones anisotropas. O sea el espacio tiempo Bianchi-I el cual es descrito en terminos de 
una metrica de Kasner. 

El argumento mas claro para usar la geometria Bianchi es que esta puede proveer una 
descripcion simplificada dc la dinamica local de un clemento de fluido dc un gas magneti- 
zado en una configuracion realista, asf este volumen desde el punto de vista termodinamico, 
intercambia parti'culas y energia con el resto del sistema. Entonces, podemos considerar a 
este volumen como parte de una distribucion gran canonica asociada a todo el gas. Ademas, 
supondremos que el elemento se encuentra lejos de la frontera del sistema donde los efectos 
de frontera no influyan en la dinamica de dicho volumen. O sea, en el caso que sea un 
objeto compacto, dicho volumen sc cncontraria en el centro del sistema, o muy cerca del 
centro donde el campo magnetico sc puede suponer que se encuentra apuntando en una 
direccion preferencial. Evidentemente, un tratamiento como este es aproximado y no se 
puede asociar a la descripcion del objeto compacto real donde el campo magnetico varia 
con el radio y la aproximada geometria esferica del objeto tambien serfa importante. No 
obstante cstc cstudio ticnc cl valor de acercarnos al estudio dinamico de un volumen central 
elemental de dichos sistemas. 

Una vez que tomamos el espacio-tiempo Bianchi-I con metrica Kasner como la metrica 
del campo asociado a un gas magnetizado, el sistema de ecuaciones de Einstein-Maxwell se 
reduce a un sistema autonomo de cuatro ecuaciones diferenciales. Este sistema representa 
un espacio de fase de cuatro dimensiones que puede ser estudiado tanto numericamente 
como cualitativamente, usando las tecnicas ya estandarizadas para sistemas dinamicos. 
Veremos que las magnitudes ffsicas que caracterizan el espacio de fase seran conveniente- 
mente adimensionalizadas. Ellas seran la componente independiente del tensor de defor- 
maciones S^n, el escalar de expansion %, el campo magnetico /3 y el potencial qui'mico 
ji normalizados respectivamente. Todas las otras cantidades pueden ser expresadas en 
terminos de estas cuatro cantidades basicas. El examen numerico del sistema aclara el tipo 
de singularidad de colapso que puede tener el volumen en dependencia de las condiciones 
iniciales. 

Las ecuaciones de estado que usaremos aquf son estrictamente validas para densidades 
del orden de (al menos) 10^° kg y 10^^ kg esperadas para objetos compactos como 
las EBs y ENs. Un gas fuertemente magnetizado y altamente degenerado en una EN esta 
en un estado muy cercano a la superfluidez con conductividad infinita (ver pag 291, [40]). 
En estas condiciones cl papcl dc la viscosidad es menor, aunque uno podrfa considerar la 
posibilidad de disipacion o fenomenos de transporte, tales como disipacion de la energfa 
rotacional en ondas electromagneticas o gravitacionales (ver, [40, 41, 42]). Sin embargo, 
incluso si la viscosidad no es significativa (al menos para ENs), la razon mas importante 
del por que la despreciaremos (asf como otros efectos disipativos) en nuestros estudios, es 

'^^Ver definiciones de tipos de singularidad en [24]. 
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por mantener la simplicidad del modclo. O sea, queremos mantener el problema tratable 
matematicamente. Nosotros pensamos que el caso de equilibrio termico es suficiente como 
primera aproximacion para este modelo y dejamos el caso de fenomenos de transporte 
disipativo para trabajos futures. 

2.2 Ecuaciones de Estado de un gas magnetizado de electrones 

Como ya hcmos visto en la introduccion, una de las mas importantes propiedades de un 
sistcma formado por un gas magnetizado de fcrmiones, es que cxhibc una anisotropi'a en 
las presiones. Esta anisotropi'a puede ser directamente asociada con la cuantizacion de los 
niveles de energia por la presencia del campo. Clasicamente, el electron orbita en el piano 
perpendicular al campo magnetico formando ci'rculos, con velocidad angular constante. 
Su movimiento pueden ser descompuesto en el movimiento de dos osciladores armonicos 
simples. Cuando este oscilador armonico simple es cuantizado, sus niveles de energia vienen 
dados por la expresion, 

£;(p^,n,r) = ±mc2[l + (^)2 + 2/3(n + r-l)]V2_ (2.I) 

Este espectro de encrgi'a para los electrones se obtiene resolviendo la ecuacion de Dirac 
en presencia de campo magnetico constante en la direccion z. En (2.1) p3, es el momen- 
tum lineal en la direccion z, P = B/Bc el cociente entre el campo magnetico y el campo 
magnetico critico Be = w?(? [eh = 4.414 x lO^T. Los signos + y - se refieren a electrones 
y a positrones respectivamente, r = l,2yn = 0,l,2..., son los mimeros cuanticos de spin 
y mimero cuantico principal de la orbita del electron. 

En el limite de n muy grande, el ultimo sumando en la expresion de la energia corre- 
sponde con la suma de los cuadrados de los momentum lineales en las direcciones x y y, es 
decir, 

2|-(n + r-l)^(^)2 + (^)^ (2.2) 
Be mc mc 

el gas se transforma en un gas de Fermi ordinario. O sea, se puede definir el parametro 
adimensional, 

Asi las propiedades magneticas del gas son clasicas si (5 >> 1 o cuanticas si 5 << 1. 
Los campos magneticos de interes para el gas de electrones describiendo una EB, estan 
cn cl rango de 10^ T y pueden estar presentes durante procesos de colapso gravitatorio de 

objetos estclares. 

Los gases magnetizados de Fermi no relativistas han sido estudiados extensivamente 
en el marco de la fi'sica estadistica y sus propiedades mas importantes estan asociadas con 
el diamagnetismo de Landau y el paramagnetismo de Pauli. 

En adicion a esto veremos la forma que toman las ecuaciones de estado para el caso 
relativista. El Gran Potencial Termodinamico que se halla con la ecuacion, 

n = kT In Z, (2.4) 
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donde k es la constante dc Boltzman, T la temperatura, Z = Tr{p) cs la funcion Gran 
Potencial Termodinamico, p = Q-i^-t^^)/^"^ ^ }j es el Hamiltoniano, p es el potencial 
qmmico y iV es el operador del mimero de parti'culas. 

El tensor de energia-momentum total, donde en la diagonal principal estan las presiones 
y la densidad de energi'a, lo obtendremos haciendo la promediacion macroscopica, 7^,/ =< 
fi^icro) ^ ^^^^g f (micro) ^ (^^/^^a J + {dC/O'^l,) , - 6^,1, entonces para la 

funcion de particion Q = —{kT)~^ln< ^fi <i^i I d^^^i^i^^) >. 14 g^^g electron-positron es 
relevante en el contexto astrofisico, y el Gran Potencial Termodinamico para el sector 
electron-positron tiene dos terminos de la forma, 

^ = nse + ^Ve, (2.5) 

donde el primer tcrmino cn cl lado derecho, es la contribucion estadi'stica y cl scgundo la 
contribucion del vacio [20]. Explfcitamente, tenemos que en la aproximacion de un lazo, 

= -r^E"- / dp3ln[/+(Me,0r(/^e,0], (2.6) 

donde ^ = l/ksT, an = 2—Son y f^il^e, = (l+e"^^^^'^^)^) representan respectivamente, 
las contribuciones para particulas y antiparticulas. El termino de vacio es dado por la 
expresion. 



= 51 / P±dpxdp3£^e, (2.7) 

el cual es divergente, pero puede ser renormalizado, y para campos de intensidad B < 10^ T, 
su contribucion es irrclevante [20] para campos de interes en los objetos astrofisicos actualcs, 
por tanto nosotros no tomaremos en cuenta este termino. Tampoco incluiremos en Q-se^ la 
contribucion clasica (o de arbol) asociada al campo magnetico, pues seguimos la corriente 
de publicaciones en donde este termino se omite (ver [21]) , aunque posteriormente en 
el capitulo 4, si lo tcndrcmos en cuenta. 

La ecuacion (2.6) pucde ser integrada facilmente para el caso degenerado (T = 0). En 
ese caso la funcion de distribucion de los electrones se convierte en la funcion paso unitario 
y la de los positrones es cero. El tensor energia-momentum asociado a un gas de Fermi en 
presencia de un campo magnetico externo y constante en el tiempo, toma la forma [21], 

-^"^ = + E ^"1^)'^" 4 + ^^"'^0,^2 (2.8) 

asi que en el limite de campo magnetico nulo obtenemos el tensor del fluido perfecto, 
T" 6 = p(5" & — (p + U)S°- 4(5^ b. Las componentes del tensor (2.8) son, 

T% = p\\ = -n = p, (2.9) 

T\ = T\=p± = -n-BM=p-BM, (2.10) 
T\ = -U = -TS-iieN-n, (2.11) 

Aquf, los corchetes angulados < > significan la integral continual de Feynman en la teoria de Matsubara. 
^^Otros autores si incluyen este termino clasico, que flnalmente introduce cambios en las direcciones de 
colapso a nivel newtoniano (ver, [43]). 
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donde = p, p± son las presiones, paralela y perpendicular al campo magnetico, U 
la densidad de energia, S la densidad dc cntropi'a, N la densidad de particulas y la 
magnetizacion. Para un objeto compacto, en especial para una estrella EB la temperatura 
T << Tp Es por ello que una descripcion de estas estrellas pueden hacerse suponiendo 
un gas degenerado de particulas. 

Todas las variables termodinamicas del sistcma pueden ser obtenidas a partir del Gran 
Potencial Termodinamico $7, y las ecuacioncs dc cstado para el limite de campo no nulo y 
degenerado (T = 0), tendran las forma expKcita, 



p±=p-BM, p = X/3Tp{fie,P), 

u = xpru{fie,/3), MB = xprM{fie,P), 

A = mc^/(47r^Ac) y las funciones T^, seran, (ver [19]), 

Tp = Y ^/Ug ~ + ^arcsinh j + ^ X] "'^(^ri - K - Cn), 



= 

Tu = 

donde, 
an = 
a„, = 



Xa„(an -bn- 2Cn), 
n=0 

«o (^l4 ~ ^) ~ ^arcsinh + X] '^ri{an + bn + Cn), 



2 -Son, n = 0, 1,... 
IJ-eV f4 - 1 - 2nP, bn = In 



(^e + an/lJ-e 



, Cn = 2nPbn, 



(2.12) 
(2.13) 



(2.14) 
(2.15) 
(2.16) 

(2.17) 
(2.18) 

(2.19) 



2/3 

y I[X] denota la parte entera de este argumento X. 
2.3 Ecuaciones de Einstein— Maxwell 

La metrica de Kasner es una de las metricas mas simples que permiten anisotropias 
en las presiones independientemente del origen de las mismas. En el caso de presiones 
anisotropicas asociadas a una fuente magnetica, ella puede ser utilizada. Esta metrica 
viene dada por, 



dr + A\t) dx^ + B\t) dy^ + C\t) dz^ 



(2.20) 



Ella esta asociada con un espacio-tiempo Bianchi-I "no-inclinado" [24]. Para una 
4-velocidad comovil n" = (5", donde u"'Ua = —1, la 4-aceleraci6n se anula y el escalar de 



^^Tf temperatura de Fermi y T/Tf ^ 10"'' para una ENs. Por lo tanto, el desorden termico no es 
rcsponsable de la presion, ni de la densidad de energia, ni de la magnetizacion del objeto compacto. 

'^Los espacios-tiempo Bianchi inclinados tienen la 4-velocidad "inclinada" , o sea no son sistemas comoviles 
con la velocidad en la direccion temporal. 
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expansion 6 y el tensor de deformaciones toma la forma, 



< = diag[(7^,c7^,a^,0], (2.21) 

A B C , , 

e = ^ + ^ + jj. (2.22) 



donde, 



'^""SA SB 3C" ^"35 3^ 3C" ^ ^ ^ 

_ 2C ^ JB , _ . . . 

"' = 30-31-3^' 

Asi, consideramos como fuente de materia para esta metrica el siguiente tensor de 
energia-momentum, 

T\ = iU + P)u''ub + Pd% + U\, P = p-^^, (2.25) 

aqui n" 5 es el tensor de presiones anisotropicas, U la densidad de energia, P es el valor 
medio dc las tres presiones diagonalcs, B cl campo magnctico y la magnetizacion, a 
todos ellos los suponemos funciones del tiempo. Note que la anisotropi'a es producida por 
el campo magnetico B. Si este campo se anula, o sea: B = 0, entonces el tensor de energi'a- 
momentum se reduce al tensor del fluido perfecto con presion isotropa P = p. En el caso 
general ;B 7^ el tensor 11" 5 tiene la forma, 

n% = diag[n,n,-2n,o], n = -^, u\ = o. (2.26) 

Las ecuaciones de campo de Einstein (ECE) asociada con la metrica de Kasner y con 
el tensor de energia-momentum quedan como, 

-G^ = -gg + - + ^ = -/^(p - BM), (2.27) 
AC' A C 

-^^ = if + 1 + 1 = 

donde el punto denota la derivacion respecto al tiempo propio de observadores fundamen- 
tales y K = SttG/c^. A partir de la ecuacion de balance T°-\ = y de las ecuaciones de 
Maxwell F"^ = y F^ab;c] = 0) tenemos que, 

U+(p + U)e- BMi^ + |) = 0, (2.31) 

A B IB , , 
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Nosotros necesitamos construir un sistema auto-consistente de ecuaciones diferenciales 
de primer orden que puedan ser resueltas numericamente, asf que sen'a muy conveniente 
eliminar primeras y segundas derivadas de las funciones metricas en las ecuaciones de 
Einstein-Maxwell y escribirlas en terminos del escalar de expansion y de las componentes 
del tensor de deformaciones. Siendo consecuentes con esta idea, combinamos las ecuaciones 
(2.21)-(2.22),(2.23)-(2.24), (2.27)-(2.30) y (2.31)-(2.32), eliminamos las funciones A,B,C 
junto con sus derivadas A, A, B, B,C,C. Despues de algunas manipulaciones algebraicas 
arribamos al siguiente vinculo, 

- (S^)2 - S^S^ + — = kU, (2.33) 
3 

y a las siguientes ecuaciones diferenciales, 

U + (U + p)@ - BMi^e - J^') = 0, (2.34) 

ty = -'^BM - S^^e, (2.35) 
2k 

= —BM - E^e, (2.36) 

e = k{BM - ^p) - ^ - lii^y + S^)' - S^'S^), (2.37) 

/3 = ^^(3E^ - 26), con P = B/Bc , (2.38) 

donde = erf que la tomamos como la componente independiente del tensor de deforma- 
ciones^^ . 

Aunque el tensor de deformaciones se puede dctcrminar por esta linica componente, 
para nosotros sera muy conveniente, desde el punto de vista de los calculos numericos, 
utilizar dos componentes de dicho tensor. 

Las funciones U,p y M. ahora vendran dadas por las ecuaciones de estado del gas^^, 

p = XTp{P,fi), BM = XprM{P,fi), U = \T^{P,fi). (2.39) 

2.4 Dinamica de cantidades covariantes 

Consideremos ahora las siguientes variables, 

& d 1 d , ^ , 

H = —, — = , (2.40) 

y las funciones adimensionales, 

s..^,s^.E.,a.^,,n^-l-. (2.4:, 

IlQ tlQ 6rlQ tlQ 

donde y estan relacionados con las componentes yy y zz del tensor de deformaciones, 
mientras que esta relacionado con el campo magnetico. El nuevo tiempo r es un tiempo 



^*Las componentes E"^ y se definen similarmente. 

^^Aquf se sobrentiende que fj, = He potencial qui'mico de los electrones, para el caso de los neutrones 
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adimensional (o tiempo "logaritmico" ) . La cantidad H{t) (debido a (2.21)-(2.22)) cambia 
las dimensiones de cm~^ y el signo dc r sc determina por el signo de H{t) (ver Apendice 
E). Sin embargo, nosotros hemos elcgido Q = (3 por tanto k\ = SHq, para eliminar la 
presencia de constantes adicionales en el sistema de ecuaciones. Insertando la ecuacion de 
estado (2.39) y las nuevas definiciones (2.40), (2.41) en (2.33) y (2.34)-(2.38), obtenemos 
el vinculo, 

-{Syf - {S'f - S^S' + 3H^ = 3Tu, 

mas el sistema, 

s^^ = -pfm - ss^n, 

^'.^/3r.-^(r. + ^^ + i^-^), 

Note que, contrario a los modelos cosmologicos discutidos en la referenda [26] donde 
Hq = 0.59 X 10~^^m~^ podria jugar el rol de la constante de Hubble, en nuestro gas 
magnetizado de Fermi tenemos que Ho = 0.86 x lO^'^'^m"^, el cual es mucho mas pequeiio. 

Esto es logico y consistente porque indica que nuestro modelo simplificado es examinado 
sobre escalas locales mucho menores que las escalas cosmicas. La escala I/Hq ~ 1.15 x 
10^'' m es del orden de la distancia de una unidad astronomica. 

Los resultados mencionados en las referencias [19, 20, 21] para el gas de electrones 
muestran que para un campo magnetico intenso, del orden del campo cn'tico Be, todos 
los electrones estan en cl nivcl basico de Landau n = 0, y consecuentemente se ticnc que 
p± = 0. Consideramos cl cstudio de como cvoluciona un gas en estos casos, en el cual las 
funciones r(/3,/i) se simplifican considerablemente, tomando la siguiente forma, 

Tp = ^(/x^ - ^) + ^arcsinh(^) + P{a, - 6J, (2.48) 

TM = iao-K), (2.49) 
11 a 

Tu = - o) - :jarcsinh(-^) + /3(ao + 6J, donde (2.50) 

Oq = AiV)"^ - 1, K = ln(/x + a^/ix), =0, = 1. 

Asi sustituyendo (2.48)-(2.50) en (2.43)-(2.47) mantenemos el sistema auto-consistente 
de cinco ecuaciones diferenciales ordinarias (2.43)-(2.47), con las funciones incognitas P, %, , 
y /X, y el vinculo (2.42), el cual solo admite una solucion numerica. 

De las ecuaciones de estado (2.39) se deduce que el potencial quimico debe satisfacer 
/X > 1, lo que es correcto para sistemas con densidades del orden de ~ lO^gm/cm'^ o 
mayores. Para EEs o ENs, el potencial quimico toma valores alrededor de \/3 ~ 1.732. 
Tambien U > Q, entonces el potencial quimico > 1 y de (2.42) obtenemos el vinculo, 

-{syf - (5^)2 - gys^ + sn^ = svu > o, (2.51) 



(2.42) 

(2.43) 
(2.44) 

(2.45) 
(2.46) 
(2.47) 
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asi nuestro espacio de fase fi'sico de 5-dimensiones, estara restringido por las relaciones, 



> {Syf + {S'f + S^S', (2.52) 
/i^ > 1 + 2/3, (2.53) 

P > (2.54) 

Para modelos Bianchi-I no-inclinados solo existe una componente independiente del 
tensor de deformacion, la cual tomamos por , sin embargo sera conveniente usar tambien 
la componente para nuestros calculos numericos. Ademas, la forma de las componentes 
del tensor de deformaciones determina la forma de los coeficientes metricos. De (2.21)- 
(2.22) y (2.23)-(2.24), tenemos que, 

^ = iS- + 'H), ^ = {Sy + n), ^ = {S'+H), (2.55) 
donde 5^ = -Sy - S\ 
2.5 Discusion y resultados 

Usando (2.21)-(2.22) y (2.40) podemos expresar el volumen local como V = ABC en 
terminos de y el tiempo adimensional r como 

V(t) = V(0) exp ^3 j ndr^ . (2.56) 

Claramente se puede notar el por que al tiempo r se le llama tiempo "logaritmico" . 
Note que el signo de T-L{t) implica expansion del volumen local, si {H{t) > 0) y colapso si 

{n{T) < 0). 

2.5.1 Singularidades 

Las ecuaciones de estado que estamos considerando estan asociadas con objetos compactos 
a muy alta densidad (al menos ~ rango de evolucion del modelo para bajas 

densidadcs o caso diluido, no es ffsicamcntc interesante para nosotros y no sera tratado. 
Esto significa que solo consideraremos las fases de colapso del modelo. Asf solo necesitamos 
examinar las condiciones iniciales en las cuales la expansion inicial Ho cs ncgativa. Nosotros 
probamos numericamente el modelo usando una gran cantidad de condiciones iniciales 
diferentes, cubriendo todo el rango de valores de interes fi'sico para los objetos compactos, 
desde EBs hasta ENs. Por ejemplo: Ho = 2 corresponde con densidades de ~ 10^° kgm~^, 
mientras que /3o = 10^^ representa campos magneticos del unos 10^ T. Junto con "Hq < 0, 
consideramos en particular: Sq = 0, ±1 y Sq = 0, ±1, que corresponden con casos de 
deformacion inicial nula, y deformaciones iniciales positivas y negativas en las direcciones 
y o z. 

Asi cuando T-Lq < el modelo exhibe un comportamiento en general de colapso H — >■ 

— oo, indcpcndientemente de los valores iniciales de las otras funciones (ver ejemplos en 
la Figura 1). En todas las configuraciones de colapso la intensidad del campo magnctico 
diverge a infinito independientemente de las condiciones iniciales. En la Figura 1 mostramos 
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varias soluciones numericas para difcrcntcs valorcs del campo magnetico inicial. Se puede 
ver como al aumentar la intensidad del campo (etiquetada por /31 < ^2 < /33 cn la 
figura), la expansion muestra un rapido decaimiento a — oo. Consecuentemente, siempre 
que aumentemos el campo magnetico inicial, el tiempo de colapso decrece. Usamos 
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Figura 1: Soluciones numericas para el campo magnetico adimensional (3{t) y para el escalar 
de expansion 'H{t). Aqui disponemos de tres condiciones inicialcs difcrcntcs, 5*^(0) = 0,5^(0) = 
-1,S'^(0) = 1,^(0) = 2, y /3i(0) = 10-^/32(0) = 5 x 10-^/33(0) = lO""^ respectivamente, para 
el campo magnetico inicial. Entonces debido a (2.51), la expansion inicial debe tomar los valores 
7^1 (0) = 'H2(0) = V-siO) = —4.82. Los rcsultados numericos para los tiempos de colapso son 
Ti = 0.1059, T2 = 0.1058, y T3 = 0.1057 respectivamente. 



condiciones iniciales con deformacion inicial no necesariamente en la direccion del campo 
magnetico, asf que las singularidades anisotropas del tipo "cigarro" emergen a lo largo 
de cualquiera de estas direcciones, o sea, el colapso puede ser paralelo o perpendicular al 
campo magnetico. 

Un analisis de todos estos casos revela que los estados de colapso dependen fuertemente 
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Figura 2: Trayectorias de las funciones para i = x,y, z. Nosotros tomamos las siguientes 

condicioncs inicialcs: 5'"^(0) = 1, 5''(0) = 0, 5^(0) = -1, /i(0) = 2, /3(0) = 0, y "^(0) = -3.43. 
Este sistema muestra la deformacion inicial en la direccion x con el campo magnetico apuntando 
en la direccion z. El colapso forma una singularidad tipo "cigarro" o "Imea" en la direccion x. 

de los valores iniciales de la deformacion. Por ejemplo, si tenemos una deformacion inicial a 
lo largo de la direccion x, se dice que : Sq ':$> Sq, Sq, Po, fiQ, entonces una singularidad tipo 
"cigarro" emerge a lo largo de la direccion x. Esto se muestra en la Figura 2, ilustrando 
(por medio del conjunto de ecuacioncs (2.55)) que la funcion metrica A ticndc a infinito, 
mientras que las otras funciones metricas B y C rapidamente decaen a cero. En gen- 
eral, la configuracion inicial del sistema y los valores iniciales del tensor de deformaciones 
determinan la direccion (privilegiada) del colapso anisotropo. 

Las pruebas numericas tambien muestran que existe un valor umbral para el campo 
magnetico inicial y dcpcndiendo de su intcnsidad, este tambien puede influir en la direccion 
que tome la singularidad tipo "cigarro". Esto se puede ilustrar apoyandonos en la propia 
Figura 2: Si aqui aumentamos el campo magnetico a Pq = 1 entonces obtendremos una 
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singularidad tipo "cigarro" pero en la direccion z para cualquier valor de la deformacion 
^^(O) = 0, ±1 (donde i = 1,2,3 direcciones espaciales). Sin embargo, /3o = 1 ~ 10^ T es 
una valor ffsicamente irreal para estrellas EBs magnetizadas, pero puede ser razonable en 
modelos de universos primer diales magnetizados. 

Una anisotropia tipo "punto" siempre puede emerger de configuraciones con defor- 
macion inicial nula, o sea con Sq = Sq = Sq = y = 0. Sin embargo, incluso 
con deformacion inicial nula existen siempre valores umbrales del campo magnetico para 
los cuales las singularidades tipo "punto" , (naturales de configuraciones sin deformacion) , 
pasan a ser singularidades tipo "cigarro" extendidas en la direccion del campo magnetico. 

2.5.2 Espacio de fase y subespacios criticos 

Como hemos mencionado previamente, el sistema de ecuaciones de Einstein-Maxwell puede 
ser escrito como un sistema autonomo asociado a un espacio de fase de 4-dimensiones, con 
las variables (S^,I3, ji^ Ti). Note que siempre puede ser encontrada si determinamos = 
S, la cual hemos supuesto como la componente independiente del tensor de deformaciones. 
Consider ando las 4-funciones mencionadas, tendremos el siguiente sistema de ecuaciones, 

tl = -(U+p-lBM)e-BMJ:, (2.57) 

S = ^kBM - es, (2.58) 

@ = k{BM + ^{U -p))-®"^, (2.59) 

P = ^^(3S - 29), (2.60) 

donde S = a^. Si pasamos a trabajar con las funciones adimensionales, definidas en (2.41) 
este sistema se convierte en, 

= 2/3r^ - 3HS^ (2.61) 

H,r = prM + l{ru-rp)-3n^ (2.62) 

= - 2H), (2.63) 

^^,r = T^[(2H - s'){rM + 2ru,f3)(3 - sniVp + Vu)]. (2.64) 

i-u,n 

Note que solo se modifica la ecuacion para 'H(r), ya que podemos Uegar a la ecuacion 
(2.62) usando (2.45) y el vinculo (2.42), por consiguiente ambos sistemas son equivalentes. 
Entonces, la ecuacion (2.43) es necesaria solamente para calcular los coeficientes metricos 
restantes. 

En la Figura 8 del Apendices F, representamos una seccion 3-dimensional del espacio de 
fase {S^, P, n) con diferentes curvas para distintas condiciones iniciales. Como se muestra en 
la figura, el valor inicial de la expansion Hq determina la evolucion global de las soluciones 
numericas. Todas las curvas que inician en r = con expansion Hq = kA/3, convergen al 
atractor estable "a" . Pero en el caso que tomcmos Hq = — \J kA/3 y empezando en r = 0, 
entonces las curvas evolucionaran tendiendo a la singularidad anisotropica. Tomando el 
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lado izquierdo del sistema de ecuaciones en (2.61)-(2.64), igualandolo a cero y resolvicndolo 
de forma algebraica, encontramos el conjunto de puntos cn'ticos. Dentro de este conjunto 
de puntos cn'ticos se encuentra el atractor estable "a" , 

a = {5^ = 0,/3 = 0,/x = l,H = 0}. (2.65) 

Para el espacio de 4-dimensiones {S^, P, ji, %) existen cuatro posibles secciones 3-dimensionales. 
En estas secciones hemos calculado tambien las soluciones numericas y los resultados son 
similares. 

2.6 Conclusiones del capftulo 

Hemos presentado un modelo basado en la descripcion dinamica de un volumen local de 
un gas magnetizado y auto-gravitante de electrones en el nivel basico de Landau n = 0. 
Consideramos las ecuaciones de estado de este sistema. Tambien, hemos trabajado con 
formas simplificadas de las ecuaciones de Einstein-Maxwell que fueron escritas asumiendo 
un espacio-tiempo Bianchi-I con metrica de Kasner. La fuente de anisotropia es el propio 
campo magnetico. Este espacio-tiempo simplificado provee un convenientc modelo juguctc 
para entender el comportamiento local de colapso del sistema que cualitativamente imita 
las fuentes de materia dentro de un objeto compacto. 

La relevancia de este capitulo emerge de nuestros estudios sobre clasificacion de sin- 
gularidades, las cuales pueden ser de tipo "punto" o de tipo "cigarro" . Las singularidades 
tipo "punto" apareccn bajo condiciones especiales del campo magnetico nulo, deformacion 
nula o ambas. Singularidades tipo "cigarro" pueden ser obtcnidas en todas las direcciones, 
dependiendo de los valores iniciales del tensor de deformaciones. Sin embargo, para in- 
tensidades suficientemente grandes del campo magnetico inicial siempre se cspcrara una 
singularidad tipo "cigarro" en la direccion del campo magnetico. 

Este resultado es importante porque el valor del campo magnetico determina el tipo 
de colapso y esto esta de acuerdo con la idea previa no-relativista en la cual se analiza el 
colapso de estos gases magnetizados dentro de un cuadro Newtoniano[19],[20], [21]. 

Por otro lado, asi como es discutido en [44] por Collins k, Ellis, los modelos Bianchi 
como el de tipo-I que estamos considerando, son globalmente hiperbolicos y solo presentan 
una singularidad, asi que las hipersuperficies de tiempo constante (ortogonales a la 4- 
velocidad) son hipersuperficies globales de Cauchy^'^ y cada punto en el espacio-tiempo 
puede ser causalmente conectado con el ultimo. La 4-velocidad es una geodesica del campo 
y la singularidad es determinada por un valor especifico de tiempo constante, asi que 
esta singularidad no es de tipo temporaloidea^^ y cada evento en el espacio puede ser 

•^"Una superficie de Cauchy es una variedad 3-dimensional espacialoidea, cuyo future y pasado de depen- 
dencia describen toda la variedad espacio-tiempo. 

^^Segiin su caracter las singularidades fisicas pueden ser: 1) Singularidades temporaloideas, como la 
que se encuentra en un agujero dc Schwarzschild en la que una partfcula deja de existir por cierto instante 
de tiempo; dependiendo de su velocidad, las partfculas rapidas tardan mas en alcanzar la singularidad 
mientras que las mas lentas desaparecen antes. Este tipo de singularidades son inevitables, ya que tarde 
o temprano todas las partfculas deben atravesar la hipersuperficie temporal singular. 2) Singularidades 
espacialoideas, como la que se encuentra en agujeros de Reissner-Nordstrom, Kerr y Kerr-Newman. Al 
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causalmente concctado con esta singularidad (en particular por geodesicas que son curvas 
integrales del campo de la 4-velocidad) . 

Nosotros solo hemos considerado regimenes de colapso a partir de una hipersuperfi- 
cie inicial de tiempo constante. Asi, cada curva temporaloidea futura (geodesica o no) 
comienza en una hipersuperficie de Cauchy de tiempo constante y termina en la singulari- 
dad de colapso. Bajo estas condiciones, la singularidad es obviamente censurada. 



ser hipersuperficies espaciales una partfcula puede escapar de ellas y por tanto se trata de singularidades 
evitables. 



3. ESTUDIO DE UN GAS MAGNETIZADO Y AUTO-GRAVITANTE 
DE NEUTRONES 



3.1 Introduccion 

En este capi'tulo extendemos el trabajo previo hecho para un gas de electrones, al caso de 
neutrones. Seguimos una similar metodologfa basada en la reescritura de las ecuaciones 
de campo y de las ecuaciones de conservacion, para el gas de neutrones en una geometria 
Bianchi-I. Como sistema dinamico, estudiamos su evolucion en un espacio de fase de 4- 
dimensiones. Este sistema es analizado cualitativamente y numericamente. 

Vale la pena resaltar las diferencias basicas entre el gas magnetizado de electrones 
examinado en [40, 22, 52] (ver capitulo anterior) y el gas de neutrones que nosotros consid- 
eraremos en este capitulo. Los electrones interactuan con el campo magnetico a traves de 
su carga electrica, contribuycndo al llamado diamagnctismo de Landau caractcrizado por 
un efecto de cuantizacion asociado con los niveles de Landau. Sin embargo, los neutrones 
interactuan con el campo magnetico mediante el Momento Magnetico Anomalo (MMA), 
en el contexto del paramagnetismo de Pauli y de las ecuaciones de Pauli-Dirac. Conse- 
cuentemente uno espera que la interaccion magnctica del neutron sea debil, aunque un 
gas de neutrones degenerado en condiciones de objeto compacto se espera que sea mas 
cn'tico (dcbido a las altas densidades) que las de un gas degenerado en condiciones de baja 
densidad. Por tanto, los efectos relativistas de la gravedad tendran que tomarse en cuenta 
y podrian dominar desde el punto de vista dinamico en el gas de neutrones. Tambien es 
importante mencionar que un gas de neutrones auto-gravitante es un modelo simplificado 
dc una fuente para un objeto compacto, donde los protones y cl potcncial qufmico de equi- 
librio deben ser incluidos, para evaluar los efectos locales [40, 50, 51]. Sin embargo, el gas 
de neutrones magnetizado [45] ya exhibe importantes diferencias, tanto cualitativas como 
cuantitativas, en comparacion con el gas de electrones previamente analizado. 

3.2 Ecuacion de Estado de un gas magnetizado de neutrones 

Propiedades importantes del gas de neutrones magnetizado y degenerado son reportadas 
en (ver [40, 20, 53]). Considerando el ensemble gran canonico, un subsistema puede ser 
pensado como un volumen local de un gas de neutrones bajo la influencia de un campo 
magnetico H asociado al resto del sistema (en un contexto astrofi'sico esta puede ser una 
buena aproximacion a un volumen local dentro de un objeto compacto). Debido a este 
campo los subsistemas se pueden polarizar, permitiendo que surja una magnetizacion que 
satisface la relacion: H = B — AttM. El campo H puede ser pensado como un campo 
"externo" al subsistema, mientras que B sen'a el "interno" de cualquier particula dentro 
del subsistema, el cual incluye (en adicion a H) la contribucion AnAi de partfculas del 
resto del subsistema. 

La ecuacion de estado para el gas de neutrones se obtiene de calcular el espectro de 
energi'a de las partfculas que forman el sistema. Nosotros podemos obtener este espectro de 
energfa a partir de las ecuaciones de Dirac para partfculas neutras con momento magnetico 
anomalo. 
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donde a^x = ^(7^ 7a — 7a 7ft) es el tensor de espm, F'^'^ es el tensor del campo electro- 
magnetico (nosotros hemos tornado /i = c = 1) y ^' cs cl campo dc Dirac. Resolviendo la 
ecuacion (3.1) se obtiene el siguiente espectro de encrgi'a [53], [54], [55], 



En{p,B,ri) 



P± + ml + riqBy, 



(3.2) 



donde ppp± son respectivamente, las componentes del momentum en la direccion paralela 
y perpendicular al campo magnetico B, rUn es la masa del neutron, q = — 1.91//Ar es el 
momento magnetico del neutron (/Ltjv = e/2mp es el magneton nuclear), 77 = ±1 son los 
autovalores correspondicntcs a las dos orientaciones posibles (paralela y antiparalela) del 
Momento Magnetico Anomalo (MMA) del neutron respecto al campo magnetico. 

Sustituyendo (3.2) en la expresion (2.6) obtcnida en el capftulo anterior podemos cal- 
cular todas las cantidades termodinamicas del sistema al igual que hicimos en ese mismo 
capftulo para el gas de electrones magnetizado. 

Tal como explicamos en dicho capftulo el termino de vacfo es divergente, y puede 
igualmente ser renormalizado, y para cl gas de neutrones su contribucion es relevante solo 
para campos de intensidad B > 10^^ T [20], por tanto nosotros no tomaremos en cuenta 
este termino^^. 

La ecuacion (2.6) puede ser integrada facilmente^^ para el caso degenerado (T = 0), y 
su forma explfcita es entonces, 



v=i,-i 



jif^ (l + 77^)(57?/3-3)//^ (1 + 77/3)^^(3 
12 24 24 



donde hemos introducido las siguientes expresiones, 

' + fv 



TT 

2 



arcsm 



1^ 



rur, 



In 



(3.3) 

(3.4) 
(3.5) 



.l + ??/3, 

con Be = run/q — 1.56 x 10^^ T el campo crftico para neutrones y A 
4.11 X 10^5 J m-3. 

Calculando la densidad de neutrones y la magnetizacion obtenemos que: = Aor^v, 
M. = MoTm, donde Nq = X/rUn, Mq = Noq, y los coeficientes Tn,Tm toman la forma. 



3 2 2 



M 



= - E ^ 



r?=l,-l ^ 



(1-277/3)///^ {l + vP)Hl-vP/2) ^ 

6 3 + 6 " 



"^"^Rcalmente el termino dc vaci'o pudiera ser relevante para cuando el sistema evoluciona hacia la sin- 
gularidad y el campo magnetico crece a infinito pero, por simplicidad no lo tendremos en cuenta en este 
trabajo y suponemos que no influye de manera determinante. 

^^Recordar que en este caso de neutrones en (2.6) no se suma por los niveles de Landau sino por —1, +1, 
los dos estados posibles del MMA. 
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Por consiguiente, dados (2.9)-(2.11) y (3.3), nosotros podemos escribir la ecuacion de es- 
tado para un gas de neutrones relativista y degenerado en presencia de un campo magnetico 
como, 

U = finN + n = \Tu(fi,n), (3.6) 

p = -n = xrpij3,fi), (3.7) 

M = BM = XPTm{P, h), (3.8) 



donde, 



j?=l,-l 



fif^ (l + r?/3)(5r?/j-3)/x/^ 
12 24 



(l + ?7^)^(3-?7^ rjl3fi' 
+ 24 ~ ~6~' 

= I^^N — Tp- 

Note que estas ecuacioncs de estado difieren de las obtenidas para el gas de electrones 
(2.12)-(2.13) y obviamcntc los rcsultados dinamicos estaran marcados por estas diferencias. 
Las ecuaciones (2.12)-(2.13) dependen de los niveles de Landau por estar cuantificado el 
espectro. Las ecuaciones (3.6)-(3.8) muestran su acoplamiento con el campo a traves de la 
presencia del MM A de los neutrones. 

Puntualicemos que en (3.3)- (3.8) nosotros estamos sumando sobre el momento magnetico 
paralelo (ry = — 1) o perpendicular (r/ = 1) al campo magnetico (o sea estamos considerando 
el conocido Paramagnetismo de Pauli). La eleccion der] = ±1, [20] es equivalente a consid- 
erar diferentes fases del sistema. La aparicion de valores umbrales para el campo magnetico 
en cada caso pueden verse si analizamos las expresiones de las funciones y en (3.4). 
Si tomamos a ^ > 0, empezando con 

> (l + r//3)2 = 1 + 277/3 + 772/32. (3.9) 
Sin embargo, como 77 = — 1, 1 entonces siempre r;^ = 1 y nosotros podemos reescribir, 

/x'>(r? + /S)2, (3.10) 

de donde, 

{p + n-fi){p + n + fi)<o. (3.11) 

Ahora tenemos dos posibilidades. Pero, es simple darse cuenta que la linica aceptable es; 

-H-ri< P <n-v- (3-12) 

Esta es exactamente la restriccion que vicnc dc la funcion, s^, 

I |< 1. (3.13) 



-41 - 



Por tanto, ya sea de (3.12) o de (3.13) los siguientes vmculos son obtenidos para el campo 
magnetico, 



Si r? = l^ -1-At</0<A*-1, (3.14) 
Si n = -l^ l-H</3<l + n. (3.15) 

Estas ultimas desigualdades restringen los valores del campo magnetico. Asi, sistemas 
de neutrones con MMA alineado en la direccion del campo magnetico solo tienen como 
valores permitidos del campo magnetico el rango dc 1 — fi < P < 1 + /i. Similarmente, para 
sistemas de neutrones con MMA orientado antiparalclamente respecto al campo magnetico, 
solo pueden tomar valores de campo en el intervalo —1 — n < /3 < /i — 1. En particular 
para el caso de = 1 tenemos que, 

Si = -1, < /3 < 2, (3.16) 
Si r? = 1, -2 < /3 < 0. (3.17) 

Esto significa que para campos menores (o iguales) que dos veces el campo magnetico 
cn'tico, los neutrones estaran alineados todos con rj = —1 (el MMA es paralelo al campo 
magnetico). O sea neutrones con t] = 1 seran forzados a invertir su sentido. 

3.3 Ecuaciones de Einstein— Maxwell 

Si deseamos estudiar la evolucion dinamica de un volumen local magnetizado y auto- 
gravitante de neutrones bajo las condiciones fi'sicas impuestas en el interior de un objeto 

compacto, entonces tenemos que contar con los efectos relativistas. Y esto implica que la 
dinamica local debe ser estudiada bajo el marco de la TGR, por medio de las ecuaciones 
de campo de Einstein, 

Gab = Rab - -Rgab = KTab, (3-18) 

junto con la ecuacion de conservacion del tensor energi'a-momentum y de las ecuaciones de 
Maxwell, 

r''^6 = 0, (3.19) 
^"';6 = 0, F^ab;c]=0, (3-20) 

dondc K = SttGn , con Gn constante de gravitacion de Newton, el corchcte cuadrado denota 
anti-simctrizacion en ab; c. El tensor energia-momentum T"^ h asociado al gas magneti- 
zado de neutrones viene dado por (2.9)-(2.11), conectado al gran potencial termodinamico 
obtenido mediante la mecanica estadi'stica, y acoplado apropiadamente a la ecuacion de 
estado. Este tensor tambien puede ser escrito en terminos de la 4-velocidad del campo 



24 

como , 



T\ = {U + P)u''ub + PS\ + U\, p = p-'^^. (3.21) 



^''Notar que el tensor (3.21), no contiene, ni el termino de vaci'o log{B/Bc), ni se adiciona el campo 
clasico B^ (de arbol). Ya que suponemos que la contribucion del vaci'o es pequena y similarmente a como 

hicimos en el capi'tulo 2 hemes seguido trabajos en donde este termino clasico no se toma en cuenta. Aunque 
debido a su importancia factible en la evolucion cercana al colapso, sf los asumiremos en el capftulo 4. 



-42 - 



Nosotros consideramos las ecuaciones de campo (3.18)-(3.20) con (3.21) como la fucntc 
de materia para un modelo Bianchi-I descrito en la representacion de la metrica de Kasner, 



ds"^ = Qiitfdx'^ + Q2{tfdy'^ + Q^(tfdz'^ - df, (3.22) 

lo cual sugiere escoger la 4-velocidad en una representacion comovil = (5" j, asi el tensor 
de presiones anisotropas en (3.21) y en las coordenadas [x,y,z,t] toma la forma, 

n% = diag[n,n,-2n,o], n = -^, u\ = o. (3.23) 

Las ecuaciones de campo (3.18) para (3. 21), (3. 22) y (3.23) toman la forma, 

-G% = §^ + 1^ + 1^ = -.(p - BM), (3.24) 

-G\ = §^ + ^ + ^ = ^U, (3.27) 

V1V3 V2V3 

donde A = A.^ u'^ = A^f De la ecuacion de conservacion del tensor energia-momentum 
(3.19) tenemos que, 

U = ^{p + U)-{^ + ^)i-BM +P + U), (3.28) 

V3 Wl W2 

mientras que de las ecuaciones de Maxwell (3.20) obtenemos que. 

Las ecuaciones de Einstein-Maxwell (3.24)-(3.27),(3.28) y (3.29) son ecuaciones difer- 
enciales ordinarias no-lineales y de segundo orden para las funciones metrica Qi,Q2,Q3 y 
U. Para tratar este sistema numericamente, es necesario introducir nuevas variables que 
transformaran el sistema, en un sistema de ecuaciones de primer orden. Sin embargo, antes 
de tomar este camino, veremos que ocurre cuando el campo magnetico es debil. 

3.4 Limite de campo magnetico debil 

Esta breve discusion para el campo magnetico debil es importantc pues ilustra la conexion 
entre nuestro campo magnetico de origen cuantico y el campo magnetico Maxweliano en 
el contexto de un tratamiento magnetohidrodinamico. En la seccion 3.2 nosotros hemos 
mostrado ecuaciones fuertemente conectadas con un campo magnetico cuantico. 

Es importante enfatizar que el tcrmino "campo magnetico cuantico" contiene la in- 
teraccion semi-clasica entre el campo magnetico y el momento magnetico anomalo. Esta 
aproximacion implica una conexion teorica entre la ecuacion de estado introducida en la 
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seccion anterior y el cuadro dc trabajo en EDC . Como una condicion de consistencia, 
este marco de trabajo debe permitir en el caso Kmite de una expansion en series alrede- 
dor de (5 = 0, Uegar al Iimite clasico Maxweliano. Los terminos Iideres en esta expansion 
deben coincidir con los terminos del tensor de energia-momentum para el campo magnetico 
Maxwelliano [26] . En general, este tipo de expansion en series puede efectuarse de la forma, 



p± 


QQn > 

n=0 






(3.30) 


p\\ 








(3.31) 


u 


n=0 


TIT 




(3.32) 



donde es facil ver que: pi = p||(/3 = 0) = p±{l3 = 0), oq = {d^U/dp"^) l^^^ /2, ai = 
{d'^p±/d/3'^) l^^g /2, as = {d'^p\\/d(5'^) \^^^ /2, Uq = U{j3 = 0), y todas estas funciones 
dependen solamente del potencial quimico adimensional ji. 

Como se muestra en los ejemplos encontrados en la literatura, los campos magneticos 
clasicos sobre un fondo geometrico tipo Bianchi-I [57, 32], los terminos lideres son terminos 
cuadraticos. Esto sugiere que las series en (3.30)-(3.32) deben ser truncadas a nivel 
cuadratico, los terminos de mayor orden como son en general contribuciones 

multipolares muy pequenas. Si asumimos que las fluctuaciones de las velocidades del 
plasma tienden a cero como promedio macroscopico, y que el medio no sufre de movimien- 
tos de volumen internos, entonces estas contribuciones pueden ser despreciadas (aunque 
tipicamente, velocidades altas podrian surgir a partir de fluctuaciones termicas o desor- 
denes cuanticos). Bajo estas suposiciones, el tensor de energia-momentum de un gas de 
neutrones mmimamente acoplado al campo magnetico, puede siempre ser escrito en la 
forma [57], 

Tab = {Uo + Umag)UaUb + (po + Pmag)hab + H^"^, (3.33) 

donde Pmag y Ujnag SOU, respectivamente, la presion magnetica y la densidad de energia del 
campo magnetico. 

Es importante tener en cuenta que po en (3.33) es la contribucion isotropica a la 
presion del sistema, la cual (cn general) puede dcpcndcr del potencial quimico y del campo 
magnetico. Por otro lado, pi es la presion para el caso de campo magnetico nulo ^ = 0. 
En nuestro caso el tensor de materia tiene la forma, 

b = diag[pi - aiP^pi - aiP\pi + a^P^ -Uq - aoP% (3.34) 

— = ao^Vs, (3.35) 
^ = ao^', (3.36) 

Electrodinamica Cuantica 
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y comparando (3.34)-(3.33) tenemos que. 



Pmag 



Ufnag 



dondc = H^Hi = 2ao/3^ y Hi son las componentes del campo magnetico, la cuales hemos 
asumido que apuntan en la direccion z. Esta suposicion es consistente con el hecho de que 
un pequeno volumen en el centre de un objeto compacto es aproximadamente homogeneo, 
y la rotacion puede efectuarse en una direccion privilegiada, que nosotros tomamos como 
la direccion del eje z, 

W = (0, 0, V2^^, 0), Hi = (0, 0, V2^/3Q3, 0). (3.37) 

El tensor H^"^ es entonces el tensor simetrico de traza nula, proyectado con hab, y 
representa las presiones anisotropas debidas al campo magnetico. Este tensor puede ser 
escrito como, 

^-^mag^a ^ ^ diag[-i(ai + a3)/3^ -^K + os)/?', ^(ai + a3)/3^ 0]. (3.38) 
Es importante puntualizar que pi en (3.34) es en general, diferente de po- Asi, 

Po = Pi - (ao + 2ai - aa) — , (3.39) 

donde pi = P(/3=o) = -f(/3=o) (3-7) y (3.21) es la presion sin campo magnetico. Solo 
cuando el campo magnetico se anula /3 = 0, las presiones coinciden: po = Pi y el tensor de 
energi'a-momentum toma la forma del tensor del fiuido perfecto con presion isotropica. En 
este caso po y Uo corresponden a la presion y la densidad de energia de un gas clasico de 
neutrones. 

3.5 Dinamica de cantidades covariantes 

Como estamos interesados en la evolucion de un elemento dc volumen para el gas de 
neutrones asociado con la fuente (3.21), similarmente a como se hizo para el caso de elec- 
trones debemos reescribir las ecuaciones dinamicas de Einstein-Maxwell en terminos de 
parametros covariantes asociados con la cinematica local de un elemento de volumen como 
es descrito por un observador con 4-velocidad Para la mctrica de Kasner en una 
representacion comovil dc geodesicas, los parametros locales cinematicos no nulos son la 
expansion @, y la deformacion a"''', dada por. 



e = u«;„, (3.40) 

e 

3 



<^ab = U{a;b) " " Kb, (3-41) 



donde hab = UaUb + Qab es el tensor de proyeccion y los parentesis denotan simetrizacion 
respccto a los indices a, b. El escalar de expansion y las componentes del tensor de defor- 
macion para la metrica de Kasner tiene la forma, 

Vl ^2 <^3 

a"^ b = diag [a^ ^, y, a' „ 0] = diag [Si, E2, S3, 0], (3.43) 
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donde, 



^a = ^-§^-§^, « 7^ /? 7^ 7, («, ^, 7 = 1, 2, 3). (3.44) 

oL^a oi^l3 ^^■y 



El tensor de deformacion a" i, es de traza nula: a" a = 0, entonces podemos eliminar 
una de las cantidades (Si, S2, S3) y escribirlas en terminos de las otras dos. En efecto, 
para el Bianchi-I con metrica Kasner, una de estas cantidades es suficiente para representar 
a cr"a, aunque para nosotros es conveniente mantener dos de estas variables, eliminar Si 
y escribirla en terminos de (S2, S3). Por medio de (3.42) y (3.44), todas las dcrivadas 
de segundo orden de las funciones metricas en (3.24)-(3.27), (3.28) y (3.29) pueden ser 
reescritas como derivadas de primer orden de 0, S2 y S3. Despues de varias manipulaciones 
podemos reescribir (3.24)-(3.27), (3.28) y (3.29) como un sistema de ecuaciones de evolucion 
de primer orden, de la siguiente forma, 

if = -(U + p~'^BM)e-BM^3, (3.45) 

52 = -^-eS2, (3.46) 

53 = ^kBM - eSs, (3.47) 
G = k{BM + ^iU-p))- e^ (3.48) 
P = ^/?(3S3 - 26), (3.49) 

junto con estas ecuaciones, tenemos el siguiente vinculo, 

- si - S2S3 + _ - Ei = (3.50) 

que proviene de (3.27). Este sistema de ecuaciones de primer orden de las variables 
U, /3, 0, S2, S3 y el vinculo (3.50) estan completamente determinados^^. La solucion de 
este sistema describe la evolucion dinamica de un volumen local en un gas de neutrones 
magnetizado y representa un modelo aproximado a un subsistema gran canonico para esta 
fuente en el centro de un objeto compacto. 

3.5.1 Ecuaciones dinamicas 

El sistema de ecuaciones de evolucion (3.45)-(3.49) y (3.50) puede ser transformado en un 
sistema dinamico, de la misma forma a como se hizo en el capitulo anterior (ver seccion 
2.4), introduciendo las siguientes variables adimensionales, 

rr ® d 1 d 



^®Notar que a diferencia del capi'tulo anterior aquf, escribimos los resultados de una forma generica, o sea 
A, B, C ahora son Qi,Q2, Qs esto nos resultara muy conveniente. 

^^Note que el sistema (3.45)-(3.49), junto con el vinculo (3.50) es identico al sistema (2.34)-(2.38), 
cambiando E^,E^ por E2,E3. 
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junto con las funciones adimensionales, 

^ = |-' ^2 = 1^, S3 = ^, P = ^, (3.52) 

-Ho -no -"0 Oc 

donde Hq es una constante inversa a la longitud de escala, la cual hemos elegido por 
conveniencia como 3Hq = kA =^ |-ffo| = 1.66 x 10~^m~^ La constante Hq, provee una 
longitud de escala de l/.ffo ~ 6 x lO^m que es adecuada con la longitud caracteristica del 
sistema que estamos considerando. Las funciones ^2 y S3 son las componentes del tensor 
de deformaciones normalizadas a esta escala, el tiempo adimensional r puede ser negative 
o positivo, dependiendo del signo de Hq = ity^ kA/3. 

Sustituycndo las variables (3.52) en el sistema (3.45)-(3.49) y en el vinculo (3.50) 
obtenemos un sistema identico al sistema (2.42) y (2.43)-(2.47)^^ , por lo que no vemos 
necesidad de reescribirlo. 

3.6 Soluciones numericas y discusion fisica 

Estamos interesados en estudiar el colapso de una configuracion dc un gas de neutrones 
magnctizados. Necesitamos considerar un elemento de volumcn local y resolver el vinculo 
(2.42). De las dos soluciones que aporta la ecuacion cuadratica del vinculo para 7i, la 
raiz negativa es la que esta relacionada con soluciones de colapso. De modo que para 
asegurarnos del colapso local debemos imponer en las soluciones del sistema (2. 43)- (2. 47), 
que la expansion inicial G, sea negativa lo cual implica que en la hipersuperficie de tiempo 
inicial r = que 'H(O) < 0. Por tanto, a partir de (3.42) y (3.52) expresadas en terminos 
del volumen local V = ^/detg^ = Q1Q2Q3 como. 



V = V(0) exp (^3 J ndrj . (3.53) 

Para invcstigar el colapso en dependencia de las direcciones en que se puede manifestar, 
podemos relacionar por medio de (3.42), (3.44), (3.51) y (3.52) a las componentes espaciales 
de la metrica y a la combinacion ^ + 5"^, ya que, 

Q0{T) = Qp{O)eMji'H + S0)dT], P = 1,2,3, (3.54) 

donde Q/3{0), son constantcs que pueden ser identificadas con los valores iniciales de Qp{t). 
Para resolver el sistema (2.43)-(2.47) usamos una gran cantidad de condiciones iniciales 
asociadas a las condiciones tipicas que existen en ENs [40, 46, 47, 48, 49, 4, 50, 51]. Por 
ejemplo: = 2 =^ p ~ lO^^kg/m^, Pq = 10"^ — 10~^, para campos magneticos entre 
10^^ T y 10^^ T. Ademas, debemos imponer en todas las pruebas numericas la condicion 
de colapso para el volumen: ^(0) < 0, junto con: ^2(0) = 0, ±1, 6*3(0) = 0, ±1, la cual 
corresponde con casos de deformacion inicial nula y deformacion inicial no nula en las 
direcciones y o z. 



^®Note que Ho, no es la constante de Hubble, dada por Ho"""^ = 0.59 x 10~^'^m"^ 

^^Claramente debemos cambiar en (2.42) y (2.43)-(2.47), S'",S'^ por 82,83 respectivamente. Ademas, 
de recordar que U se reemplazo por fi debido a que de (3.45) tenemos que: U = U{p, jj.) =^ Z7,r = 
A(r[/,|u/Li,T- + Tu,i3l3,r), lo cual nos permite obtener /x,r a partir de U^t- 
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La solucion numerica para la funcion %, cs ploteada en la Figura 3, para diferentes 
condiciones iniciales, muestra que % — )■ — oo, independientemente de las condiciones ini- 
ciales. El campo magnetico tiende ha aumentar, pero siempre se mantiene por debajo del 
campo critico B = Be- Este comportamiento es mostrado en la Figura 3, para todo un 
rango de condiciones iniciales. 
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Figura 3: En el grafico de la izquierda se muestra el comportamiento de vs t para diferentes condiciones 
iniciales. En el grafico de la derecha se muestra la intensidad del campo magnetico adimensional {P = B/Bc), 
este tiene una tendencia a crecer, pero siempre se mantiene por debajo del campo cn'tico. 

Es evidente concluir de estas graficas que las cantidades Si + H ^ — oo, asi como que 
los coeficientes metricos tiende a cero {Qi,Q2,Q3 — > 0), lo cual muestra claramente que el 
elemento de volumen colapsa a una singularidad isotropica tipo "punto" . 

3.6.1 Espacio de fase 

Del mismo modo a como se hizo en el capitulo 2, podemos usar el vinculo (2.42) para 
transformar el sistema (2.43)-(2.47) a un sistema de ecuaciones en funcion de las variables 
{S3, P, IjIjH). Si hacemos esto llegaremos al sistema (2.61)-(2.64) de la seccion 2.5.2, donde 
solo tenemos que intercambiar por ^3. 

Las trayectorias en la subseccion 3-dimensional del espacio de fase, parametrizada 
por (S'3,/3,^), son mostradas en la Figura 9, del Apcndice F. La evolucion del sistema 
es determinada por el signo de Hq. Para r < Hq = — -^/kA/3, y ^ > el sistema 
evoluciona tendiendo al atractor estable "a", mientras que si r > i^o = v^kA/3, y 
H < las trayectorias evolucionan hacia una singularidad "puntual" . Un estudio similar 
fue hecho en las restantes subsecciones 3-dimensionales del espacio de fase, obteniendo 
resultados similares. 
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Figura 4: Es ploteado el comportamiento de {Si +"}{) y {S2 + W) versus r. En estos graficos se puede 
apreciar como estas cantidades tienden a —00, ademas de tener diferentes tiempos de colapso para diferentes 
condiciones iniciales. 
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3.7 Conclusiones del capitulo 

Hemos usado en el modelo Bianchi-I para 
estudiar la evolucion de un elemento de vol- 
umen de un gas magnetizado de neutrones, 
caracterizado fisicamente por una ecuacion 
de estado completamente relativista. Por 
otro lado, en nuestros estudios teoricos ar- 
gumentamos la simplicidad de la geometn'a 
Bianclii como una aproximacion a un subsis- 
tema gran canonico formado por una fuente 
magnetizada de neutrones y sometida a las 
condiciones existentes en el centro de un ob- 
jeto compacto. Por este motivo, nuestro 
tratamiento debe ser concebido como una 
primera aproximacion al problcma, ya que 
cerca del centro del objeto compacto se es- 
pera que las metricas tengan un compor- 
tamiento aproximado, y esta descripcion seria 
litil para entender la evolucion local de elementos de volumen bajo estas condiciones. 

Sin embargo, un apropiado examen sobre las limitaciones especificas dc la dinamica 
de este modelo de juguetc y su conexion concreta con los estudios astroffsicos actuales de 
los objetos compactos van mas alia del alcance de este capi'tulo. En un futuro cercano, se 



a= (S2(0)=-l, 63(0) =1,JK0) =2, P(0) =0,01, Jf(0) =-2,46 , ,^ 
a= (82(0)= -I, 83(0) =-1,(1(0) =2, p(0) =0,001, Jf(0) =-2,59 * 
a= (Sj(0)= 0, Sj(0) =0, n(0) =2, p(0) =0,0001, Jf (0) ^2,39 
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Figura 5: El grafico muestra el comportamiento 
de {S3 + T-L) versus t. Este comportamiento es simi- 
lar a las trayectorias cn la Fig 4, por consiguiente la 
cantidad {S3 +'H) ^ —00. 
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trabajara con metodos perturbativos, metodos numericos mas elaborados y fuentes menos 
idealizadas. 

Las ecuaciones de campo de Einstein-Maxwell para un gas magnetizado de neutrones 
en una geometria Bianchi-I fueron transformadas en un sistema de ecuaciones no-lineales 
de evolucion, las cuales fueron resueltas numericamente para una variedad de condiciones 
iniciales de colapso. Estas soluciones fueron analizadas cualitativamente como un sistema 
dinamico propio. Asi los resultados que encontramos fueron: 

• El estado final en la evolucion de colapso de un elemento de volumen del sistema 

es una singularidad tipo "punto". Este estado final ocurre para un ancho rango 
de condiciones iniciales asociadas con valores que tfpicamente existen en los objetos 
compactos. 

• El campo magnctico aumenta rapidamente a medida que el sistema evoluciona ha- 
cia la singularidad. Aunque su intensidad siempre se mantiene por debajo del campo 
magnetico critico. Este resultado es consistente con los valores obtenidos numericamente 
para el campo magnetico maximo en modelos con configuraciones magnetizadas y en 
rotacion estable (ver [46]). 

• El estudio del espacio de fase asociado con las ecuaciones dinamicas muestran que el 
sistema evoluciona, para H > 0, aun punto de equilibrio, (o sea, a una configuracion 

estable). Es posible introducir una dependencia con la temperatura en la ecuacion 
dc estado. Al introducir la temperatura, la evolucion del gas de neutrones podrfa 
asociarse con fuentes de neutrones a altas temperaturas en el contexto de universos 
magnetizados primordiales en modelos cosmologicos [58]. 

Es importante tener en cuenta que, al contrario del estudio dinamico para un gas de 
electrones en el capftulo anterior (ver tambien [52]), las singularidades anisotropas tipo 
"cigarro" no ocurren para todas las condiciones iniciales impuestas. Sin embargo, el efecto 
dinamico del campo magnetico bajo condiciones crfticas es un estado final de singularidad 
anisotropa alineada en la direccion del campo magnetico. La exclusiva aparicion de sin- 
gularidades tipo "punto" para el gas dc neutrones magnetizado sugiere que cl estado final 
de evolucion de este gas es dominado mas intensamente por la fuerte gravedad en com- 
paracion con el gas de electrones. Esto es consistente con el hecho que los electrones estan 
fuertemente acoplados al campo magnetico mediante su carga electrica, sin embargo los 
neutrones tiene un acoplamiento mas debil asociado con su momento magnetico anomalo. 

Los modelos con geometria Bianchi I en una representacion comovil tiene una severa 
limitacion relacionada con los efectos dinamicos asociados con el campo magnetico. Esto 
es importante cuando se consideran a los neutrones como fuente para los cuales la carga 
electrica es cero, no asi su momento magnetico anomalo. 

Siendo estos espacio-tiempos espacialmente homogeneos con 4-velocidad ortonormal a 
las hipersuperficies de simetrfa maximal, la fuerza de Lorentz es cero: = qUaF"'^ = 0. 
Tambien, en efecto los modelos Bianchi-I son espacialmente pianos y esto hacc al modelo 
inadecuado para estudiar (incluso como modelo juguete) el colapso local de dicho volumen 
elemental magnetizado dentro de objetos localizados. O sea en tales objetos compactos la 
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curvatura espacial cs fucrtcmcntc positiva y no plana. Sin embargo, en nuestro caso tales 
inconsistencias pueden ser sanadas (al menos parcialmente) si consideramos perturbaciones 
generales al modelo Bianchi-I, en los cuales la curvatura escalar y la 4-aceleraci6n no son 
nulas a orden perturbativo. En estos casos es posible examinar los efectos de acoplamiento 
tipo magneto-curvatura asociados con fuerza de Lorentz y parametro de desaceleracion 
no nulos en la ecuacion de Raychaudhuri (ver [56]). El uso de tales modelos Bianchi-I 
perturbados para la descripcion de un gas con fuente de materia en general sera presentado 
en el proximo capitulo. 

A la par de introducir posibles perturbaciones al modelo Bianchi-I, otra posibilidad 
que mejoran'a el modelo esta sobre la idea de una mejor descripcion dinamica de la fuente. 
Nucvas consideraciones para mcjorar esto, podn'a ser suponer modclos Bianchi I, V, VII 
6 IX con 4-velocidad inclinada, la cual introduce mas grados de libcrtad dinamicos, in- 
cluycndo hasta la posibilidad dc rotacion (ver [25] y referencias incluidas). Estos modelos 
podn'an permitirnos un mejor estudio(menos restrictivo) de los efectos dinamicos. Estos 
efectos son reportados en [56], en cuyo trabajo se muestra que la tension magnetica y 
el colapso gravitatorio pueden presentar un acoplamiento no nulo con la 4-aceleraci6n y 
vorticidad distinta de cero (incluso no pcrturbativa) en presencia de campo magnctico. 

Finalmente, un gas magnctizado contcniendo solamente neutrones puede ser tcoricamente 
intcrcsantc, pcro cs muy idcalizado como fuente para un objcto compacto. Asi, sc puede 
considerar como una extension a este trabajo, adicionar ademas de los grados de libertad 
dinamicos (mencionados arriba), mezclas de neutrones, electrones y protones, acompanados 
con sus respectivas ecuaciones de balance y potenciales qui'micos adecuados, en comparacion 
con otras ecuaciones de estado [46, 50, 51]. Todas estas extensiones estan bajo consideracion 
para trabajos futuros. 
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4. CONDICIONES PARA FRENAR EL COLAPSO 



Es posible mostrar, a partir de la ecuacion de Raychaudhuri (ecuacion (1.47) del Capitulo 
1), como el acoplamiento relativista entre magnetismo y geometn'a genera una tension 
magnetica que puede "enrollar" las Imeas de fuerza del campo magnctico. Este efecto 
tambien se identifica con una "elasticidad" de estas Imeas de fuerza del campo, las cuales 
tienden a mantenerse rectas. Concretamente, esta tension magnetica genera fuerzas de 
magneto-curvatura cuya manifestacion cinematica es la posibilidad de prevenir la conver- 
gencia de Imeas de universo (no necesariamente geodesicas) que inicialmente convergen. 
Evidentemente, la consccuencia principal de este efecto es la posibilidad de frenar el co- 
lapso gravitacional dc las fuentes asociadas a estas Imeas de universo, lo cual (en caso de 
suceder) ocurrirfa sin violar las condicioncs dc cncrgfa cstandarcs.'^^ 

Un resultado interesante que trata sobre el colapso gravitacional de un fluido magneti- 
zado sin simetria esferica fue obtenido por K. Thorne en el llamado "Universo magnetizado 
de Melvin" [59, 60, 61, 62, 23]. En dicho modelo el colapso del sistema se frena antes de 
alcanzar la singular idad, lo cual sucede debido al cambio de la densidad de energia de 
la materia magnetica. Mas recientemente, C.G. Tsagas ha encontrado [56] que existe un 
acoplamiento intrfnscco entre el magnetismo y la gcometrfa que surge dc las propicdades 
de la tension de las Imeas de fuerza magneticas (la "elasticidad"), y cuya consecuencia es 
el incremento a la resistencia al colapso. Esta tension magneto-geometrica es la general- 
izacion relativista del mismo efecto en la fisica newtoniana, apareciendo algunas veces en 
la litcratura como "frenado magnetico". 

En el presente capftulo examinamos, haciendo uso dc la ecuacion de Raychaudhuri, el 
efecto del campo magnetico sobre el colapso gravitatorio. Para llevar a cabo este objetivo 
tomamos en cuenta el trabajo de Tsagas en [56] y extendemos, al contexto astroffsico, el 
formalismo desarrollado por C.G.Tsagas y R.Martens [56, 57, 62, 63, 64]. Retomando las 
limitaciones planteadas en los dos capitulos anteriores, pasaremos a calcular las cantidades 
dinamicas covariantes que aparecen en la ecuacion de Raychaudhuri (ver Capftulo 1) para 
un modelo Bianchi-I perturbado. Trataremos la cuestion del colapso de la fuente magneti- 
zada en forma cualitativa, sin resolver las ecuaciones de Einstein-Maxwell. El objetivo de 
dicho estudio cualitativo es obtener informacion sobre el rol de las tensiones magneticas y 
entender las implicaciones potenciales de estas en modelar la evolucion de fuentes magne- 
tizadas en situaciones menos idealizadas. 

4.1 Espacio-tiempo Bianchi-I perturbado 

Usando 1-formas diferenciales^^ , el elemento de Imea de un espacio-tiempo Bianchi I per- 
turbado puede ser escrito como, 

3°Ademas de la tension magnetica, el tensor de deformacion distorsiona la distribucion de campos gradi- 
entes y, por lo tanto, tambien podn'a contribuir a resistir o favorecer al colapso. 

3^ Las 1-formas difcrenciales constituyen la base de cualquier tensor n- veces covariante o sea T = 
Taia2...anW°'^w'^^ ...w°'" , para mayor informacion ver [26] 
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donde, 



tiP= ^/ir^dx"^, (4.2) 

nf' =Qi^/T+^idx\ (4.3) 

nf' = Q2^/T+^dx'^ , (4.4) 

vf' = Q3^/l + ^3dx^, (4.5) 



y las = q^{t, x), = 0, 1, 2, 3 son las funciones que representan las perturbaciones a la 
metrica de Kasner, por lo que satisfacen \q^\ <C 1, mientras que las Qj = Qj{t),j = 1, 2, 3.^^ 
son las funciones no perturbadas, de modo que recuperamos la metrica de Kasner si q^ = 0. 
Elegimos la 4-velocidad ortogonal a las hipersuperficies espaciales, entonces, 

u^ = {l + qo)-'^^6^, u"uo^ = -l, (4.6) 

por lo que la 4-aceleraci6n iia = u^Ua;p toma la forma, 

1 



[0,5., go, 5., 90, 5.3 go]. (4.7) 



2(1 + qo) 

4.2 Calculo de cantidades dinamicas y ecuacion de Raychaudhuri 

Tomando en cucnta la metrica del modclo Bianclii I perturbado, las variables dinamicas 
asociadas a la ecuacion dc RaychaTidhuri (vcr capftulo 1) pucdcn ser directamente calcu- 
ladas. El escalar de expansion y el tensor dc deformaciones seran, 




® ^/TTq^ ^\Qi ' 2(1 + qi) 



1=1 

3 



^ [ln(g,yrT^)] , (4.8) 



< = diag[Si, S2, S3, 0], 
1 



(4.9) 



2Qi Qj Qk , q, 
+ 



qk 



' 3vTT^ \^ Qi Qj Qk 1 + qi 2(1 + qj) 2(1 + qk) j ' 
i ^j^k{i,j,k = 1,2,3). 

El parametro de deceleracion A = -i y el escalar de deformacion toman la forma, 



A 



1 ^ 1 r 

= 2(TT^ g Q|(i + g,) p'^" 



2(1 + go) 



+ 



^kHk ,sr^ 2 

l+Qk ^ 1 



dkQl 



+ Qk 1 + Ql 



(4.10) 



^^En este capftulo, a diferencia de los otros capftulos, usamos indices griegos para el espacio-tiempo e 
indices latinos para la parte espacial. Ademas, usaremos indistintamente para la derivada covariante el 
punto y coma o el si'mbolo "V" . 
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2^ 3 -3(l + go)^\l^Q. + 2(l + g,)^ 



-lyfg^+ ) f^ + ^^U. (4.11) 

2^^^^^ 2(l + gfe)yll^Qi 2(l + gOylj ^ ^ 

Por otra parte, nos sera litil definir el siguiente parametro adimensional de impacto de la 
deformacion, 

/ 3 / / ■ . \ 2 

^ I \ Qk , Qk 1 



^ e 1^10. + 2(1 



Qk 2(1 + Qk) 



k=l 



_ _ \" ^ J 

2f-igfc 2(l + gfc) M ' 2(l + g,); ( ( ysfOi+^i] 

1=1 \ '/ \ 'J )) Z^t=l yQi^ 2(l+qi) J 

Dada la existencia de una fuente magnetizada (con campo electrico nulo), el tensor de 
energi'a-momento viene dado por,^^ 

T^iu = (eo + ^)^i/iW,. + {po + -^)hnu + n^i/, (4.12) 

el cual corresponde al tipo de fuentes conocido como el fluido ideal "conductor perfecto" 
o de "conductividad infinita" [65] . Como PqJ jJ, son las presiones y la densidad de energia 
del gas isotropo, y H"^ la densidad de energi'a del campo magnetico [Jm~^], entonces los 
factores que miden las anisotropi'as son, 

w = — =po{t,x)/eo, h = — = B\t,il)/eo. (4.13) 
fj, fj, 

Sin embargo, la existencia de las perturbaciones implica un entorno anisotropico e 
inhomogeneo que generaliza la anisotropia del modelo Bianchi I. Por lo tanto, la ccuacion 
de conservacion asociada al tensor energia-momento que hemos definido en (4.12) viene a 
ser (ver [63]), 

(eo + Po + B^)u^^ + h^^ipo + ^B%, - h'' ,B'' .^B^^ - B'^B'^ = 0, (4.14) 

la cual tambien se puede escribir como, 

2 2 

(eo +Po + -^B = -D/^po -a^- Un^u", (4.15) 

donde D^po = h'j^VuPo denota el gradiente de presion medido por observadores en reposo 
comoviles con el fluido. Es importante enfatizar que todos los terminos en (4.15) se deben 

exclusivamcntc a la inhomogeneidad introducida por las perturbaciones. 

En cl modelo Bianchi I "puro" sin perturbaciones, al ser homogeneo y al no haber 
flujos termicos disipativos o efectos de rotacion que puedan generar perturbaciones fuera 



^^Note que aquf, a diferencia de los tratamientos anteriores, si hemos incluido la aproximacion de arbol, 
o sea la contribucion clasica del campo magnetico B^. 
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dc las condicioncs dc cquilibrio, es evidente que la anisotropia de un fluido ideal debe ser 
causada solo por el campo magnetico. Sin embargo, la existencia de las perturbaciones 
implica que las desviaciones de la homogeneidad y la isotropi'a no solo se deben al campo 
magnetico, por lo que podemos estudiar el acoplamiento de la geometn'a con el campo 
magnetico en un contexto mucho mas general que incorpora tambien a la inhomogeneidad 
(ademas de la anisotropia). 

Si llevamos a cabo un analisis gravito-magnetico, vemos que es el vector que deter- 
mina la aceleracion asociada a la fuerza de Lorentz, la cual siempre es normal a la direccion 
de las Imeas del campo a^Bf^ = (ademas se cumplc a^u" = 0), y surge cuando el modelo 
magnetico distorsiona las condiciones de equilibrio local. Este vector aceleracion tiene la 
forma. 



1 



-D^B' + B^D.B^. 

El primer termino del lado derecho de (4.16) viene dado por la "presion magnetica" y el 
segundo viene de la "tension del campo". Si este ultimo no es balanceado aparecera una 
fuerza neta sobre las particulas del fluido. Para el tensor energia-momento (4.12), el cual 
es anisotropico y corresponde a un fluido conductor infinito, todos los terminos en (4.15) 
son no nulos. La ecuacion cinematica de propagacion, es la ecuacion de Raychaudhuri vista 
en el capi'tulo 1 con A = 0, la cual reescribimos como, 



G 



e + ^ + 2a' + !^{l + 3w + h)-A = 0, 



(4.17) 



G 



1 



VI + ?o 



fe=i 



G^ _ 1 
3 " 3(l + (/o) 



qk 



+ 



Qk Ql 1 



ql 



qo 



qk 



Qk 

2(1 + qk) ' Qk Ql 2(1 + Qk)^ 2(1 + qo) \ 2(1 + q^) ' Qk 



+ 



E 

k=l 



Ql 



+ 



ql 



qkqj 



Ql 4(l + gfe)2 j^\QkQj 4(l + %)(l + g,)^ 



(4.18) 
(4.19) 



Ql . ql I Qkqk 



3(l + '?o)^lQ^ 4(l+gfe)2 Qk{l + qk) 



_ 1 / QkQi qkm qkQi 

2^VQfe<3i 4(l + gfe)(l + gz) {l + qk)Ql^ 
+ |(1 + Sw; + /i) = ^(eo + 3po + B"^) = 



j qk "^Qk 1 ql 1 ^0 qk 
I (^ -\-n,\ n, (^ -\- nA2 o 



2(l + go)frti (^ + ^^) 2(l + gfe)2 2 (1 + go) (1 + gfe) 

Qfc / go 2gfc \ 1 ^ 

QkW + qo) (i + gfc)yj 

Note que tanto el escalar de expansion como el tensor de deformacioncs son cantidades 
positivas, y si se cumple la condicion de energia fuerte, entonces el escalar de Raychaudhuri 
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{TlabU°"u^) debe ser mayor que cero, o sea, se cumple que IJ,{1 + 3w + h) > 0. Entonces, el 
linico termino que se opone al colapso gravito-magnetico es el coeficiente de deceleracion 
A (si es que resulta ser positive). 

De las ecuaciones de Einstein obtenemos para este caso, que las componentes = 
r| = Tg toman la forma, 

El termino de la ecuacion de Einstein sera, 



2' '^A 2(l + (/fc)(l + g;)Q2 s{l + qo) {I + qk){l + m) 



{l + qo)2Qk{l + qi) S {I + qk){l + qifQl ^' ' 



1 dkqkdkQi 1 QfcQz 

4(l + gjt)2(l + g,)Q2 2{1 + qo)QkQi 



l^k^m r-i / o o O 



rn= 



AQlil + qk){l + qm) \ (1 + 9/) (1 + 9fe) (l + ^m). 
y la curvatura escalar 3-dimensional toma la forma no nula. 



f- (1 + gfe)(l + g,)Q2 ] ^■^'^ + 2(1 + q^) ^ 2(1 + g,) 2^ 2(1 + g^) 

(4'.22) 

Note que la curvatura 3-dimensional (al igual que A) solo depende de las perturbaciones, 
y sus derivadas espaciales. Si pasamos a calcular el escalar de Ricci 4-dimensional, ob- 
tendremos expresiones similares, pero tambien dependeran de las funciones Qk y de sus 
derivadas, tanto espaciales como temporales. Para el caso no perturbado, se obtiene la 
simple expresion, 

4.3 Conclusiones del capi'tulo 

Considerando perturbaciones a la metrica de Kasner (modelo Bianchi I), hemos calculado la 
serie completa de parametros cinematicos que contiene la ecuacion dc Raychaudhuri, para 
una 4-velocidad comovil. Este resultado es importante desde el punto de vista teorico, 
pues nos permite examinar las condiciones de colapso para fuentes magnetizadas en el 
contexto de un espacio-tiempo inhomogeneo (aunque la inhomogeneidad es perturbativa) . 
Como se puede apreciar, al introducir las perturbaciones se complica cn gran medida la 
forma expli'cita de los resultados analiticos, pero estos nos aportan informacion cualitativa 
importante. 



-56- 



Se puede notar inmediatamente que las perturbaciones aparecen en todas las ex- 
presiones, pero una dependencia directa de las perturbaciones solo se manifiesta en el 
parametro de deceleracion A y en la curvatura 3-dimensional (^^TZ. O sea, estas son las 
variables que nos traen informacion de los efectos perturbadores del campo magnetico sobre 
las fuentes de materia-energia del sistema. 

El solo hecho de tener una curvatura S-dimensional no nula (aunque sea perturbativa) 
ya nos dice que este modelo perturbado (como supomamos desde las conclusiones del 
capi'tulo anterior) es mas realista que un Bianchi-I puro, el cual como ya sabcmos es piano 
en sus hipersuperficies tridimensionales. Si exigimos que la cantidad (4.22) sea positiva 
entonces logramos una mejor aproximacion a una descripcion "realista" de un volumen 
elemental dentro de un objeto compacto. Por otra parte, si la cantidad (4.10) resulta ser 
positiva, entonces los efectos de las perturbaciones tienden a frenar el colapso, y en caso 
contrario a ayudar a colapsar al sistema. 

Si proponemos formas analiticas expKcitas de las perturbaciones q^, podrfamos usar 
las formas generales de las soluciones analiticas que hemos dado para explorar los efectos 
cinematicos y dinamicos mencionados. Esto constituye un tema interesante a tomar en 
cuenta para investigaciones futuras, pero antes debemos encontrar (o al menos suponer) 
cual es la forma explicita de las perturbaciones que mejor corresponda a una condicion 
fisica menos idealizada del colapso gravitational de un sistema magnetizado. 

Un punto en el cual no podemos ofrecer prediccion alguna es en el tipo de singularidad 
que pudiera resultar del colapso gravitacional: singularidades tipo "cigarro" o "punto". 
El escalar de Ricci 4-dimensional tambien depende de las perturbaciones, de sus primeras 
derivadas espaciales, y de sus primeras y segundas derivadas temporales, y tambien depende 
de derivadas primeras y segundas temporales de los coeficientes metricos. Por tanto, para 
analizar su comportamiento en condiciones singulares (cuando este y los demas escalares 
de curvatura divergen) tendriamos que conocer la forma expli'cita de las perturbaciones 
(o al menos proponer formas especificas) . El analisis que hemos realizado cumple con el 
objetivo de este capitulo: encontrar resultados analiticos y cualitativos, lo mas generales 
posibles, que contribuyan al entendimiento de la interaccion entre la geometri'a y el campo 
magnetico, asi como faciliten la elaboracion de estudios numericos para casos particulares 
que podn'an ser llevados a cabo en el futuro. 

La limitacion principal del modelo Bianchi I perturbado es que las funciones qk deben 
ser muctio menores que la unidad, de modo que desde el inicio ya estamos imponiendo 
una fuerte limitacion a la geometri'a y a la dinamica de las fuentes (como sucede en todo 
tipo de analisis perturbativo) . Si las perturbaciones infiuyen sobre la intensidad del campo 

magnetico, entonces dicho campo nunca llegarfa a veneer a la intensidad del campo gravita- 
torio. Por esta razon, para lograr una mejor descripcion y entendimiento de la dinamica del 
colapso gravito-magnetico nos vemos en la necesidad de considerar para trabajos futuros, 
modelos Bianchi con mas grados de libertad que el modelo I. 
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Part II 

PEQUENOS AGUJEROS NEGROS 
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5. ENTROPIA DE AGUJEROS NEGROS 



5.1 Entropia de agujeros negros: pelos clasicos 

Existe una simpatica expresion sobre los ANs, que soli'a repetir John Wheeler, "los ANs 
no tienen pelos" . Con esto queria decir que, excepto por unas pocas caracten'sticas que los 
distinguen, todos los ANs resultan parecidos no exhiben ningun "peinado" personal, algo 
caractcn'stico que nos permita diferenciar un AN de otro. Lo linico que los caracteriza son 
su masa, su carga electrica, y su momento angular. Nadie podn'a distinguir dos ANs con 
estos mismos valores caracteristicos, y, precisamente, esta similitud de los rasgos definitorios 
ha hecho creer, a lo largo de los anos, la "extraha suposicion" de que los ANs podrian ser 
gigantescas parti'culas elementales. 

La Teoria de Cuerdas nos ha permitido avanzar en esta dircccion pues variedades 
como la dc Calabi-Yau^'^ presentan interacciones intcrnas cntrc sus subvariedades, que 
apantalladas por Branas^^ envolventes, permiten obtener soluciones de AN sin masa, esto 
es lo mas similar a un foton o sea pura radiacion electromagnetica en nuestro mundo. 

Por otro lado y menos especulativo, los ANs deberian estar rodeados por un anillo 
dc luz que revela sus propiedades fundamentales. Muy pocos astronomos ponen en duda 
la existencia de los ANs, aunquc la evidencia de su existencia es indirccta, a partir del 
comportamicnto dc otros objetos ccrcanos a los ANs, talcs como estrellas brillantes. 
Estos objetos literalmente podrfan ser engullidos por el AN y en caso de dar sehales medi- 
bles de tal comportamiento, entonces esto constituin'a una prueba de su existencia. Obvia- 
mente, para la comunidad de Astronomos y para los Fi'sicos, la deteccion directa contimia 
siendo un desafio. 

Por ejemplo, en los dos liltimos ahos, sc ha utilizado la interferometrfa para crear 
imagcncs del AN central dc nucstra galaxia. Estas imagcncs junto con las que sc cspcran 
obtener en los proximos afios deberian ser capaces de resolver detalles sobre la escala del 
horizonte de sucesos del AN. 

Retornemos entonces al conocido "teorema del no pelo" que naturalmente impone unos 
Ifmites estrictos a la influencia que el AN pueda ejerccr en su entorno. Imagenes de un AN 
o mas bien de la materia que lo circunda, podrfan refutar cl teorema, si se prueba que dicha 
materia se distorsiona de alguna manera. Tim Johannsen y Psaltis Dimitrios [66], de la 
Universidad de Arizona, en Tucson, dicen que hay una manera facil de probar el teorema. 
En la actualidad, calculan que los ANs deben estar rodeados por un anillo de luz. 

Esta luz proviene de los fotones que han quedado atrapados en una orbita circular 
alrededor del AN, en las afueras del horizonte dc sucesos, que luego se dispersa por el gas y 
cl polvo que cae en el AN. Este anillo debe ser mucho mas brillante que el gas y el polvo que 
lo rodea, debe tener un diametro que es unas diez veces el tamaiio del AN, lo que significa 
que deben'a ser visible en las imagenes que pronto estaran disponibles y que proporcionaran 
una medida directa de la masa del AN. La forma de este anillo, depende de las propiedades 
del AN, y no de la estructura de la materia que cae en el AN. Eso significa que la forma del 

^*Es una variedad con 6-dimensiones enrolladas (compactificadas), y presenta dos esferas empotradas una 
de 2-dimensiones y otra de 3-dimensiones que dinamicamente pueden interactuar. 
^^Es un objeto extendido, sobre el cual viven las cuerdas. 
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anillo cs la medida de las propicdades del AN y cualquier asimetria en el anillo sera una 
violacion directa del teorema del no pelo. Los astronomos no tienen que mirar muy lejos 
para encontrar tal anillo. "El AN del centre de la Via Lactea" , es el candidato ideal para 
probar el teorema del no pelo debido a su brillo, a su gran tamano, momento angular, y 
relativa cercani'a. 

La solucion estacionaria de Schwarzschild (1916) la discutimos en el capitulo 1 y 
la ecuacion (1.4), muestra cl cstado final que debc tcncr todo colapso gravitacional con 
simetria central. Historicamente, despues surgieron otras dos familias de soluciones, la 
solucion con carga electrica "e" de Reissner (1916)-Nordstrom (1918) y posteriormente, 
la solucion de Kerr (1963) con momento angular "J = aM". Una generalizacion de la 
familia de Kerr, cargada fue obtenida por Newman en 1965. Esta ultima solucion de Kerr 
cargada esta formada por tres-parametros, y la mctrica del espacio-tiempo junto con el 
vector potencial electromagnetico son dados matematicamcnte por la ecuacion siguiente, 

(fs = -(^ J dt^ ^— ^-dtdcl) + 



+ 



sin^edcj)^ + ^dr'^ + T.de'^, (5.1) 



Aa = --[{dt)a -asin^9{dct>)a], (5.2) 



donde. 



S = + cos'^e , (5.3) 
A = + + - 2Mr, (5.4) 

y e, a, y M son los tres parametros de la familia. Cuando e = 0, tenemos a Aa = 0, la 
metrica del espacio-tiempo se reduce a la familia de soluciones de Kerr. Cuando a = se 
recupera la solucion de Reissner-Nordstrom y cuando e = a = 0, la ecuacion (5.1) se reduce 
a la solucion de Schwarzschild. Asf, todas las soluciones de ANs, estacionarias conocidas 
estan contenidas dentro de la familia 3-parametrica, segiin la TGR. 

En todas las soluciones que tienen una interpretacion fi'sica directa, aparecen los 3- 
parametros e, a y M. Para cualquier 2-esfera, S, en la region asintotica, tenemos que, 

l [ eabcdF"' dx'' A dx^ = 47re, (5.5) 
f ^ahcdy"^'^ dx^ Adx'' = M, (5.6) 

OTT Js 

f eabcd^'^P^'dx'' Adx' = Ma, (5.7) 

lOTT Js 

donde, = (d/dt)", y ^Z'" = {d/d^f son los vectores de Killing (ver [27])^^. 
5.2 Mecanismo de Sen 

Es conocido que la Tcorfa dc Cucrdas cn cl Ifmite de bajas energfas surgen acoplamientos 
entre la gravcdad y otros campos. Estas tcorfas tambien contienen soluciones tipo AN. Asf 



^Para mas informacion sobre la notacion usada en (5.5)-(5.7), ver pag 151, [26] o Apendice B de [27]. 
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la teoria de cucrdas brinda un cuadro de trabajo para estudiar las propiedades de los ANs, 
sean clasicos o cuanticos [67]. 

Clasicamente los ANs son soluciones de las ecuaciones de Einstein con propiedades 
especiales. Todos los ANs tienen una superficie hipotetica que los rodea, llamada hori- 
zonte de eventos de donde ningun objeto puede escapar despues de atravesar su frontera. 
En la teoria cuantica los ANs se comportan como cuerpos negros con temperatura finita 
conocida como "temperatura de Hawking". Asi mismo emiten radiacion Hawking en cor- 
respondencia con las leyes de radiacion del cuerpo negro, c interactuan con la materia 
como un sistema termodinamico. Si nos basamos en leycs aproximadas donde se ignoran 
terminos de derivacion superior a segundo orden, la entropia de los ANs vendra dada por 
la expresion de Bekenstein-Hawking, 

Sbh = A/{AGn), (5.8) 

donde A es el area del horizonte de eventos y Gjv es la constante de Newton. Una impor- 
tante pregunta es si: "^Podremos entender esta entropia desde el punto de vista estadistico, 
o sea como el logaritmo del mimero de microestados cuanticos asociados al AN ?" . 

Aunque no se tiene una respuesta definitiva a esta pregunta, para una clase especial de 
AN en teorfa de cuerdas conocido como "AN extremal" , esta pregunta ha sido contestada 
y en forma afirmativa. Estos ANs tienen temperatura cero por lo que no emiten radiacion 
Hawking y son usualmente estables. Algunas veces (no siempre) ellos son invariantes ante 
cierto mimero de transformaciones de supersimetria y en ese caso se llaman "ANs BPS" . 
Debido a la estabilidad y a las propiedades supersimetricas, se tiene cierto control sobre 
la dinamica de la configuracion microscopica (esencialmente envuelven D-branas, cuerdas 
fundamentales y otros objetos solitonicos) de estos ANs. Esto permite calcular la degen- 
eracion de tales estados con un acoplamiento debil donde los efectos dc "backreaction" 
del sistema pueden ser ignorados. La supersimetria nos permite extender el resultado al 
caso de acoplamiento fuerte donde el "backreaction" gravitational pasa a ser importante y 
el sistema puede ser descrito como un AN. En teoria de cuerdas nos encontramos con una 
gran variedad de ANs extremes tipo BPS, que en el li'mite de AN grande [68], la entropia 
toma la forma, 

Sbh{Q) = SestiQ), (5.9) 

donde Sbh{Q) es la entropia de Bekenstein-Hawking de un AN extremal debido a las 
cargas "Q", y SestiQ) es definida como, 

SestiQ) = In d{Q), (5-10) 

donde d{Q) es la degeneracion de los estados BPS en la teorfa debido al conjunto de cargas 
"Q" . Esto nos brinda un buen entendimiento de la entropia de Bekenstein-Hawking desde 
el punto de vista microcopico. 

^^En Hsica teorica la "backreaction" es frecuentemente utilizada para calcular el comportamiento de una 
parti'cula de prueba u objeto, en un campo externo. Por ejemplo, si la particula de prueba tiene masa o 
carga infinitesimal entonces se ignora la "backreaction" , sin embargo dicha particula posee carga y masa no 
nula lo que altera a su vez el asumir campo extorno fijo. Por esto, en ciertos casos se debe tomar en cuenta 
su efecto, entonces los efectos de "backreaction" no se podran despreciar. 
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Las comparaciones entre Sbh Y Sest fueron hechas primeramente en el li'mite de cargas 
grandes. En este li'mite el horizonte es grande de modo que la curvatura y otros campos de 
fuerza en el son pequenos y podemos calcular la entropia usando (5.9) sin preocuparnos por 
correcciones con terminos de orden superior de derivacion^^ a la accion efectiva de cuerdas. 
Ademas el calculo de SestiQ) se simplifica porque la dinamica del sistema microscopico es 
frecuentemente descrita por una Teoria de Campo Conforme(TCC), 1 + 1 dimensional con 
coordenada espacial compactificada en un circulo. Un AN extremo fuertemente cargado, 
se corresponde con estados en esta TCC con los autovalores Lq (o Lq) y cero Lq (o Lq). La 
degeneracion de tales estados puede ser calculada usando la formula de Cardy en terminos 
de la carga central izquierda y derecha {cl,cr), sin conocer los detalles de la teon'a TCC, 
o sea, 



para Lq = 0, (5-11) 



para Lq = 0. (5.12) 

Entonces constituye una sorpresa que dos calculos completamente diferentes (uno para 
Sbh{Q) y otro para SestiQ)), ^cn resultados iguales. Dc cstc modo cuando nos mova- 
mos lejos del li'mite de cargas grandes, la curvatura y los restantes campos de fuerza en 
el horizonte no seran despreciables. Asi tendremos que tomar en cuenta el efecto de los 
terminos de derivacion superior en la accion efectiva de la teon'a. Terminos de derivacion 
superior son por ejemplo potencias cuadraticas y de mayor orden del tensor de Riemann. 
Nosotros esperamos que para los ANs grandes (no infinites) los efectos de estos terminos 
de derivacion superior sobrc el horizonte sean pequenos pero no nulos, estos efectos intro- 
duciran pequehas modificaciones a la geometn'a del horizonte y a la entropi'a de los ANs. 
Por otro lado, la formula de Cardy para la entropi'a estadistica, si se calcula fuera del li'mite 
de AN grande, tambien recibira modificaciones inversas a potencias de la carga del AN, 
acentuandose sorpresivamente cada vez mas la similitud entre Sbh y Sest, incluso para 
cuando tenemos en cuenta terminos de derivacion superior. 

La entropia Sbh de Bekenstein-Hawking fue calculada originalmente en TGR. Sin em- 
bargo, pretendemos aqui' incorporarle al resultado terminos de derivacion superior. Varias 
dualidades^^ en teorias de cuerdas pueden mapear contribuciones clasicas a la teoria efec- 
tiva cuantica, por tanto, no tiene sentido restringir nuestro analisis a la teon'a clasica. Una 
eleccion es elegir la accion efectiva que genera los diagramas irreducibles de una parti'cula 
(IPX), esta respeta todas las simetrias de dualidad. Sin embargo, la teon'a de cuerdas tiene 
particulas no masivas, y si vamos a ordenes mayores en terminos de derivacion superior, 
la accion IPX contendra contribuciones no-locales. Como hemos visto aqui, para cualquier 
teon'a de gravedad con terminos de derivacion superior y con una densidad Lagrangiana 
local existe un buen algoritmo para calcular la entropia de ANs, hasta el presente no exis- 
ten tecnicas para tratar con teon'as con accion no-local. Esto causa un potencial problema 

^*Se les conoce por "terminos de orden superior de derivacion", porque introducen derivadas de la metrica 

superiores al segundo orden, cn las ccuacioncs dc campo. 

^^Las dualidades son transformaciones matematicas que relacionan las 5-teon'as de cuerdas existentes. 



SestiQ) - 27r 
~ 27r 



clLo 



crLo 
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para definir entropia dc un AN en teorias de cuerdas, mas alia de los primeros ordenes. 
Se podn'a esquivar este problema usando una accion efectiva local de Wilson, pero esta no 
respecta todas las simetrias de dualidad de la teoria. 

Tambien se presentan ambiguedades en el lado estadistico, donde las correcciones a la 
entropi'a dependen del inverse de la carga del AN, y se sabe que tales correcciones depen- 
den del ensemble que elegimos para definir entropi'a. Por ejemplo podemos usar ensembles 
micro-canonicos o gran-canonicos invariantes de dualidad, o podemos usar ensembles mix- 
tos no invariantes de dualidad [69]. 

Esperamos que los ejemplos explicitos para el calculo de las correcciones a la entropi'a 
de ANs, puedan resolver las ambiguedades anter lores y hacer una formulacion mas pre- 
cisa sobre la relacion entre las dos entropi'as. Nosotros nos dedicaremos al calculo de la 
entropia Sbh de ANs teniendo en cuenta terminos de derivacion superior, basandonos en 
el formalismo cntropico [70] que es una adaptacion de un formalismo mas general [71] al 
caso especial de ANs extremales. La scgunda li'nca de trabajo o sea calcular la entropi'a 
estadistica se puede encontrar en [72, 73, 74], y [67] donde incluso se comparan las dos vias 
para el caso de un AN en una teon'a de cuerda supersimetrica en 4-dimensiones. 



5.3 Definicion de AN extremal 

Analizemos un AN de Reissner- Nordstrom (RN) , que describe un AN esfericamente simetrico 
con carga electrica en la teoria usual de Einstein-Maxwell, 

S = j d^x./^d^jr, jr = Y^^R - (5.13) 

La solucion de RN en esta teon'a viene dada por, 

ds^ = -(1 - a/p)(l - b/p)dT^ + -^7- — - + p\de^ + sin^9d4?), (5.14) 



{l-a/p){l-b/p) 



^pr = Fep = {-sine, (5.15) 



donde p,6,(p y t son las coordenadas del espacio-tiempo, a y b dos constantes dadas por, 

a + 6 = 2GjvM, ab = ^{q"^ + p"^), (5.16) 

47r 

donde q, py M que denotan la carga electrica, magnetica y la masa del AN. Cuando a> b 
entonces el horizonte inter no estaria en r = 6 y el externo en r = a. El Kmite extremo sera 
cuando r = a = b, entonces. 



Ahora si definimos, 

t = \T/a^, r = X-^{p-a), (5.18) 
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donde A es una constante arbitraria, y reescribimos la solucion extrema en estas nuevas 
coordenadas, obtenemos que, 

ds'^ = ''^t ,o dt^+ ^ ^^"^^"^ dr^ + (a + \rf{de'^ + sm^e d(P'^), (5.19) 
(a + ArJ^ 

y en el limite muy "cercano al horizonte" A — > 0, la solucion tomara la forma, 

ds^ = a^(^-r^dt^ + ^)+a^{de^ + sin^edct^^), (5.21) 
Frt = Fo^ = £-sin9. (5.22) 

La entropia que se obtiene dividiendo el area del horizonte por 4Gjv es, 

SBH = \{q^+P% (5.23) 

y tiene las propiedades siguientes: 

• En el limite A — t- manteniendo a r fija, la coordenada original p se aproxima a a. 
Asi (5.21) describe el campo de configuracion de un AN cercano al horizonte. 

• El espacio de (5.21) se desdobla en dos espacios. Uno etiquetado por {0,(p) descrito 
por una esfera S"^ de 2-dimensiones, y otro etiquetado por (r, t) descrito por un 
espacio-tiempo Anti dc Sitter (AdS) de 2-dimensiones. Este ultimo rcsulta de la 
solucion de las ecuaciones de Einstein en 2-dimensiones con constante cosmologica 
negativa. 

• El reflejo de la simetria esferica del AN original hace que en (5.21) se presente una 
isometria SO{3) que actua sobre la esfera S^. En (5.21) tambien se presenta una 
isometn'a -SO (2, 1) actuando sobre el espacio AdS2, esta es generada por, 

Li = dt, Lo = tdt - rdr, L_i = l{\+ t^)dt - trdr. (5.24) 

Asi, tanto la metrica como los campos de calibracion son invariantes bajo una trans- 
formacion del tipo 50(2, 1) x 50(3). 

Entonces si postulamos que la geometn'a del AN no cambiara si afiadimos tcrminos 
de derivacion de orden superior, y en el case dc D-dimensiones se puede generalizar a que 
los ANs tengan isometn'as del tipo SO{2, 1) x SO{D — 1), entonces todo esto constituye 
nuestra definicion de "AN extremo"^°. 



■'Para 4 y 5-dimensiones, estos postulados han sido probados [77]. 
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5.4 Metodo de Sen 

Siguiendo [75, 76, 70], supongamos una teoria acoplada a campos de calibracion A^^^ y a 
campos escalares (ps- La densidad Lagrangiana sera expresada en funcion de la mctrica, de 
estos campos, y de las derivadas covariantes de estos campos. Ademas, suponicndo ANs 
extremes, esfericamente simetricos con invariancia tipo SO{2, 1) x SO{D — 1) cerca del 
horizonte, la configuracion de campos mas general y consistente con esta isometn'a es, 

-r^ dt^ + —J + V2 [de"^ + sin^e dcp'^), (5.25) 

ct^s = Us, F^f = a, F^l = ^sinO, (5.26) 

donde ui,f2,Us,ej y Pi son constantes. Para esta metrica las componentes del tensor de 
Riemann son, 

Rabcd = -Vi^iQacQbd " gad9bc), a, b,c,d = r, t, (5.27) 

Rmnpq = {gmpQuq " 9mq9np), m, n,p,q = 9, (p. (5.28) 

(i) 

Las derivadas covariantes de los campos escalares 0^, los campos dc calibracion 
y del tensor de Riemann Rabcd se anulan cerca del horizonte de geometn'a. Debido a con- 
sideraciones generales de simetn'a, todas las contribuciones a las ecuaciones de movimiento 
de cualquier termino en C que contenga derivadas covariantes de los campos en general 
se anulan para esta metrica, entonces podemos concentrarnos en aquellos terminos que no 
contengan derivadas covariantes de los campos. 

Dcfinamos la densidad de lagrange / (u, v, e, p) evaluada cerca de la geometria del 
horizonte (5.25) e integrada por las variables angulares, 

f{u,v,e,p) = J dOd^^-detgC. (5.29) 

Las ecuaciones de campo cerca del horizonte corresponden a extremar / respecto a las 
variables u y v, 

^ = 0, 1^ = 0. (5.30) 

OUs OVi 

Las componentes no trivialcs dc las ecuaciones del campo dc calibracion y de las 
identidades de Bianchi para la solucion completa del AN toman la forma, 

5r(^)=0, a,F» = 0, (5.31) 

donde S = J d'^x \/ —del g C cs la acccion. Estas ecuaciones son satisfechas cerca del hori- 
zonte de (5.25) y podemos extraer mas informacion de ellas, pues evaluando estas integrales, 



[ dOdcP-^ = ai, [ ded<l)F^l = hi, (5.32) 
cerca del horizonte (5.25), nos da que. 



ai = ^, bi=pi. (5.33) 
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Ellas son justamente la carga electrica a^, y la carga magnetica bi. Denotemos por, 

como la carga electrica del AN. Ahora el conjunto (5.30) y (5.34), forma un sistema de 
ecuaciones con igual niimero de incognitas u,v y e. En la mayon'a de los casos este sistema 
se puede resolver y quedar determinado en funcion de q y p. Esto es consistente con el 
mecanismo atractor en supersimetria el cual enuncia que cerca de la configuracion de un 
AN, las magnitudes dependen solamente de la carga electrica y magnetica del AN y no del 
valor de sus campos escalares. 
Definamos ahora, 

£{u,v,e,q,p) = 27r(eigi - f{u,v,e,p)). (5.35) 

Las ecuaciones (5.30) y (5.34), determinando u,vy e pueden ser escritas como, 

d£ d£ d£ ^ d£ ^ 

O sea, todos los parametros cerca del horizonte pueden ser calculados extremando 
la funcion £. Por otro lado, la formula general planteada en [71] y escrita para ANs 
esfericamente simetricos, tiene la forma, 

Sbh = -Svr / d9d(f)— — V-grrgtt, (5.37) 

JH Ottrtrt 

donde H denota el horizonte del AN. Despues de una serie de pasos a partir de (5.37) se 
Uega a que (ver [67], pag 14), 

SBH = lT^[ei^-f). (5.38) 

Comparando (5.38) con (5.35) se obtiene que, 

Sbh = £iu,v,e,q,p). (5.39) 

O sea Sbh{q,P)/'^'^ puede ser interpretado como la transformacion de Legendre de la 
funcion f{u,v,e,p) con respecto a la variable e^, despues de eliminar u y v mediante el 
sistema (5.30). 

Asi el formalismo de la funcion entropica, reduce el problema del calculo de la entropi'a 
de un AN extremal a un problema de solucionar un sistema de ecuaciones. Veamos un caso 
muy simple, el AN de RN, 

Frt = e, Fe^=psin9/iTr. (5.41) 
Usando (5.13), (5.29) tenemos que, 

nv,,v,,e,p) = 4-1^2 [^(-^ + ^) + IvT^e^ - ^v^H^f] , (5.42) 
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y que, 

£{vi,V2, e,q,p) = 27r(ge- /) = 27r 



qe- 77^(2^1 - 2v2) - 2'kv2V^ ^e^ + 2'kviV2 '^{-j-)'^ 
4Gjv 47r 



Aqui £ es extrema si, 



Vi=V2 = Gn 



4tt ' 



e = 



47r' 



Sustituyendo en la expresion para £, tenemos que, 

SBH=£=\iq^+P% 



(5.43) 
(5.44) 

(5.45) 



que coincide con al entropfa clasica de Bekenstein-Hawking para un AN de RN. 
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6. ENTROPIA DE AGUJEROS NEGROS EXTREMOS 



6.1 Nuevos resultados en la aplicacion del formalismo de Sen 

Motivados ya sea por lograr la total unificacion de las fuerzas fundamentales de la Nat- 
uraleza o por tratar de explicar fenomenos que hoy en di'a permanecen sin resolver en la 
fisica teorica, la TGR ya ha sido generalizada de muchas formas. Unas de estas vi'as es 
apoyarse en el contenido geometrico que posee toda teoria relativista y de esta forma lograr 
dicha generalizacion. Esto es especi'ficamente lo que hacemos en este capi'tulo aplicado al 
caso particular de ANs extremos. 

En la fisica es muy importante conocer los invariantes que posee la teon'a si se piensa 
proponerla como una teoria nueva, asi se hace por ejemplo en electrodinamica, EDC o 
CDC^^ donde el invariante fundamental es, 

-F'^'Fah, (6.1) 
a 

donde a es la masa del gap.^^ Si sobre ellos aplicamos el metodo extremal se obtienen las 
ecuaciones de campo dc la teon'a. Ya sea las ecuaciones de Maxwell para la EDC o las de 
Yang Mills en CDC. Similarmente en la TGR el invariante principal es el escalar de Ricci 
R y las ecuaciones de Einstein se obtienen extremando su accion. Aunque este invariante 
es indispensable para la formulacion de toda la teoria de la TGR, no constituye el linico 
invariante posible, existen muchos mas que se pueden construir. Por ejemplo, usando el 
tensor de Ricci Rab, J el tensor de Riemann Rabcd, se obtienen los invariantes, 

R^, RabR°'^i RabcdR"'^'^'^ ■ (6.2) 

Al incorporarlos a ellos se define el sector geometrico (o gravitatorio) de la teoria. 
Claramente esto mismo se puede hacer con la accion de la EDC y en el marco de las teorias 
de campo se define como trabajar en el sector de calibracion. Ademas cuando decimos 
que trabajamos con invariantes de or den superior (o de or den superior de derivacion) en el 
sector geometrico, significa que dicho invariante es de mayor orden que R, la interpretacion 
es similar en el sector de calibrado para invariantes de ordenes mayores^^ que F^. Hoy en 
dia la incorporacion de invariantes de orden superior es un tenia abierto. 

En este capitulo calculamos la entropia de ANs extremos y estacionarios teniendo en 
cuenta terminos de derivacion superior. La inclusion de estos terminos en la teon'a de alta 
gravedad puede hacerse usando por ejemplo los terminos de Gauss-Bonnet (GB) [88, 89, 90]. 
Estos terminos pueden aparecer en varios escenarios teoricos como; Teorias de Cuerdas, 
Teon'as de Branas [91, 92, 93] y Gravedad Cuantica Semiclasica. En particular, terminos 
de derivacion superior aparecen en la teon'a de cuerdas cuando se evaliia el h'mite efectivo 
de baja cnergi'a. 

*^EDC (Electrodinamica Cuantica) y CDC (Cromodinamica Cuantica). 

^^Determinar la masa del gap es un problema abierto, y esta dentro de los 23 problemas de Hilbert 
(Axiomatizar toda la fisica), tambien se encuentra entre los problemas del milenio. Esta masa se corresponde 

con la primera excitacion o estado basico fundamental con respecto al vacfo. 

*^Asf mismo como = RabR'^^ , el termino = FabF'^^ , se define similarmente. 
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Por otro lado, se conoce el problema de no-renormalizacion que sufren las teorias de 
gravedad. En especial la teon'a de GB si puede ser renormalizada pero a su vez aparecen 
otros problemas de igual importancia, emergen fantasmas (particulas de spin = 2 masivas) 
o gravitones masivos [86, 87]. 

En este capitulo presentaremos una serie de soluciones aproximadas para la entropia 
de ANs extremalcs ccrca del horizonte de geometn'a para cuatro y cinco dimensiones. Asi 
para alcanzar esta mcta, usaremos el formalismo do Sen. La teoria de partida sera detallada 
en el epi'grafe 6.2, pero en efecto ella esta compuesta por una teon'a de Einstein-Maxwell 
con constante cosmologica, ademas del conjunto completo de los invariantes de Riemann. 
Es importante senalar que en aquellos casos en que las ecuaciones de movimiento (en 
la cercania del horizonte), no se puedan resolver analfticamente, entonces efectuaremos 
Tina aproximacion por series de potencias basandonos en un parametro inverso a la carga 
electrica del AN. 

Por otro lado, si queremos incluir el conjunto completo de los invariantes de Riemann 
como terminos de derivacion superior de la teon'a, debemos considerar la version mas actual 
de este conjunto. O sea, el conjunto definido por Carminati y McLenagham (CM) [30], 
mas cl invariantc mg introducido por Zakhary y Mcintosh [31]. 

Los invariantes CM son escalares construidos a partir del tensor de Riemann, Rabcdi el 
tensor de Weyl Cabcd ij su dual) y cl tensor dc Ricci dc traza mila, definido como Sab = 
Rab—{^/d)Rgab (en d-dimensiones) . Comenzando por seis escalares reales R, ri, r2, ra, ma, m4 
y cinco escalares complejos wi,W2,'mi,m2, m^, ellos hacen un total de dieciseis invariantes 
escalares (Note que los si'mbolos r, it; y m se asocian respectivamente a invariantes de Rie- 
mann, Weyl y su mezcla). El conjunto de invariantes de CM contiene el mimero requerido 
de invariantes para los casos de Einstein-Maxwell y el fluido perfecto. Ademas, con la in- 
clusion del invariante niQ, se ha probado que el conjunto de CM se convierte en un conjunto 
completo. De esta forma este conjunto cubre los 90 casos posibles (6 tipo Petrov x 15 tipo 
Segre, ver [31]). En otras palabras, un conjunto completo de invariantes debe contener (al 
menos para el espacio-tiempo d = 4) los bien conocidos invariantes fi'sicos, asi' como los 
invariantes geometricos. 

En el Cuadro 2, (ver Apendice G) mostramos las definiciones y resultados obtenidos 
para el conjunto de invariantes no nulos de Riemann sobre una geometn'a AdS2 x 3"^^^ 
para d = 4 y d = 5. Es importante senalar que para el caso de d = 5 no es posible calcular 
los invariantes complejos de Riemann, ellos son desconocidos. Para esta dimension podria 
requerirse de mas invariantes y la pregunta seria, cuantos invariantes se necesitan para que 
estos conjuntos con d > 4 sean completos. Este todavi'a es un problema abierto. 

Igualmente en el Cuadro 3, (ver Apendice G), se muestran los invariantes organizados 
por su grado. Siguiendo las definiciones dadas en [31], de un conjunto dado de invariantes, 
el j— esimo invariante escrito de la forma 7j, es Uamado invariante de orden p. Si existe 
otro invariante tal que pueda ser escrito de la forma IjI^If, entonces se dice que este es un 
invariante de orden p—q—r y el resultado de la suma;j-|-g-|-r sera su grado. Por ejemplo, en 
el Cuadro 2, el invariantc 9ie(mi), dcnota la parte real del invariante mi, tienc orden 1 — 2 
y es de 1 + 2 = 3, tcxccx-grado. Ademas, el invariante sera independiente si el no puede 
ser escrito en terminos de otro invariante, ya sea de igual o menor grado. Tambien, dos 
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invariantes se dice que son equivalentes si ellos pueden ser escritos en terminos de otros, 
o como el producto de otros de menor grado. Las relaciones de equivalencias o (syzygies) 
son dadas en el Cuadro 3. Es facil de ver que en 4 y 5 dimensiones todos los invariantes 
pueden ser escritos en terminos de los invariantes independientes R, ri y r2. Es importante 
senalar que todos los invariantes no nulos fueron tornados en cuenta en nuestros calculos, 
los restantes nulos no son incluidos en los Cuadros 2-3. 

Naturalmente, una pregunta sale a flote, y es: ipoi que no estan presentcs los invari- 
antes de orden superior de la teorfa de calibracion?, dado que en efecto pensamos incluir la 
teoria de Einstein-Maxwell donde los campos de calibracion estan presentes. La respuesta 
es, que lo hacemos por simplicidad, o sea solo trabajaremos en el sector de gravedad pura, 
para los terminos dc ordcn superior dc la teorfa, quedando fuera invariantes de orden 
superior como: {FabF^'^ f, F" ^ cF'' dF'^ a, RF^,..., los covariantes y los in- 

variantes de forma. Es decir solamente el sector gravitational sera tratado para el conjunto 
de CM. 

Todos los calculos fueron hechos con el paquete para tensores GRTensor sobre el 
programa de computo algebraico Maple. 

6.2 Generalizacion de la Teorfa de Einstein usando los invariantes de Riemann 

Consideraremos una teoria de gravedad de ordcn superior en el sector geometrico, intro- 
duciendo el conjunto completo de los invariantes de Riemann como terminos de derivacion 
superior en la accion. Ademas incluimos constante cosmologica, y campos electromagneticos. 
Asf nuestra accion toma la forma, 



?2 



1 /■ , ^ . . F 



donde, 



Cf-^ = 0,2 + 62i?2 + a^R^ + h^R i?2 + a^R"^ + b^Rl + c^R^ R2 + a^R^ + hR^ R2 , 
Cf=^ = a2R^ + 62i?2 + a^R^ + bsRRi + C3R3 + a^R"^ + b^Rj + c^R^ R2 + e^RRs , 

siendo Gd la constante de Newton d-dimensional, R el escalar de Ricci, A la constante cos- 
mologica. Fab el tensor electromagnetico, y = FabF"-^ . R2 y R3 son los dos primeros in- 
variantes de Riemann reales definidos en el Cuadro 2, para la metrica (6.4). Los parametros 
ai,bi,Ci,ei... son las constantes de acoplamiento para cada termino de derivacion superior 
de i-esimo grado. Cf^^^ denota los terminos de derivacion superior que se tendran en cuenta. 
En ambos casos, en cuatro y cinco dimensiones, los invariantes forman un conjunto com- 
pleto, y se tiene en cuenta todos los invariantes de alto orden dentro del conjunto de 
invariantes de Carminati-McLenagham. 

El espacio tiempo mas general [67] que se puede escribir para una configuracion de AN 
estatico y extremal con topologia AdS2 x S'^~^ cerca del horizonte de geometria es, 

ds'^ = vii-r'^dt'^ + 7^) + V2dnd-2, (6.4) 
e'^'^'lH = ui, F^, = e, F^T* = psinO, Flf = 0, (6.5) 



-70- 



d-2 i-1 

dnl_2 = d9f + ^Yl sin^ ^i^^i^ < 0i < TT, < 0d_2 < 27r, (1 < i < d - 3), 

i=2 j=l 

donde e y p son funciones relacionadas con la carga electrica y magnetica, mientras que 
vi y V2 son funciones conectadas con la garganta del AN. En general las constantes Uj se 
relacionan con los valores de los campos escalares sobre cl horizontc, pcro en nuestro 
caso ellos no estan presentes por tanto Ui = 0. En los epi'grafes siguientes usaremos el 
formalismo de Sen [67] , en el cual la funcion entropica es definida come, 

£{u,v,e,q,p) = 27r{eiqi - f{u,v,e,p)), (6.6) 

y f{u,v,e,p) cs la dcnsidad Lagrangiana \/—det gC, cvaluada cerca del horizontc dc ge- 
ometria. Todos estos parametros pueden ser determinados extremando la funcion entropica. 

Las ecuaciones (6.7) representan el conjunto de ecuaciones de movimiento del sistema 
escritas cerca del horizonte de geometria del fondo extremal (6.4). Asi, la entropia del AN en 
el limite extremal se obtiene resolviendo el sistema de ecuaciones (6.7) y sustituyendo estos 
parametros en la funcion entropica. Esto demuestra que Sbh/'^t^ puede ser considerada 
como la transformada de Legendre de la funcion /(l^, 1^, , p^), con respecto a las variables 

6.2.1 Soluciones analiticas para invariantes de 2'^° orden en d = 4 

Los calculos para contribuciones de segundo orden a la entropi'a de ANs extremales cor- 
responden a considerar solamente los invariantes de segundo orden. Sobre la base del 
principio de correspondencia, la solucion de Gauss-Bonnet(GB) debe estar contenida en 
estos resultados [88], asf como la solucion de Reissner-Nordstrom (RN). Por tanto, en este 
caso tomamos Oj = 6j = Cj = para i > 3 en (6.3), entonces la funcion f{y, e,p) sera, 

f{vi,V2,e,p) = j y^[R + K-^+a2R^ + h2R^ded(t), (6.8) 

mientras la funcion entropica tomara la forma, 

= {2,eG, + 2 (2- ^1 - I) (-|L _ _ 1) 




y el sistema de ecuaciones de movimiento basandonos en (6.7) sen'a, 

lSA_i + 25i-2!lf + |L-^l^ + ii + ^4 = 0, (6.11) 

2 vi V2^ 8vi 8v2-^ 4^1 4i;2^ 

A ^ ^ „ V2 Q2 „ ^2 , ^^2 b2 b2 t;2 e2 p2 

-V2 — - 1 + 2 — 5- - 2 h - — 2 - h - — 2 + 1 — = 0- (6-12) 

2 vi^ V2 8vi^ 8v2 Avi^ Av2 
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Note que cuando rcsolvemos este sistema de ecuaciones, todas estas soluciones pueden ser 
escritas en terminos de la funcion V2 como sigue, 

f = ^/8 {Av2 + 2)v2 - 2(16 02 + 62) + 2)v2A - 4.p^ , 

y la entropia para un AN extremal y estatico (teniendo en cuenta invariantes de segundo 
orden) sera, 

....(._a^-±MA)-_|^. (e.M, 

La solucion de GB que fue encontrada por Morales y Samtleben en [88] puede ser 
obtenida sustituyendo 62 = —8a;, y 02 = a/2 en (6.14), asi, 

Sgb = {v2 + H-^, (6-15) 

donde el parametro a es la constante de acoplamiento de GB. De igual manera, si tomamos 
0^2 = ^'2 = en (6.14), entonces la entropia se convierte en la conocida entropia del AN de 
RN en el li'mite extremal y con carga magnetica p, 

SBH\a2=b2=0 = = 'S'rjv • (6.16) 

Las variables vi,e y q ser an entonces, 

V2 ^2 {kv2 + 2)V2 - p2 72 {AV2 + 2)^2 - 

"^ = ^;^aTi' ^ = 4G^ ' ' = — ^^aTi — ' ^'-''^ 

donde directamente de (6.14) se puede ver que la constante cosmologica por si misma no 
modifica la ley del area de Bekestein-Hawking. Por tanto ella necesita estar acompaiiada 
de terminos de derivacion superior, de al menos segundo orden, para poder modificarla y 

obtener un resultado nuevo. O sea la constante cosmologica A solo cambia la geometria 
de la garganta como se ve en (6.13) y (6.17). Sin embargo, las constantcs 02 y 62 en 
(6.14) representan las desviaciones de dicha ley. En la siguiente seccion se obtendran solu- 
ciones aproximadas para gravedad de orden superior, teniendo en cuenta las contribuciones 
restantes R^,B!^y R^. 

6.2.2 Soluciones aproximadas para el conjunto complete en d = 4 

Debido a la no linealidad de las ecuaciones de Einstein-Maxwell es muy diffcil encontrar 
soluciones exactas, mas aun si incluimos terminos de alta derivacion. En muchos casos 
se emplean metodos aproximados o soluciones numericas. Por tanto, si consideramos el 
conjunto completo de invariantes de Riemann entonces no podremos encontrar soluciones 
expKcitas del sistema de ecuaciones de movimiento (6.7). Nosotros resolveremos cstc prob- 
lema introduciendo un parametro w para hacer una expansion alrcdcdor dc cl. En cfccto 
en (6.3), w puede sicmprc ser extrafdo de las constantcs de acoplamiento ai,bi,Ci,ei por 
re-escalamiento, teniendo en cuenta que debe ser extrafdo con el orden apropiado. La 
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solucion de RN (caso Cmv = 0) es bien conocida analiticamente, y las primeras constantes 
de acoplamiento para invariantes de orden superior son 02 y 62- Consecuentemente para 
este caso: a2, 62 ^ w a2,w b2, ■■■ y en general, tendremos que, 

ai,bm,Cn,ep =^ vJ-~^auw^~^bm,w^~^Cn,w^~^ep, l,m,n,p = 2,3,4,5, ... (6.18) 

El parametro de expansion w puede ser considerado proportional al inverso de la carga 
del AN, la cual es cero a primer orden. Esta es la forma en que w debe aparecer en 
la funcion entropica 6. Ademas, como estamos interesados en soluciones aproximadas, 
entonces debemos construir series de expansion de las funciones e,q y vi. En general, 
cualquiera de estas funciones puede ser expandida alrededor de w, 




(6.19) 



Esta expansion nos permite escribir, 

e 2± eo + wAei + w'^Ae2 + u;^Ae3 + 
g ~ go + wAqi + w'^Aq2 + w^Aq^ + 

vi c:^ vio + wAvii + w'^Avi2 + w^Avis + (6.20) 
S c^So + wASi + v?AS2 + w^AS-i + 

donde en (6.20) hemos tornado A% = {l/k\){d^h/ dw^)^^^. Asi, w fija el nivel de aprox- 
imacion, en el caso que w = Q tenemos soluciones de orden cero que seran de tipo RN 
(6.16)-(6.17). Finalmente incluyendo el conjunto completo de invariantes (o sea en (6.3) 
tomando todos los terminos) , o contribuciones de quinto orden para la entropi'a de un AN 
extreme, y despues de resolver el conjunto de ecuaciones iteradas que aparecen al extremar 
la funcion entropica, en este trabajo se obtienen la siguientes expresiones aproximadas para 
la entropi'a Sbh-, 

Sbh = StN + + S'iljw + s'^IjW^ + s'^IjW^) w + Oiw'), (6.21) 

donde, 

S'bli = - ((16a2 + 62) jV2 + I) ^, (6.22) 
4% = f 12A2i;2 03 + dv2A'' + 2A + -)b3) ^, (6.23) 

V 4 V2 J (Jd 



■.(4) *3 (v2A + 2fbi 2A (U2A + 2) (V2A + 1) C4 

41 = -< 32^2 A^ 04 + ^ ^ / — - + — — — + (16a2 + 62) X 

I 0V2 V2 



A {V2A + 2)2 ba 

16^2 



+ V2a3A-' ] } — ^ 

Lr4 



TT 



(6.24) 
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^BH = \ 80^2 A^a5 + {5v2A^ + 16A + — )65 + 96^2 A^af + 4(— + 3v2A^ + 8A)a363 + (6.25) 

I V2 V2 

, (3i;2A + 4)(^2A + 2)\ 2 , na ^u^(ak2 , {v2k + 2fc^ , (i;2A + 2)% ^^ | tt A^ 
+ + (''"^ + ('"^ "^"^ + 4^;^ + 64A^.3 J I 

En caso que estudiemos una teon'a efectiva donde solo invariantcs dc tercer orden son 
necesarios, entonces dcbcmos toniar en (6.3), (6.35) las constantes de acoplamiento con los 
valores aj = bj = cj = ej = para j > 4, y la solucion para la entropia del AN extreme, 
sera simplemente: Sbh = Srn + {S^bh + '^bh'^)''^ + 0{w^). Esto mismo se puede hacer 
para cuarto orden. Esta aproximacion es mas que suficiente, pues se puede ver facilmente 
que la solucion anali'tica de segundo orden (6.14) es exactamente reproducida en (6.22). 
Note que cada contribucion cntropica de invariantes de i— esimo grado cs ctiquctada como 
S^*) , entonces el super-mdice (i) solamente indica el grado del invariante que produce esto 
y no el orden de aproximacion. 

En el conjunto de soluciones aproximadas los terminos no lineales con respecto a la con- 
stante de acoplamiento aparecen en S^^lj y S^^lj con contribuciones como 0203, 0263, 6203, ?'2&3 
y a|, 0363, 0204, 0204, 0264, 6204, b2C4, h2bi respectivamente. 

6.2.3 Soluciones analiticas para invariantes de 2'^° orden en d = 5 

En este caso, la funcion f{v,e,'p) toma la misma forma que (6.8), pero la integracion esta 
vez se realiza sobre S^, entonces la funcion entropica es dada por, 

7r2(3t;i -?;2)^a2 3 tt^ {2 vi + V2f b2 



£{vi,V2, e,q,p) = 27rge - 



G5 Vl s/v^ 40 G5 vx ^Jni 

n2{2vxV2{h.vx - 2) + e^V2 + 12 vi^) tt^ 



(6.26) 



8G5 Vl 

y el sistema de ecuaciones dc movimicnto (6.7) ccrca del horizontc scran, 

Q _ 7r(3t;i ^V2) (3t;i - ^2)^2 ^ 37r (2^i - V2) {2vi + ^2)^2 _ ,^ 27) 

\fv^{^-V2v\ - — 2kv\'V2 + e}v2) 
Q ^ 37r {v\ + V2) (3t>i - t;2)a2 ^ 37r {^vx + V2){2vx -3^2)^2 _ 

^JV2V\^ 4S)^/V2V\^ 

_ 37r ^2 {2vl{2 + A^2) + ^2(e^ - ^Vl)) 

3 

El sistema puede ser resuelto expli'citamente, y entonces las funciones e, y g pueden ser 
escritas en terminos de V2 como, 

= ^2 (5t;2 + 20a2 - ?>2) „„x 
60a2 + 2 62 + 5?;2(4 + i;2A)' ^ ' ' 

^ (6.31) 



8 {hv2 + 2002 - &2) G5 ' 20u2 + 60a2 + 262 + hv\k ' 
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donde, 



/ = V5 {v2^K + 2002 + 2 {iQvl - 600262 - 40a2uiA - 2hl - 362^'iA - IO62W2) 



Por consiguiente, la entropia de un AN estatico y extreme en d = 5, teniendo en cuenta 
terminos de derivacion superior de segundo grado, tiene la forma, 

1 

^ Ti'^V2^ (4002 + 352) A Stt^^I {-V2^ + IOQ262 + 262^^2) , . 

4(-5?;2-20o2 + 62)G5 2 (-5^;2 - 2O02 + 62) G5 ' ^' ' 



De forma analoga, la solucion extremal de Gauss-Bonnet mostrada en [88] puede ser 

I' 



obtenida con las sustituciones, 62 = —8a, y 02 = la, en la solucion principal (6.32), 



asi, 

2 i 

SGB = {v2^Vla)—^. (6.33) 

Esto no pasa para una contribucion generica tipo GB. En general, el establecimiento de 

una rclacion entre las soluciones de GB y las solucioncs dc segundo grado requieren de un 
vinculo (ver seccion 6.2.5). Por otro lado, si 02 = 62 = en (6.32), la conocida solucion de 
RN extremal tambien es obtenida, 

3 

^2y 2 

SBH\ai=hi=Q = = ^RN- (6-34) 

Asf como acurrio en (6.14), la constante cosmologica no cambia la entropfa del AN por 
si misma: Ella podria necesitar de mayores terminos de derivacion superior. Una prueba 
de esto sigue de tomar A = en (6.32), donde se ve claramente que las contribuciones de 
02 y 62 no desaparecen. Estas contribuciones tambien producen una marcada desviacion 
de la ley del area. Note que, contrariamente al caso de d = 4, no solo la constante cos- 
mologica modifica la geometn'a de la garganta del AN (ver (6.30)), sino que las constantes 
de acoplamiento asociadas a los invariantes de segundo orden, tambien tienen su efecto 
sobre la topologfa de la garganta. 

6.2.4 Soluciones aproximadas para el conjunto completo en d = 5 

Como hicimos en la seccion 6.2.2, podemos construir aproximadamente soluciones para 
cinco dimensiones. Asf, considerando un AN extremo en el cual el conjunto completo de 
invariantes de Riemann (o sea en (6.3) se deben tomar las constantes de acoplamiento 
aj , bj , Cj , Cj 7^ 0) a diferentes niveles de aproximacion tiene la siguiente forma, 

Sbh = 4f + (S^BH + S'iljw + S'iljw'')w + Oiw^), (6.35) 



donde. 



S^BH = - (2 (1 + At;2) 02 + ^ (6 + A^;2) 62) ^ , (6.36) 
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S^BH = i^6{l + Av2fa3 + ^{6 + Av2)iSAv2 + 8)b3-^{6 + Av2fc3+ (6.37) 
+ 8(1 + Av2)al + (-9 + V2A)^ - ^(^^A + 6)6i 1 , 

Sbh = -^G{l + Av2fai + ^{& + Av2fhi + ^{& + Av2){2Av2 + 7){l + Av2)ci- 

- ^{4.Av2 + ^){<o + Av2fe^-?,2{l + V2A)al + Ull + V2A)b2al+ (6.38) 

+ (9 + Av2A)^ - 8(10 + t;2A)(l + V2Afa2a3 - (270 + 248i;2A + ?,vlA^ + bQvlA^)^ + 
25 5 

+ ]|q(^2A + &)(vlA^ + 16t;2A + 30)a2C3 - ^{v2A + 6)6i - ^(1 + ^^2A)(t;2A + 5) x 



X (^;2A - 4)0362 - ^(^2A + Q){2,vlA^ + 8^;2A - 20)6263 + ^(^sA + 6)^6203 



1^ 



J G5V2 

Es correcto decir que si la ecuacion (6.32) se expande en series de segundo orden de 
V2^, entonces logramos reproducir la ecuacion (6.36). Si nosotros tomamos 02 = |q; y 
62 = —8a en (6.36) se obtiene la solucion de GB (6.33). En este caso, las contribuciones 
no-lineales de las constantcs dc acoplamicnto aparcccn cn S^^\j y S^^\j como combinaciones 
de las constantes de acoplamiento de pequeiio orden como 02, 0262, 62, y a2,.--etc. Note que 
la ausencia de terminos 03,0262 y 63 en (6.23) en contraste con (6.37) es simplemente un 
resultado del calculo. Lo mismo pasa en (6.24) y (6.38) con los terminos a^, 620^, a26|, 63. 

6.2.5 Invariantes de Riemann y gravedad generica de Gauss-Bonnet 

En esta seccion examinaremos la teorfa de Einstein-Maxwell con una modificacion a la cual 
Uamaremos Gauss Bonnet generica (o sea con tres constantes de acoplamiento diferentes) . 
Esto es algo que debe hacerse con cuidado pues se sabe que en general una teoria asi no 
es renormalizable, pero en nuestro caso solo lo hacemos a modo de generalizar nuestro 
resultado. Ademas de todas formas se sabe que las teorias de Gauss-Bonnet, aunque 
superan el problema de la renormalizacion, tienen problemas con la existcncias de parti'culas 
fantasmas. Asf la accion para una teorfa de GB generica en d-dimensiones puede ser escrita 
como. 



(6.39) 



de este modo, cuando se aplica el mecanismo de la funcion entropica, la cntropfa de un AN 
extremo para d = 4 toma la forma. 



,(rf=4) _ 27r (t;2^A^ + 2i;2A + 2) xi ^ 47r (2+ z;2^A^ + 2 1;2 A) X2 



G4 {V2 A + 1) 

2tT V2'^A^XZ T[V2 
G4 {V2 A + 1) G4 ' 



G4 {V2 A + 1) 



(6.40) 
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Las soluciones para las funciones vi, e, y q son, 



V2 fo fo rc: A-,\ 

fo = I -^XivW - 16X1^^2 A + lQx2vl^^ + 32x2t'2A - Sxa^iA' 



1/2 



16 X3^^2 A + 2t;2^ A + Av2 - 



entonces si Xi = X2 = Xs = X) sn (6.40), se obtienen las soluciones (6.15) con (x = ct). 
Sin embargo, si queremos encontrar un resultado general para GB (6.40) a partir de las 
soluciones analiticas (6.14), (o viceversa) entonces los siguientes vmculos deben exigirse 
sobre ambos resultados, 

362 62 
XI - 2x2 + X3 = 0, XI = 2X2 + ^ + a2, X3 = 02 - ■ (6.42) 

Igualmente para d = 5 la entropia resultante viene dada por, 
2 

^Gb'^ = r^^^(-8^^i + Xi(48x2V^ - 2^{4v2 + 16x3 + ^v^)) - 

[V2 + 4X3 jCrS 

5 

-64V^xi + 4 V^( 20x3 + + 4i;2)X2 - ^(^^sA - 1)0^.43) 
para las funciones vi, e y q, dadas por, 

VI = , ;2(.2 + 4x3) 

4^2 + Av^ - 16x2 + 12X3 + 4X1 ^ ' 

q = q{v2,A,Xi,X2,X3), e = e(t;2, A, xi, X2, Xa) • (6.45) 

Tambien, la relacion entre las entropi'as (6.43) y (6.32) vienen dadas para cuando se cumplan 
los siguientes vmculos, 

Xi-2X2 + X3 = 0, a2 = ^(2x2-Xi), 62 = 8(xi - 2x2) • (6.46) 

Si entonces tomamos Xi = X2 = X3 = Xi ^n (6.43), se obtienc la ccuacion (6.33). 

Estos resultados son interesantes porque ellos se corresponden con casos donde los 
invariantes de Riemann no pueden reproducir los resultados de la teon'a generica de GB. 
Esto nos sugiere que a nuestra teoria siempre le faltaran invariantes por adicionar. Aunque, 
esta situacion es irrclcvante si se tiene en cuenta que es bien conocido que la teon'a de GB 
sobre una gravedad AdS es en general inconsistente debido al principio variacional [98] y 
al problema de la regularizacion que presentan. Una discusion extensiva de este problema 
puede verse en [99], [100], [101]. 
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6.3 Conclusiones y perspectivas del capitulo 

Nosotros hemos calculado la entropia de ANs extremes en el caso ded = 4:yd = 5 
dimensiones, teniendo en cuenta terminos de derivacion superior construidos a partir del 
conjunto completo de invariantes de Riemann, en una forma, que hasta donde sabemos, no 
se ha realizado con anterioridad. 

Aunque el caso de GB fue estudiado con anterioridad en [88], en nuestro trabajo, 
hemos cncontrado invariantes de segundo orden que han gcncralizado los rcsultados del 
caso GB cncontrando como caso particular el caso dc RN. Tambien se han obtenido los 
terminos principalcs para las aproximaciones en el caso de terminos de derivacion superior. 
Por consiguiente, el conjunto de soluciones obtenido muestra un ejemplo concreto de que 
el formalismo de Sen funciona correctamente y es un proceso menos complicado que el uso 
de la formula de Wald [94]- [97]. Para poder aplicar la formula de Wald se debe conocer 
la solucion exacta de las ecuaciones de campo del problema y esto solo se tiene para un 
mimero muy reducido de cases. 
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7. CONCLUSIONES Y PERSPECTIVAS 



Conclusiones 

En esta tesis nos plantcamos dos objetivos inspirados en el estudio de la Fisica de Objetos 
Compactos en presencia de campos magneticos intensos. El primer objetivo, cl estudio del 
colapso de un gas denso autogravitante de materia fermionica: electrones y neutrones sobre 
la base de un modelo construido para describir, al menos cualitativamente, el interior de 
una EB o una EN. Bajo las condiciones asumidas en el modelo, arribamos a las siguientes 
conclusiones: 

1. El gas de electrones magnetizado dependiendo de las condiciones iniciales del prob- 
lema colapsa a una singularidad tipo "punto" o "cigarro" . 

2. Para el gas de electrones magnetizado siempre existe un valor de campo magnetico 
para el cual el colapso tipo "cigarro" se extiende en la direccion del campo magnetico. 

3. Un gas dc neutrones siempre colapsa a una singularidad tipo "punto"' . 

4. El que el gas de electrones pueda tener una singularidad tipo "cigarro", y no asi para 
el gas de neutrones, se debe a que el campo magnetico se acopla de manera mas 
fuerte (mediante su carga electrica) con los electrones, a diferencia de los neutrones, 
ya que ellos no poseen carga electrica y se acoplan al campo magnetico mediante su 
MMA. 

5. Hemos introducido perturbaciones enriqueciendo la dinamica del espacio-tiempo Bianchi 
demostrando que este espacio-tiempo puede aproximarse mas a situaciones reales. 

Del segundo objetivo podemos destacar a manera de conclusiones lo siguiente: 

1. El metodo de Sen nos permitio encontrar nuevos resultados para la entropi'a de ANs 
extremos en 4 y 5— dimensiones con topologia AdS y teniendo en cuenta terminos 
de derivacion superior. Se mostro asi que estos resultados generalizan la formula de 
Area encontrada por Bekenstein-Hawking (S = A/ AG) para la TGR. 

2. El metodo de Sen nos demostro que es una herramienta realmente potente y funciona 
bien para teorfas con terminos de derivacion superior. Calculos que anteriormente 
eran imposibles de efectuar, ahora son posibles de realizarlos. 
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Perspectivas 



Como perspectivas podemos destacar varias relacionadas con los dos objetivos. 

1. El tenia de estudiar la evolucion dinamica de un gas degencrado y magnctizado 
autogravitante de fermiones, podria estudiarse en condiciones mas realistas. Se podn'a 
enriquecer la dinamica del espacio-tiempo, ya sea introduciendo perturbaciones como 
se intenta hacer en el capitulo 4 o cambiando a otros Bianchis de ordenes superiores 
en donde se presentan condiciones de 4-velocidad inclinada con vorticidad no nula, 
asi como curvatura positiva que es la que naturalmente posee un objeto compacto. 

El hecho dc incorporar fuentes de materia de distinto tipo, o sea materia mezclada, 
clcctroncs, protones, neutroncs, tambien incrementarfa la validez del modelo. Incor- 
porando las condiciones de neutralidad de carga y equilibrio beta que estan presentes 
en el interior de una EB y de una EN, ademas entender si en esas condiciones prevale- 
cen la singularidad de tipo punto o si las de tipo cigarro aparecen. 

Examinar la posibilidad de incorporar la aproximacion de arbol en el potencial ter- 
modinamico y las de un lazo para las ecuaciones de Maxwell. 

Tambien podn'a generalizarse este estudio al contexto cosmologico anadiendo los efec- 
tos de la temperatura. Ello podn'a describir la dinamica de la nucleosmtesis en pres- 
encia de campos magneticos. 

2. Para el segundo tema, se podn'a incorporar hasta orden cinco los invariantes rela- 
cionados con la teoria de calibracion o sea el mismo orden maximo que posccn los 
invariantes de Riemann y encontrar alH invariantes independientes, despucs con cl- 
los calcular la entropia de ANs extremes con topologi'a AdS de forma aproximada. 
Tambien, se podrfan generalizar estos resultados a 6,7,8,...n-dimensiones. 
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Anexo A Convenios y Notaciones 

Constantes fi'sicas mas usadas, 

Velocidad de la luz c = 2.998 x 10*^ms"\ 

Carga electrica del electron e = 1.60 x lO^^'^ C, 

Constante de Dirac h = 1.054 x 10^'^'' J s, 

Masa en rcposo del electron rrie = 9.109 x lO^'^"*^ kg, 

Constante de gravitacion G = 6.6742 x 10^^"*^ N m'^ kg~'^, 

SttG 

Acoplamiento gravitatorio k = — ^ = 2.075 x 10~"^^ m J~ 

Masa del Sol M© = 1.99 x 10^° kg, 
Radio del Sol Rq = 6.96 x 10^ km, 
Temperatura del Sol T© = 15 x 10^ K, 
Campo Magnetico Bq = 1G, 
Densidad del Sol p© = 1410 kg m~^. 

Matrices de Pauli, 



^2 = -z , = I 1 I . (A.l) 





Matrices gamma de Dirac, 




7^' = I I . (A.2) 

Si se usa el convenio h = c = 1, en este sistema, 

[longitud] = [tiempo] = [masa]~^ = [energia]"^ , (-^-3) 

se deduce que si se escribe una m sola, entonces se entendera como energia en reposo mc^, 
asi como su inverso se entendera, como la longitud de onda Compton {h/mc), 

me = 9.109 X 10-^^% = 0.511 MeV = (3.862 x 10-^^cm)-^ (A.4) 

En el contexto de la teoria cuantica de campos, se introduce una derivada covariante de la 
forma, 

D^ = d^ + leA^ , (A.5) 

donde es el 4-vector potencial del campo electromagnetico y se toma el convenio para 

la metrica del modo [H ]. Asi, = (x'^,x), y = g^u^'^ las que para el caso de 

una metrica plana de Lorentz, toman la forma = {x^, — x). Asi mismo, para el operador 
derivada se define como es usual, 
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Sin embargo, si estamos en la TGR, generalmente tomamos el convenio [- H — I — 1-] y 
la relacion cntrc las constantes Gjc = 1. Usamos indices latinos para el espacio-tiempo 
{a,b,c,d, ...i,j,k, ... = 1,2,3,4) e indices griegos para el espacio 3-dimensional (a, ^,7,... = 
1,2,3), con la excepcion del Capitulo-4 donde invertimos esta seleccion. 

Tambien, empleamos la notacion de coma "," para la derivada normal, 

T"^ u = ^TT- (A.7) 

y punto y coma (;) o V(antes del ente que se desea derivar) para la derivada covariantc. 
Ejemplo para un tensor T"^ dos veces contravariante, 

.J = ^ + r;.„r"*^ + r^^r"^ = VJ•r^ (a.s) 

donde T^.^ son los indices de Christoffel. 
Ademas se recuerda el conjunto de definiciones: 

• "2, b2,..6q Componentes de un tensor "p" veces contravariante y "g" veces 
covariante. 

• u"' Componentes de un vector contravariante en el espacio-tiempo. 

• u Componente de un vector del espacio 3-dimensional. 

• 1-forma diferencial. 

• Cj base covariante. 

En general un tensor cualquiera se puede escribir rigurosamente de la forma: 

T = r«i«2,-,ap ^^ j^^ j^y^ w''^ ® ...It;^" ® Ca, ® Ca,... ® Ca,, (A.9) 

donde las bases deben cumplir con las siguientes relaciones: 

<ei,ej > = Qij, (A. 10) 

< w\w^ > = g'^, (A.ll) 
<w^,ei>= 6^i, (A.12) 

donde <, > denota el producto escalar. 
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Anexo B Diagramas del Agujero Negro de Schwarzschild 




Figura 6: Sc muestra el diagrama del AN de Schwarzschild en las coordenadas de Kruskal-Szekeres. 
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Figura 7: Se muestra el diagrama del AN de Schwarzschild en las coordenadas de Kruskal-Szekeres. 
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Anexo C Fuentes de Materia 



C.l Fluido Termodinamico 

La ecuacion de estado mas general para un fluido termodinamico mono-componente puede 
ser escrita como, 

p = p(n,r), p = p{n,T), (C.l) 

donde n y T son la densidad del niimero de particulas y la temperatura. En la geometn'a 
RW este fluido cumple con la ecuacion de conservacion, 

Va(rau") = 0, => - + ^ = 0, noc R-'\ (C.2) 

n R 

y la ecuacion de equilibrio de Gibbs, 

donde S es la entropi'a por parti'cula. Comparando (C.2) y (C.3), con (1.42) obtenemos 
S = 0, de modo que el fluido es isoentropico (la entropi'a S se conserva a lo largo de u"), lo 
cual es de esperar para un fluido en equilibrio termico (bajo un proceso reversible y cuasi- 
estatico). Usando (1.42) y (C.2), la condicion de integrabilidad (C.3) permite obtener la 
ley de evolucion. 



dp/dT 



1=0, (C.4) 



T [dp/dT \^ R 

donde el mdice „ indica evaluacion a n constante. La integracion de (C.4) para cualquier 
ecuacion de estado como (C.l) permite una relacion del tipo T = T{R,n), pero n oc R~^, 
la relacion T = T{R,n) transforma a (C.l) en p = p{R), permitiendo la integracion de la 
ecuacion de Priedmann (1.43). 

C.2 Gas ideal relativista no-degenerado. 

Un importante ejemplo de un sistema termodinamico es el gas ideal relativista de una 
componente, asociado con la distribucion de Juttner, la cual es la generalizacion relativista 
de la distribucion de Maxwell-Boltzmann en la Teoria Cinetica Relativista. Una descripcion 
hidrodinamica para este tipo de gas viene dada por la ecuacion de estado, 

p = mc^nT{l3)-nk^T, p = nkj^T, (C.5) 

donde, 

aqui m es la masa de la particula, kg es la constante de Boltzmann y K2,K^ son funciones de 
Bessel modificadas de segundo tipo, de segundo y tercer orden. La ecuacion de estado (C.5) 

es un poco complicada, aunque se simplifica para dos situaciones extremas de temperatura 
y energia: el regimen ultra-relativista (UR) y el no-relativista (NR), caractcrizados por 
^ <C 1 y /3 3> 1 respectivamente. Es litil apreciar estos rangos de temperaturas para estos 
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casos li'mites. Considerando masas de orden nuclear (~ 10^^'' kg), con ~ 10"^^ J/K, 
nosotros tenemos que ^ « 1 para T « 10^^ K. De modo que podemos modelar materia 
cosmica basandonos en un gas ideal no relativista. El caso ultra-relativista podn'a ser 
valido para eras cosmicas muy calientes. Considerando el comportamiento de r(/3) para 
estos casos limites, nosotros tenemos, 

^<1, r ^ ^ + ^ + 0{P^), p^SnkgT UR, (C.6) 

P^l, T ^1 + ^ + 0{p-^), p^mc^n+^nkgT NR. (C.7) 

Identificando la ecuacion de estado para el caso del gas ideal ultra-relativista (UR), 

p = 3nkgT, p = nkgT = ^p, (C.8) 

y para el gas clasico no-relativista (NR), 

3 

p = m(?n + -nkgT, p = nkgT, (C9) 



de las ecuaciones de conservacion (C.2), y (C.3) se obtiene inmediatamente para el gas NR, 

p{R) = nokgU , (C.ll) 



p{R) = mc^no ( ^ 



2 mc^ 



(C.12) 



donde Rq, no, Tq son R, n, T evaluados en el valor fijo t = to- La ecuacion (C.12) es la 
relacion necesaria p = p{R) que transforma a (1.43) en una cuadratura integrable. 
La ecuacion de estado del gas ideal relativista (C.5) es aplicable a gases diluidos en los 
que las distancias medias entre las parti'culas (entre colisioncs) son mucho mayores que 
la longitud de onda termica. Por lo tanto, los efectos cuanticos del tipo de estadi'stica, 
pueden ser despreciados (gases no-degenerados) . En los capi'tulos 2 y 3 se tratara el caso 
de gases de fermiones, en los cuales estas suposiciones no son validas. Aunque las ecua- 
ciones de estado resultantes seran mucho mas gcncrales, el tratamiento dinamico atraves 
de las ecuaciones de Einstein para estos gases sigue la misma mctodologfa que la descrita 
anteriormente. Sin embargo, la inclusion de un campo magnetico introduce la existencia 
de una direccion privilegiada, y es por lo tanto incompatible con la metrica FLRW. Para 
este caso se utilizaran modelos Bianchi I, que son homogeneos, pero anisotropicos. 

C.3 Polvo 

Es evidente a partir del sistema de ecuaciones (C.2), que expandiendo al Ifmite i? S> Ro, la 
presion y la energfa dccaen rapidamente como 0{R~^), y el resto de la masa encrgia como 
0{R~^). Por otro lado el rango k^T/mc^ = 1/ j3 es despreciable (<C 10~^) para masas de 
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nucleones o clectrones y para un ancho rango de temperaturas asociadas con condiciones 
clasicas (T < 10^ K), entonces es evidente a partir de (C.12) que practicamente para 
toda su evolucion, el gas ideal puede ser razonablemente aproximado suponiendo p y 
p m(?n. Puentes de materia que satisfacen la ecuacion de estado, 

p = 0, p{R) oc R-^, (C.13) 

son llamadas "polvo", y reducen a (1.38) a la forma Tah = pUa^b- Aunque (C.13) puede ser 
justificado como limite de temperatura cero de un gas ideal NR, dejando a p = mcPn (o sea 
toda la densidad de materia-energia seria igual a la energia en reposo) . Esta ecuacion posee 
una convcnicnte descripcion aproximada de ciertos tipos de materia cosmica no-rclativista. 
Por ejemplo, materia fn'a no colisionada, o materia oscura fn'a no-barionica (aunque esta 
ultima relacion no es del todo segura para materia oscura). 

C.4 Radiacion 

Para el gas ideal UR, de las ecuaciones de conservacion (C.2), y (C.3) se obtiene que, 

n.nJ!^)\ r.r„(f), (C.14, 



asi que, 

p{R) = Snok^To j =MR), (C.15) 

permiten la integracion de (1.42) y (1.43). La ecuacion (C.8) puede describir aproxi- 
madamente un gas de particulas masivas (m > 0) a muy alta temperatura, pero esta es 
exactamente la ecuacion de estado para un gas de particulas no masivas, tales como fotones 
o neutrinos. Una relacion equivalente a (C.15) para fotones es la ley de Steffan-Boltzmann, 



p^ = a,T' = a,T^i^^y , (C.16) 

dondc aj cs la constante de Steffan-Boltzmann, asi que la conservacion dc la densidad 
del numcro dc fotones es uq = Tq /3kg . El termino radiacion es aplicado a un gas de 
fotones caracterizado por el espectro de energi'a de un cuerpo negro. La Radiacion Cosmica 
de Fondo o CMB es la reliquia de eras primordiales en el Universe donde dominaba la 
radiacion, es la mas importante manifestacion de este tipo de gas en el contexto cosmologico. 



*CMB viene del ingles Cosmic Microwave Background 
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Anexo D Ecuacion de Raychaudhuri 



Antes de deducir la ecuacion de Raychaudhuri definamos primero la congruencia de geodesicas: 

• Sea M una variedad y sea O C M un abierto. Una congruencia en O es una familia 
de curvas tales que a traves de cada p e O pasa precisamente una curva de esta 
familia. 

• Una congruencia se dice que es suave si el correspondiente vector del campo de los 
vectores tangentes es suave. 

Consideremos una suave congruencia de geodesicas temporaloideas . Suponemos que las 
geodesica son parametrizables por el tiempo propio r, asi que el vector del campo u", 
tangente a las Imeas dc materia, es normalizado a la unidad de la forma u"'Ua = —1. 
Entonces se puede definir cierto tensor del campo B^b, 

Bab = ^bUa , (D.l) 

que sera puramente espacial, o sea, 

BabU" = BabU^ = . (D.2) 

Consideremos una suave subfamilia 7s (r) uni-parametrica de geodesicas en la congruencia. 
Sea 1]°- el vector de desviacion ortogonal de 70 para esta subvariedad. Entonces r] representa 
el desplazamiento infinitesimal espacial desde 70 hasta una geodesica cercana, 

^urf = 0, (D.3) 

donde, la anterior derivada de Lie quiere decir que, 

txVbJy" = ri'^Vbu" = BiriK (D.4) 

Entonces B puede ser interpretado ffsicamente como quien mide la desviacion de 77" al ser 
transportado paralelamente. Un observador sobre las geodesicas de la familia 70 podra 
encontrar geodesicas cercanas a el rodeandolas estiradas o rotadas por el mapa lineal B"^^. 
Continuemos ahora con el siguiente conjunto de definiciones necesarias: 
La metrica espacial puede ser definida como: 

Kb = gab + UaUb , (D.5) 

donde el tensor de proyeccion ortogonal sobre la hiper-superficie ortogonal a es dado 
por, 

h% = g''''Kb. (D.6) 

La expansion, se define como: 

e = B-'Kb . (D.7) 
El tensor de deformaciones (o cizallamiento) se define segiin: 

o- = B^ab) - oKb , (D-8) 
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y el tensor de vorticidad como: 

De esta forma Bab puede ser descompuesto como: 

Q 

Bab = -^hab + CTab + ^ab- (D.IO) 

Todos ellos son piiramente espaciales ya que, aabU^ = UabU^ = 0. Ahora los respectivos 
significados fisicos de estas magnitudes son: 

• 9: Mide el promedio de expansion de geodesicas circundantes e infinitesimamente 
cercanas. 

• LOab- Mide la rotacion de una geodesica circundante e infinitesimamente cercana re- 
specto a la otra. 

• aab- Mide la dcformacion. O sea, sea una esfera inicial cn cl cspacio tangente la cual 
es transportada segiin Lie (transporte paralelo) a lo largo dc y se distorsiona en 
un elipsoide con los ejes principales en la direccion de los eigenvectores de a" 

De la ecuacion de la desviacion de la geodesica es facil derivar lo siguiente: 

U^VcBab = u'^VcVbUa = U%VbVcUa + i?e6a«d) (D.ll) 
= -B%Ba^ + RjaU^Ud. (D.12) 

Tomando la traza, obtenemos, 

fI0 f)^ 

u^Vce = — = --- aaba"'' + a;„ba;«^ - RcdU^u'^. (D.13) 
dT 6 

La ecuacion (D.13) es la conocida ecuacion de Raychaudhuri y es la piedra angular para 
demostrar los teoremas de singularidad. Usando las ecuaciones de Einstein's y la normal- 
izacion de en el ultimo termino de la ecuacion de Raychaudhuri (D.13) tenemos: 

Rabu'^u'' = 87r[r„6 - ^gabKu'' = 87r[Tabu''u^ + |] , (D.14) 

donde TabU°'u^ representa la densidad de energia de la materia medida por un observador 
cuya 4-velocidad es (temporaloidea). De la expresion (D.14) se pueden definir las 
siguientes desigualdades que constituyen las condiciones de energia para la materia: 

(Condicion de energia debil) Tabu"-u'' > 0, (D.15) 

T 

(Condicion dc cncrgi'a fuerte) Tabu"'u^ + — > 0, (D.16) 
(Condicion de energia dominante) — Tj^v!' = vector temporaloideo o nulo. (D.17) 

• Esta ultima condicion (D.17) debe ser un vector temporaloideo o nulo, el representa 
la densidad de 4-corriente del tensor energi'a-momentum de la materia. Medida por 
un observador con 4-velocidad Puede ser interpretado como la velocidad del flujo 
de energia de la materia con velocidad menor que c. 
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• Si V^'Tab = entonces la condicion de energia dominate sc satisface. Es decir, se 
desvanece sobre un conjunto cerrado acronal S, entonces tambien se desvanece en 

Note que la condicion de energia dominante implica la condicion de energia debil, pero en 
otro caso todas estas condiciones seran matematicamente independientes. 



Esto constituye un lema, o sea el Lema 4.3.1, enunciado por Hawking-Ellis en 1973. 



Anexo E Significado del tiempo adimensional (r) 

El tiempo adimensional (r) se define a partir de la ecuacion (3.51): 



Hot, 



(E.l) 



A. - J_A 

dr Hq dt 

en la cual podemos notar que el signo de r depende del signo de f y del signo de Hq. Para 
comprender el sentido fi'sico del tiempo adimensional r observemos que de (3.42), (3.51) y 
(3.52) obtenemos: 



V = ABC^ ^ = e^foHdt^ 



y despejando 7i de (2.42) tenemos: 



n = ±=;\/sru + si + si + S2S3 . 



Si fijamos t ^ 0, en (E.3) tenemos la siguiente interpretacion fi'sica: 



(E.2) 
(E.3) 

(E.4) 



7 = 3 J Hdt 




= Si/>0=^y>Vo (expansion) , 
= Si7<0^F<Vb (colapso). 



(E.5) 
(E.6) 



Las posibles combinaciones de signo se dan en el cuadro (1), en el cual aparecen 
resaltados los convenios que hemos escogido en esta tesis. 



Casos 


Expansion 


Volumen elemental 


T 


H>0 


n>Qy Ho>0 


Expansion 
V>Vq 


r > 


UKOyHoKQ 


r < 


H <0 


HkOj Ho>(} 


Colapso 

V <Vq 


r > 


n>Oy Ho<0 


r < 


H = 


n = o 


V = Vq 





Tabla 1: Posibles combinaciones de signos entre H, Hq y t para t > 0. 

Notemos que el sistema comienza a evolucionar a partir de t = r = 0, para r > el 
sistema colapsa y para r < el sistema tiende hacia un estado diluido como se puede ver 
en la Figura 5. 
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Anexo F Espacios de fase 




4 



Figura 8: Traycctorias cn la scccion del cspacio dc fasc {S^',l3,ii) para cuatro condicioncs inicialcs 
diferentes de un gas magnetizado y auto-gravitante de electrones degenerados. El piano sombreado 
/X = 1, esta acotado por — 4<S'^<4yO</3< 0.2. El atractor "estable" es el punto marcado 
por a, y las curvas b,c,dy e son soluciones numericas sobre la seccion 3-dimensional del espacio de 
fase. Todas las traycctorias caen desde r = tcndicndo al atractor estable a. Similarmente, desde 
T = (esfera vacia) todas escapan tendiendo a la singularidad. 
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Figura 9: El grafico muestra varies caminos en una subseccion S-dimensional (S'3,/3, /i) del es- 
pacio de fase para un gas magnetizado y auto-gravitante de neutrones degenerados. Los puntos 
oscuros representan las condiciones iniciales a "r = 0" . El punto "a" representa un "atractor" con 
coordenadas {S3 = 0, f3 = 0, fi — 1,71 — 0). 
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Anexo G Invariantes de Riemann para ^4^52 x 



Invariantes 


Definiciones 


Soluciones 
d = 4 


Soluciones 
d = 5 


R 


9°''^ d"^ Rahcd 


271 


273 


ri = R2 


1 Q b Q a 


4/2 


10 '4 


r2 = R3 


1 Q b Q CQ a 







rs = Ra 


1 Q b Q c Q d ca 
16 '-'a '-'b '-'c '^d 


64 '2 


^7i 


ine(u;i) = 'Kt{W2) 


I ri C abed 


6 '1 




^Z{W2) = ^t{W^) 


1 1^ cdri e//^ ab 
lE^ab '^cd ^ef 


36 '1 




IHe(mi) = lHe(M3) 


1 cab Qcd/-i 

'^acdb 


T5-7i7i 




lHe(m2) = d\z{M^) 




36 '1 '2 




r?i3 = M+ 




1 , 2, 2 
36 'i '2 




«Hc(m5) = ^He(M5) 


±S'"^S''-l'C"-3'^''{CacdbCgefh + Cl^db^lefh) 


108 '1 '2 




donde : 


71 = ("Wl - V2)/viV2, 72 = ('(-'1 + V2)/viV2, 
73 = (3t'l - V2)/viV2, 74 = (2^1 + V2)/viV2 







Tabla 2: Cuadro de definiciones del conjunto complcto dc invariantes de Riemann y sus resTiltados 
para la geometria tipo AdS2 x S''-~^ cerca del horizonte de un AN end = 4y(i = 5 dimensiones 
(el resto de los invariantes son nulos). Todas las definiciones son basadas sobre el tensor de Ricci 
de traza nula Sab, el tensor de Weyl Cabcd y (?1 tc^nsor dc Riemann Rabcd- Aunque, estos tensores 
tambien pucdcn scr dcfinidos sobrc una base cspinorial [30]. Aquf, "-" significa que no existe 
una definicion para los invariantes complejos, para el caso de cinco dimensiones. Nosotros usamos 
aquf letras miniisculas para los invariantes de Carminati, pero en el Lagrangeano usamos letras 
mayiisculas. El sfmbolo 9le significa la parte real del invariante complejo, y para re-marcar el 
grado del tensor (en lugar de los subfndices de Carminati) hemos escrito el grado del tensor como 
un subfndice (ej. ^Rt{W2) significa la parte real del invariante complejo de Weyl de segundo orden 

Wl). 
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d 


Grado 


Invariantes 


Equivalencias 


d = 4: 


-yst 


i? = 27i 


71 


ond 


= 472, i?2 = hi DKt{W2) = hi 


7?, 7| 


grd 


B? = 87?, RxR2 = hill R X MW2) = hi 
9te(VFs) = w7?, JHcfMs) = w7i7l 


7i, 7172 


4th 


= I674, M{W2) = ^ i?2 = 7f72', 

^2 = ^72', X lHe(H/2) = §71', 
i22 X lRe(l^2) = l^lhl = M+ = ^72^1, 
X 9^e(M/3) = X lRe(M3) = ^^4^, 

R* = mTs' 


7^ 7I, ihl 




= 327f , R^xR2 = 27372, i?3 X 9le(W2) = §71, 
i?2x$Ke(W3) = -|, i?ixi?=^, 
i?2 X ^e(M3) = ^172 ^ ^ X 5He2(W2) = ll , 

RXR2X fnciW2) = 4?, 

i?2 X 9le(W^3) = ~fj]^ , i?2 X fHe(M3) - , 
lHe(W2) X 9^e(M3) = fHe(W^2) x lHe(1^3) = 

4 3 2 

K X K4 = -32^, -TL X inC(M4j = it X 71^4 = 


iliiitihl 


d = 5 


\§L 


i? = 273 


73 




d2 ^^,2 E) 3 ^,2 
^ = 473, R2 = io74 


^,2 ^,2 

73, 74 




i?3 = 87I, Rx R2 = I7374', i?3 = Tio7l 


73> 7374 > 74 


4th 


= l673^ ^2 = mil r:'x R2 = hhl 

RxR3 = ^7374^ R4 = #7l 


4 4 2 2 
73, 74> 7374> 

7374 




donc/e : 71 = {vi — V2)/viV2, 72 = {vi + V2)/viV2, 

73 = (3i;i — V2)/viV2, 74 = (2^1 + V2)/viV2 





Tabla 3: Conjunto completo de invariantes no nulos, ellos han sido organizados por grados para 
d = 4 y d = 5, basandonos en la geometria AdS2 x S'^~'^. 
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